MIT/WHOI  2005-13 


Massachusetts  Institute  of  Technology 
Woods  Hole  Oceanographic  Institution 


OfTtC^ 


Joint  Program 
in  Oceanography/ 
Applied  Ocean  Science 
and  Engineering 


1930 


DOCTORAL  DISSERTATION 

Adult  Demography  and  Larval  Processes  in  Coastal 
Benthic  Populations:  Intertidal  Barnacles  in  Southern 
California  and  Baja  California 

by 

Fabian  Tapia 
September  2005 


DISTRIBUTION  STATEMENT  A 

Approved  for  Public  Release 
Distribution  Unlimited 


MIT/WHOI 

2005-13 


Adult  Demography  and  Larval  Processes  in  Coastal  Benthic  Populations:  Intertidal 
Barnacles  in  Southern  California  and  Baja  California 


by 

Fabian  Tapia 

Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139 

and 

Woods  Hole  Oceanographic  Institution 
Woods  Hole,  Massachusetts  02543 

September  2005 
DOCTORAL  DISSERTATION 


Funding  was  provided  by  National  Science  Foundation  grants  OCE-0083976  and  OCE-9986627, 
through  the  WHOI  Academic  Programs  Office,  and  through  a  Presidential  Fellowship  (Beca 
Presidente  de  la  Republica)  from  the  Chilean  Ministerio  de  Planificacion  y  Cooperation. 


Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose  of  the  United  States  Government. 

This  thesis  should  be  cited  as:  Fabian  Tapia,  2005.  Adult  Demography  and  Larval  Processes  in 
Coastal  Benthic  Populations:  Intertidal  Barnacles  in  Southern  California  and  Baja  California.  Ph.D. 

Thesis.  MIT/WHOI.  2005-13. 

Approved  for  publication;  distribution  unlimited. 


Approved  for  Distribution: 


n  J.  Stegeman,  Chair 
Department  of  Biology 


Paola  Malanotte-Rizzoli 


MIT  Director  of  Joint  Program 


WHOI  Dean  of  Graduate 
Studies 


Adult  Demography  and  Larval  Processes  in  Coastal  Benthic 
Populations:  Intertidal  Barnacles  in  Southern  California 
and  Baja  California 

by 

Fabian  J.  Tapia 

B.S.,  University  of  Concepcion,  Chile,  1 994 
Submitted  in  partial  fulfillment  of  the  requirements  for  the  degree  of 


Doctor  of  Philosophy 
at  the 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 

and  the 

WOODS  HOLE  OCEANOGRAPHIC  INSTITUTION 

September  2005 

©  2005  Fabian  J.  Tapia 
All  rights  reserved. 


The  author  hereby  grants  to  MIT  and  WHOI  permission  to  reproduce  p £per  and  electronic  copies  of  this 
thesis  in  whole  or  in  part  and  to  distribute  thpfn  mrolicly 


Signature  of  Author _ _  _ 

Joint  Progran>rrrT5ceari(^graphy/Ap^H4«4^0cean  Science  and  Engineering 

Massachusetts  Institute  of  Technology 
and  Woods  Hole  Oceanographic  Institution 

September  2005 


Certified  by 


Accepted  by 


Jesus  Pineda 
Thesis  Supervisor 


John  Waterbury 

Chair,  Joint^cyhrnittee  for  Biological  Oceanography 
oods  Hole  Oceanographic  Institution 


2 


To  my  parents  Olga  and  Luis  Alberto. 

Their  believing  in  me  taught  me  how  to  believe  in  myself. 


3 


Adult  Demography  and  Larval  Processes  in  Coastal  Benthic  Populations: 
Intertidal  Barnacles  in  Southern  California  and  Baja  California 


by 


Fabian  J.  Tapia 

Joint  Program  in  Oceanography/Applied  Ocean  Science  and  Engineering 
Massachusetts  Institute  of  Technology  and  Woods  Hole  Oceanographic  Institution 

September  2005 


ABSTRACT 


The  geographic  distribution  and  dynamics  of  coastal  benthic  populations  are  shaped  by  physical  - 
biological  interactions  affecting  larval  dispersal  and  the  demography  of  juvenile  and  adult 
individuals.  This  thesis  focused  on  nearshore  patterns  of  larval  distribution  and  regional  patterns 
in  demography  of  intertidal  barnacles  in  Southern  and  Baja  California.  Horizontal  and  vertical 
distributions,  and  the  mortality  rates  of  larvae,  were  assessed  from  short  term  (i.e.  days)  small- 
scale  observations  (0.1-1  km)  in  nearshore  waters.  Observations  on  spatial  variability  of  adult 
barnacle  demography  were  gathered  over  1.5  years  at  scales  of  hundreds  of  kilometers. 

Stage-specific  horizontal  distributions  and  nearshore  current  measurements  suggested  that  larvae 
of  Balanus  giandula  and  Chthamalus  spp.  may  experience  limited  dispersal.  High  mortality  rates 
could  further  limit  travel  distances  and  the  exchange  of  individuals  among  disjunct  populations. 
Data  on  vertical  distributions  indicated  that  nauplii  and  cyprids  of  Balanus  nubilus  and  Pollicipes 
polymerus  occur  at  different  depths.  Nauplii  remained  near  the  surface  at  all  times,  whereas 
cyprids  occurred  in  the  bottom  half  of  the  water  column.  Such  distributions,  combined  with 
vertical  variability  in  horizontal  flows,  might  cause  the  observed  horizontal  segregation  of  nauplii 
and  cyprids. 

Differences  in  survival,  growth  rate,  size  structure,  and  per  capita  fertility  of  adult  Balanus 
giandula  were  observed  between  Dana  Point  (Southern  California)  and  Punta  Baja  (Baja 
California),  a  site  located  near  the  species’  southern  limit  of  distribution.  Effects  of  spatial 
differences  in  demography  on  population  persistence  were  assessed  with  a  stage-structured  matrix 
model.  Model  analyses  indicated  that  the  Punta  Baja  population  is  more  susceptible  to 
environmental  stochasticity  and  more  prone  to  local  extinction  than  populations  located  further 
north. 

This  thesis  emphasizes  the  importance  of  characterizing  factors  that  affect  the  dynamics  of 
benthic  populations  at  multiple  spatial-temporal  scales,  and  the  usefulness  of  small  scale  high- 
frequency  observations  of  nearshore  phenomena,  especially  in  relation  with  the  dispersal  of 
larvae. 


Thesis  Supervisor  :  Jesus  Pineda,  Associate  Scientist 

Biology  Department 
Woods  Hole  Oceanographic  Institution 
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CHAPTER  1 


GENERAL  INTRODUCTION 


1.1.  Background  and  motivation  of  this  study 


Coastal  ecosystems  have  been  modified  by  human  activities  over  the  past 
centuries  (Vitousek  et  al.  1997;  Jackson  et  al.  2001).  A  large  fraction  of  the  world 
population  is  concentrated  around  coastal  areas,  which  are  heavily  exploited  for  food, 
recreation,  and  waste  disposal.  Marine  biological  communities  have  been  directly 
affected  by  these  uses,  especially  by  overfishing  and  habitat  degradation  (Botsford  et  al. 
1997;  Jackson  et  ah  2001),  and  will  likely  be  impacted  by  future  changes  in  global 
climate  (Sagarin  et  al.  1999).  Such  changes  may  affect  the  distribution  and  abundance  of 
species,  altering  the  composition  of  coastal  communities  and  making  some  species  more 
prone  to  local  extinction,  while  facilitating  the  colonization  of  others.  Many  species  have 
disappeared  from  coastal  habitats  due  to  natural  and  human-induced  changes  in  the 
physical  and  biological  environment  (Jackson  et  al.  2001 ;  Lotze  and  Milewski  2004),  and 
reports  of  invasive  or  exotic  species  in  coastal  habitats  have  recently  increased  (Ruiz  et  al. 
2000;  Orensanz  et  al.  2002;  Castilla  et  al.  2005).  Marine  reserves  have  been  established 
over  the  last  decade  as  a  means  to  protect  habitats  and  biodiversity,  as  well  as  to  manage 
marine  resources  under  exploitation  (Dayton  et  al.  2000;  Lubchenco  et  al.  2003). 
Underlying  the  concept  of  a  marine  reserve  is  the  principle  that  an  area  can  be  defined 
such  that  it  encompasses  the  spatial  scale  at  which  population  renewal  takes  place.  The 
size  of  this  “home  range”  is  given  by  the  scale  of  dispersal  of  its  individuals,  which  in 
turn  is  a  function  of  a  species’  life  cycle  and  the  scale  of  dominant  physical  processes 
(Shanks  et  al.  2003). 

Many  species  of  coastal  invertebrates  and  fish  produce  larvae  that  spend  minutes 
to  months  in  the  plankton  (Levin  and  Bridges  1995),  where  they  undergo  morphological, 
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behavioral  and  physiological  changes  before  reaching  a  stage  that  is  competent  to 
metamorphose  into  a  juvenile  individual.  During  their  time  in  the  plankton,  larvae  may  be 
transported  by  currents  and  disperse  over  long  distances  (Scheltema  1968,  1986),  joining 
the  adult  population  at  localities  that  are  often  far  from  their  parental  sites.  Thus,  the 
supply  and  settlement  of  larvae  onto  a  given  population  may  be  decoupled  from  its 
reproductive  output.  To  assess  the  extent  to  which  population  renewal  is  decoupled  from 
local  reproductive  output,  and  to  estimate  the  spatial  scales  at  which  disjunct  populations 
exchange  individuals  (i.e.  scales  of  connectivity),  it  is  necessary  to  characterize  the 
interaction  between  larval  biology  and  hydrodynamic  processes  affecting  their  dispersal. 
For  instance,  increases  in  the  rate  of  larval  mortality  with  respect  to  residual  current 
velocities  may  greatly  reduce  dispersal  distances  (Ellien  et  al.  2004).  On  the  other  hand, 
differences  in  vertical  distribution  and  swimming  behavior  may  result  in  distinct  patterns 
of  dispersal  for  different  stages  of  larval  development  under  similar  advective  regimes 
(c.g.  Bousfield  1955;  Rothlisberg  and  Miller  1983). 

Invertebrate  larvae  vary  in  size  and  morphology  (Levin  and  Bridges  1995).  Size 
ranges  from  tens  to  thousands  of  microns,  whereas  morphology  and  body  architecture 
include  a  variety  of  swimming  appendages,  calcareous  shells,  chitinous  carapaces,  spines, 
as  well  as  structures  that  increase  buoyancy,  facilitate  locomotion,  or  might  protect 
against  predators  (Morgan  1995).  A  similarly  wide  range  of  variability  is  found  for 
swimming  speeds,  which  may  range  from  less  than  1  cm  s'1  in  ciliated  larvae  to  2-8  cm  s’ 
1  in  some  crustacean  larvae  (Chia  et  al.  1984).  Horizontal  current  velocities  in  coastal 
waters  are  typically  greater  than  these  swimming  speeds  (Largier  2003,  p.  S75  and 
references  therein),  hence  most  invertebrate  larvae  are  unable  to  determine  their  position 
by  horizontal  swimming.  However,  vertical  shear  in  horizontal  currents  may  determine 
specific  patterns  of  circulation  and  generate  differences  in  the  dispersal  patterns  of  larvae 
occurring  at  different  depths  in  the  water  column  (Hannah  et  al.  1997;  Hill  1998).  For 
instance,  it  has  been  suggested  that  a  number  of  coastal  and  estuarine  larvae  may  take 
advantage  of  vertical  shear  and  control  their  horizontal  displacements  by  moving  up  and 
down  the  water  column  at  various  times  of  the  tidal  or  diurnal  cycle  (e.g.  Bousfield  1955; 
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Forward  1988;  DiBacco  et  al.  2001;  Poulin  et  al.  2002).  Neustonic  larvae  may  be 
transported  by  onshore-moving  warm  fronts  (Pineda  1994),  whereas  the  larvae  of  other 
taxa  may  be  transported  by  responding  to  transient  hydrodynamic  features  such  as 
internal  bores  (Pineda  1999;  Helfrich  and  Pineda  2003).  Thus,  there  are  many  cases 
where  larvae  do  not  behave  as  passive  particles,  and  dispersal  cannot  be  predicted  simply 
from  mean  flows  without  information  on  the  behavior  and  vertical  distribution  of  larvae. 

It  has  been  suggested  that  variability  in  mortality  rates  may  greatly  constrain  the 
dispersal  distances  of  invertebrate  larvae  (e.g.  Ellien  et  al.  2004).  However,  despite  its 
importance  for  dispersal  and  connectivity,  available  estimates  of  larval  mortality  remain 
scarce  and  highly  uncertain  (Rumrill  1990;  Morgan  1995).  Larval  mortality  can  be 
estimated  directly  when  larval  aggregations  are  tracked  over  time,  using  the  observed 
change  in  larval  abundance  to  compute  a  mortality  rate  averaged  over  that  period. 
Although  such  tracking  may  be  feasible  in  closed  or  semi-enclosed  environments  (e.g. 
Lamare  and  Barker  1999;  Arnold  et  al.  2005),  in  open  coastal  waters  it  requires 
observations  at  spatial  and  temporal  scales  that  are  seldom  achieved  (e.g.  Natunewicz  et 
al.  2001).  Thus,  laboratory  (Johnson  and  Brink  1998)  and  indirect  field  methods  are 
usually  employed  to  estimate  mortality  from  ratios  of  local  larval  production  to 
recruitment  (e.g.  Thorson  1950;  Connell  1970)  or  from  information  on  the  stage 
composition  of  larvae  over  a  given  area  (Aksnes  et  al.  1997). 

While  the  physical  and  biological  processes  that  affect  larval  dispersal  and 
mortality  may  define  the  potential  limits  of  a  species  geographic  range,  the  factors  that 
ultimately  determine  its  persistence  and  demography  at  a  given  site  are  those  affecting 
post-settlement  survival  (see  Gosselin  and  Qian  1997)  and  the  subsequent  recruitment  of 
settlers  into  reproductive  adults.  Regional  differences  in  these  factors  must  affect  the 
spatial  variability  of  a  population’s  structure  and  vital  rates  and,  by  extension,  its 
geographic  distribution. 
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1.2.  Synopsis 


The  motivation  of  this  thesis  is  to  understand  the  processes  that  affect  the 
nearshore  distribution  of  larval  invertebrates  and  their  potential  effect  on  distribution, 
connectivity,  and  persistence  of  coastal  benthic  populations.  Intertidal  barnacles  were 
chosen  as  model  species  because  of  their  abundance,  accessibility  and  wide  distribution 
in  coastal  environments,  and  the  extensive  body  of  knowledge  on  their  morphology, 
reproductive  biology,  larval  development  and  population  dynamics  accumulated  in  the 
literature  for  over  100  years  (e.g.  Darwin  1854;  Hatton  1938;  Barnes  and  Barnes  1956; 
Crisp  1962;  Hines  1979;  Crisp  and  Bourget  1985;  Barnes  1989;  Crisp  et  al.  1991;  Barnes 
1999;  Wethey  2002).  Results  presented  in  following  chapters  focus  on  the  larvae  and 
adults  of  the  acorn  barnacle  Balanus  glandula  Darwin,  1854,  and  on  the  larvae  of 
Chthamalus  fissus  Darwin  1854,  Pollicipes  polymerns  Sowerby,  1833,  and  Balanus 
nubilus  Darwin,  1854.  The  first  three  species  are  numerically  dominant  on  intertidal 
habitats  in  Southern  California  and  Baja  California.  Balanus  glandula  occurs  in  the  high 
to  middle  intertidal  zone  along  the  west  coast  of  North  America  from  the  Aleutian  Islands 
to  San  Quintin,  Baja  California  (Newman  and  Abbott  1980).  Chthamalus  fissus  is  found 
in  the  high  and  upper  middle  intertidal  zone  between  San  Francisco  and  Baja  California 
(Newman  and  Abbott  1980).  Its  distribution  occasionally  overlaps  with  that  of 
Chthamalus  dalli  Pilsbry  1916,  a  northern  species  found  from  Alaska  to  a  southern  limit 
that  fluctuates  between  Point  Conception  (Wares  2001)  and  San  Diego  (Newman  and 
Abbott  1980).  Pollicipes  polymerus  is  found  in  the  middle  intertidal  zone  of  wave-swept 
rocky  shores  from  British  Columbia  to  Punta  Abreojos,  Baja  California  (Newman  and 
Abbott  1980),  whereas  Balanus  nubilus  occurs  from  the  lower  intertidal  to  90  m  depths, 
and  from  La  Jolla,  Southern  California,  to  Alaska  (Newman  and  Abbott  1980). 

Observations  on  the  adults  and  larvae  of  these  species  were  gathered  at  regional 
(10-100  km)  and  local  scales  (0.1-1  km),  and  aimed  at  addressing  questions  regarding  the 
effect  of  nearshore  physical  processes  on  larval  distribution,  dispersal  and  mortality,  as 
well  as  the  geographic  distribution  and  persistence  of  adults. 
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In  Chapter  2, 1  attempt  to  characterize  the  small-scale  distribution  and  mortality  of 
larval  barnacles  in  nearshore  waters.  Stage-specific  spatial  distributions  of  Balanus 
glandula  and  Chthamalus  spp.  larvae  were  observed  daily  over  seven  days.  The  results 
are  used  to  discuss  the  potential  for  dispersal  of  larvae  in  nearshore  waters.  Chapter  3 
presents  observations  on  the  vertical  distribution  of  larvae  in  nearshore  waters  over  a  48 
hour  period,  showing  that  different  stages  of  development  occur  at  different  depths  in  the 
water  column.  These  patterns  do  not  change  between  day  and  night  and,  combined  with 
vertical  variability  in  horizontal  flows,  may  explain  the  horizontal  segregation  of  nauplii 
and  cyprids  observed  in  Chapter  2.  The  role  of  onshore  winds  as  a  physical  mechanism 
forcing  the  nearshore  distribution  of  surface  material,  and  possibly  of  neustonic  larvae,  is 
discussed  in  the  publication  included  as  Appendix  A. 

Chapter  4  focuses  on  a  regional  analysis  of  the  population  structure  and 
demography  of  Balanus  glandula  in  Southern  California  and  Baja  California.  Field 
observations  on  size  structure,  growth  rates,  survival  and  per  capita  fertility  were 
combined  with  a  stage-structured  matrix  population  model  to  assess  the  effect  that  spatial 
changes  in  demography  and  environmental  variability  may  have  on  population 
persistence  and  geographic  distribution. 

Finally,  general  conclusions  from  the  results  of  this  thesis  and  recommendations 
for  future  research  are  presented  on  Chapter  5. 
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2.  CHAPTER  2 

LIMITED  DISPERSAL  AND  HIGH  MORTALITY  OF  BARNACLE 
LARVAE  IN  THE  NEARSHORE 


2.1.  Abstract 


The  stage-specific  spatial  distribution  and  mortality  of  Balanns  glandula  and 
Chthamalus  spp.  larvae  were  assessed  from  a  series  of  daily  vertical  plankton  tows 
collected  off  La  Jolla,  Southern  California,  on  6-12  March  2003.  Sampling  stations  were 
located  within  1.1  km  of  the  shoreline,  at  depths  of  10  -  45  meters.  Vertical  distributions 
of  temperature  and  current  velocities  were  recorded  at  an  average  depth  of  16  m.  For  both 
species,  a  spatial  segregation  of  naupliar  stages  and  cyprids  was  observed,  although  the 
differences  were  statistically  significant  for  Chthamalus  spp.  only.  Earlier  naupliar  stages 
(NII-NIII)  were  more  abundant  at  the  inshore  stations,  whereas  later  stages  (NIV-NVI) 
were  more  abundant  at  the  offshore.  Cyprid  concentrations  were  higher  at  the  inshore 
stations.  These  structured  spatial  distributions  suggest  limited  dispersal  of  intertidal 
barnacle  larvae  in  the  nearshore.  Trajectories  computed  from  current  measurements,  and 
a  set  of  assumptions  regarding  the  vertical  distribution  of  larvae,  predicted  short  dispersal 
distances  (2-8  km)  over  a  7-day  period.  Vertical  life  tables  were  used  to  estimate  naupliar 
mortality  from  daily  stage  distributions.  Average  estimates  for  the  instantaneous  rate  of 
larval  mortality  in  B.  glandula  and  Chthamalus  spp.  were  0.329  d'1  (SD=0.045  d'1)  and 
0.232  d1  (SD=0.033  d"1),  respectively.  These  mortality  estimates  are  substantially  higher 
than  previously  assumed  for  these  species.  I  discuss  the  limitations  of  these  data,  as  well 
as  the  implications  of  limited  dispersal  and  high  larval  mortality  for  the  connectivity 
among  populations  of  intertidal  barnacles  and  other  coastal  species  with  a  similar  life 
cycle. 


19 


2.2.  Introduction 


Many  coastal  benthic  invertebrates  have  complex  life  cycles.  Sedentary  adults 
release  larvae  into  the  pelagic  environment,  where  development  is  completed  within  days 
to  months  (Levin  and  Bridges  1995).  Larvae  may  be  transported  over  long  distances 
during  their  development  (Scheltema  1968,  1986),  which  often  results  in  recruitment 
being  decoupled  from  adult  abundance  at  any  given  site.  Coastal  populations  are  thus 
considered  open  in  terms  of  their  reproductive  output  and  local  demography.  To  assess 
how  much  of  the  local  recaiitment  is  determined  by  local  reproductive  output,  the  rates  of 
dispersal  and  mortality  of  early  larval  stages  must  be  known.  The  speed,  trajectory  and 
vertical  variability  of  horizontal  currents  play  an  important  role  on  larval  dispersal, 
especially  for  species  whose  vertical  distribution  changes  over  the  course  of  development 
(e.g.  Bousfield  1955). 

Intertidal  barnacles  are  conspicuous  members  of  coastal  communities  around  the 
world.  Their  high  abundance,  wide  geographical  distribution,  and  sessile  adult  stage 
make  them  good  model  species  for  the  study  of  interactions  between  larval  dispersal  and 
the  openness  of  a  population.  Past  studies  have  focused  on  mesoscale  patterns  of  adult 
and  larval  distribution  (i.e.  tens  to  hundreds  of  km),  as  well  as  their  dependence  on 
physical  variability  at  similar  scales.  Mathematical  models  that  couple  benthic  and 
oceanic  processes  have  been  used  in  a  series  of  influential  studies  (Roughgarden  et  al. 
1988;  Possingham  and  Roughgarden  1990;  Alexander  and  Roughgarden  1996;  Connolly 
and  Roughgarden  1998;  Gaylord  and  Gaines  2000)  to  assess  the  effect  of  biological  and 
physical  variables  on  mesoscale  patterns  of  barnacle  distribution.  Central  to  all  of  these 
studies  is  the  idea  that  larvae  are  swept  away  from  their  parental,  intertidal  populations  by 
along-shore  and  cross-shore  advection,  determined  by  mesoscale  coastal  circulation  and 
Ekman  transport,  respectively  (Roughgarden  et  al.  1988).  In  this  framework,  the 
openness  and  connectivity  of  intertidal  barnacle  populations  are  determined  by  the 


20 


survival  of  larvae  in  the  plankton,  the  intensity  of  onshore  transport,  and  substrate 
availability. 

Extensive  research  has  been  done  on  the  influence  of  substrate  type  and 
availability  on  larval  settlement  and  recruitment  (e.g.  Raimondi  1988;  Pineda  1994b; 
Pineda  and  Caswell  1997),  as  well  as  the  physical  mechanisms  that  mediate  the  onshore 
transport  of  larvae.  Patterns  of  larval  settlement  have  been  linked  to  the  occurrence  of 
internal  waves  (Shanks  1983),  internal  tidal  bores  (Pineda  1991,  1994a),  persistent 
onshore  winds  (Bertness  et  al.  1996),  and  the  cross-shore  displacement  of  surface  waters 
during  relaxation  of  wind-induced  upwelling  (Farrel  et  al.  1991;  Wing  et  al.  1995). 
Contrastingly,  the  survival  and  dispersal  of  early  larval  stages  in  the  nearshore  have  not 
been  investigated  extensively  (but  see  Natunewicz  et  al.  2001),  although  retention  and 
self-recruitment  might  be  more  prevalent  than  previously  thought  (Cowen  et  al.  2000; 
Sponaugle  et  al.  2002;  Warner  and  Cowen  2002).  Assessing  the  scales  of  dispersal  and 
the  mortality  of  larvae  in  the  nearshore  is  thus  necessary  for  understanding  the 
interactions  between  nearshore  hydrodynamics  and  the  openness  of  coastal  populations. 

In  this  contribution  1  examine  the  small-scale  (~1  km  from  shore)  horizontal 
distribution  and  estimate  the  mortality  rate  of  larva]  Balanus  glcindula  and  Chthamalus 
spp.,  the  two  most  abundant  species  of  intertidal  barnacles  on  the  US  West  Coast. 
Persistent  patterns  of  spatial  distribution  were  observed,  with  a  horizontal  segregation  of 
late  naupliar  stages  and  cyprids.  These  observations  and  a  series  of  particle  trajectories 
predicted  from  nearshore  current  velocities  suggest  a  limited  larval  dispersal  larvae. 
Larval  mortality  estimates  for  both  species  were  substantially  higher  than  the  rates 
previously  assumed  (see  Gaylord  and  Gaines  2000,  p.  774).  I  discuss  the  limitations  of 
these  results  in  relation  to  local  physical  variability  and  larval  behavior,  as  well  as  the 
assumptions  underlying  the  estimation  of  larval  mortality  rates. 
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2.3.  Materials  and  Methods 


2.3.7.  Study  site  and  species 

The  survey  was  conducted  in  nearshore  waters  off  the  Scripps  Institution  of 
Oceanography  (SIO)  pier  in  La  Jolla,  Southern  California,  between  6-13  March  2003 
(hereafter  Days  1-8  of  the  survey).  Intertidal  barnacles  on  which  this  study  was  based  are 
two  of  the  most  abundant  species  found  on  the  US  West  Coast.  The  acorn  barnacle 
Balanus  glandula  Darwin,  1854  occurs  in  the  high  to  middle  intertidal  zone  along  the 
west  coast  of  North  America  from  the  Aleutian  Islands  to  San  Quintin,  Baja  California 
(Newman  and  Abbott  1980).  Chthamalus  fissus  is  found  between  San  Francisco  and  Baja 
California  (Newman  and  Abbott  1980),  and  dominates  the  high  to  upper  middle  intertidal 
in  La  Jolla.  Its  distribution,  however,  occasionally  overlaps  with  that  of  C.  dalli ,  a 
northern  species  found  from  Alaska  to  a  southern  limit  that  fluctuates  between  Point 
Conception  (Wares  2001)  and  San  Diego  (Newman  and  Abbott  1980).  Although  adult  C. 
fissus  and  C.  dalli  can  be  identified  by  dissection  and  microscopic  examination  (Newman 
and  Abbott  1980),  their  larvae  are  morphologically  identical  (Miller  et  al.  1989).  Thus, 
chthamalid  larvae  will  be  referred  to  as  Chthamalus  spp.,  as  larval  identification  was 
based  solely  on  morphological  criteria  (see  below). 

Timing  of  the  survey  was  based  on  the  known  reproductive  cycle  of  B .  glandula 
and  C.  fissus ,  and  aimed  at  finding  their  larvae  in  nearshore  waters.  Chthamalus  settlers 
can  be  found  year  round  at  a  wide  range  of  tidal  heights  in  Central  and  Southern 
California  (Hines  1979;  Pineda  1994b  and  unpublished  observations).  Individuals  reach 
sexual  maturity  ca.  2  months  after  settlement,  after  which  up  to  16  small  broods  are 
continuously  produced  each  year  (Hines  1976).  Adult  B.  glandula ,  in  contrast,  produce  2 
to  6  broods  over  a  reproductive  season  that  extends  from  early  winter  to  late  spring 
(Barnes  and  Barnes  1956;  Hines  1979;  Newman  and  Abbott  1980).  Settlement  takes 
place  between  late  winter  and  early  summer  (Hines  1979  and  author’s  unpublished 
observations). 
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2. 3. 2.  Sampling  procedure 


Plankton  samples  were  collected  daily  over  7  days  from  3  stations  located  within 
0.3,  0.6  and  1.1  km  of  the  shoreline,  at  depths  of  10  -  40  m  (Fig.  2.1).  Plankton  was 
collected  by  vertically  towing  a  1 10  pm  mesh  net  (0.75  m  diameter,  2.75  m  length)  from 
bottom  to  surface.  The  water  volume  sampled  at  each  station  was  estimated  using  a 
mechanical  flowmeter  (General  Oceanics  2030).  Samples  were  taken  into  the  laboratory 
and  preserved  in  95%  ethanol. 

Concurrently  with  the  plankton  survey,  I  monitored  settlement  of  B.  glandula  and 
Chthamalus  spp.  larvae  on  collectors  deployed  within  the  vertical  range  of  adult 
distribution  at  a  nearby  intertidal  site  (Dike  Rock,  Fig.  2.1).  Collectors  were  made  out  of 
11 -cm  long  pieces  of  white  PVC  pipe  (2.54  cm  diameter),  which  were  cut  in  half 
lengthwise  and  grooved  in  three  areas.  A  hole  was  drilled  in  the  center,  so  that  the 
collectors  could  be  attached  to  the  substrate  using  stainless  steel  screws  that  were 
cemented  into  holes  drilled  in  the  rock.  Settlement  collectors  were  first  deployed  on  Day 
1  and  recovered/replaced  daily  during  the  daytime  low  tide  until  Day  8.  Tidal  height  and 
wave  action  precluded  the  recovery  of  collectors  on  Day  4.  Barnacles  on  recovered 
collectors  were  identified  and  enumerated  using  a  dissecting  microscope.  Only  larvae  that 
settled  within  ca.  0.16  mm  of  the  intersection  between  the  groove  bottom  and  wall  were 
counted,  since  settlement  rarely  occurs  outside  this  area.  Hence,  the  area  of  substrate  per 
collector  was  -1.9  cm2  (Pineda  1994b).  Settlement  rates  were  computed  as  number  of 
settlers  per  day  of  deployment. 


2.3.3.  Identification  of  larvae  and  analysis  of  spatial  patterns 

The  classification  of  larvae  by  species  and  stage  was  based  on  morphological 
descriptions  for  Balanus  glandula  (Branscomb  and  Vedder  1982;  Brown  and 
Roughgarden  1985)  and  Chthamalus  spp.  larvae  (Miller  et  al.  1989),  as  well  as  for  other 
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species  that  occur  in  the  study  area  (Lewis  1975;  Miller  and  Roughgarden  1994). 
Abundances  of  naupliar  stages  Nil  through  NVI  are  reported  only,  as  NI  larvae  molt  into 
Nil  within  hours  of  hatching  at  the  temperatures  recorded  in  the  survey  area  (Brown  and 
Roughgarden  1985;  Miller  et  al.  1989). 

Spatial  patterns  in  the  distribution  of  different  larval  stages  were  analyzed  using  a 
two-way  Analysis  of  Variance  (ANOVA)  without  replicates,  in  which  sampling  stations 
and  days  were  used  as  factors.  To  achieve  homogeneity  of  variances  prior  to  the  analysis, 
concentrations  of  nauplii  and  cyprids  were  transformed  as  log(x+l)  and  -l/(x-H), 
respectively. 


2.3.4.  Environmental  variables 

From  Day  1,  temperature  in  the  water  column  was  recorded  at  1  minute  intervals 
with  Onset  Stowaway  XTI  loggers  (response  time  <  1 5  sec),  located  at  the  surface  and  at 
1,  2,  4,  6,  8,  10,  12  and  14  meters  above  the  bottom  (mab).  The  string  of  loggers  was 
deployed  at  a  depth  of  ca.  18  m  (Fig.  2,1).  Florizontal  and  vertical  currents  were 
measured  with  a  bottom-mounted  1,200  kHz  Acoustic  Doppler  Current  Profiler  (ADCP, 
RD  Instruments,  USA),  deployed  on  Day  1  near  Station  2  at  a  depth  of  ca.  16  m  (Fig.  2.1). 
Measurements  were  recorded  in  1  meter  bins  at  1-min  intervals  (60  pings  per  ensemble). 
Due  to  the  temporal  variability  in  water  depth  introduced  by  tidal  fluctuations,  surface 
wave  action,  and  side  lobbing,  data  collected  from  the  surface  bins  had  to  be  discarded. 
The  uppermost  bin  with  reliable  data  was  at  14  mab,  which  was  on  average  2  m  below 
the  surface. 

Data  on  hourly  wind  velocity  and  direction  were  obtained  from  the  SIO  Coastal 
Data  Information  Program  website  (http://cdip.ucsd.edu).  The  wind  gauge  (Qualimetrics 
Skyvane  anemometer)  was  located  at  20.2  m  above  Mean  Lower  Low  Water  (MLLW), 
on  the  west  end  of  the  SIO  Pier.  Wind  data  were  rotated  and  aligned  to  the  average 
shoreline  orientation  (21°  with  respect  to  true  north)  and  then  decomposed  into  along- 
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shore  and  cross-shore  components.  ADCP  data  were  aligned  to  the  main  axis  of 
variability  for  horizontal  currents  (24.7°  with  respect  to  true  north),  which  roughly 
paralleled  the  shoreline  orientation.  The  orientation  of  this  axis  was  given  by  the  major 
eigenvector  of  a  covariance  matrix  computed  for  depth-averaged  east  and  north  velocities. 

2.3.5.  Assessment  of  potential  dispersal  -  Progressive  Vector  Diagrams 

Hypotheses  regarding  the  nearshore  distribution  of  barnacle  larvae  were 
formulated  using  Progressive  Vector  Diagrams  (PVD).  The  construction  of  these 
diagrams  is  based  on  the  strong  assumption  that  horizontal  currents  measured  at  a  single 
point  are  representative  of  a  larger  area.  PVDs  were  used  as  a  first-order  approach  to 
project  the  trajectory  of  particles  found  in  nearshore  waters  at  the  study  area.  The  goal  of 
this  analysis  was  to  assess  how  likely  it  is  for  a  particle  to  remain  in  nearshore  waters 
within  the  time  scale  of  this  study  (7  days).  Near-surface  current  velocities  from  moored 
current  meters  might  fail  to  predict  the  trajectory  of  particles  found  in  the  uppermost 
layer  of  the  water  column  (Tapia  et  al.  2004,  Appendix  A).  However,  it  was  assumed  that 
the  depth  range  for  which  reliable  current  meter  data  were  available  (2-14  mab)  was 
representative  of  the  vertical  range  of  larval  distribution. 

PVDs  were  computed  as  predicted  (x,y)  positions  at  time  t„  : 


n 


(2.1) 


1=1 


where  u(t)  and  v(f)  are  the  east  and  north  current  speeds  measured  by  the  current  meter  at 
a  given  depth  at  time  t ,  and  At  is  the  sampling  interval  (Emery  and  Thomson  1998,  p. 
1 65).  PVDs  were  computed  for  the  duration  of  the  study,  and  for  six  different  patterns  of 
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larval  vertical  distribution.  First,  it  was  assumed  that  larvae  are  uniformly  distributed  in 
the  water  column,  and  their  trajectory  was  predicted  using  depth-averaged  current 
velocities.  For  the  second  through  fourth  alternative  patterns,  it  was  assumed  that  larvae 
remain  at  mid  depths  (7-9  mab),  near  the  bottom  (2-4  mab),  and  near  the  surface  (12-14 
mab)  at  all  times.  In  the  last  two  cases,  larvae  were  assigned  behavior  in  the  form  of 
diurnal  vertical  migrations  (DVM),  so  that  they  occur  near  the  bottom  during  the  day 
(defined  as  6:00-18:00  PST)  and  near  the  surface  at  night.  In  one  case  (Two-layer  DVM  ) 
larvae  were  assumed  to  move  between  the  upper  and  lower  half  of  the  water  column, 
whereas  in  the  final  example  (Extreme  DVM)  larvae  moved  between  a  3-m  thick  layer 
near  the  surface  (12-14  mab)  and  a  layer  of  similar  thickness  near  the  bottom  (2-4  mab). 


2.3.6 .  Estimation  of  larval  mortality 

Vertical  life  tables  (Aksnes  and  Ohman  1996)  were  used  to  estimate  the  mortality 
of  nauplius  larvae.  The  method  yields  estimates  of  mortality  at  the  transitions  between 
contiguous  larval  stages,  and  is  recommended  for  cases  where  advection  is  expected  to 
affect  the  horizontal  distribution  of  larvae  (Aksnes  et  al.  1997).  Mortality  is  estimated 
from  instantaneous  stage  distributions  (i.e.  obtained  at  one  point  in  time),  hence  it  is  not 
necessary  to  track  a  particular  larval  aggregation.  Instead,  the  stage  distribution  must  be 
representative  of  the  population  under  study,  and  spatial  coverage  and  resolution  should 
be  large  enough  to  compensate  for  the  effect  of  small-scale  patchiness. 

Derivations  of  the  method  presented  by  Aksnes  and  Ohman  (1996,  p.  1462)  are 
summarized  below.  Three  main  assumptions  are  made  :  (1)  daily  recruitment  (p,)  to  a 
given  stage  i  is  constant  over  the  duration  of  that  stage,  (2)  the  duration  of  a  stage  (cc,)  is 
constant  and  equal  for  all  individuals  in  the  same  stage,  and  (3)  the  mortality  over  that 
period  (0,)  is  constant.  Thus,  the  number  of  nauplii  in  stage  i  at  time  t  =  x  can  be 
expressed  as: 
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(2.2) 


/?,.  =  p,  Jexp[-  Or(x  -  t)]dt  =  ^[l-exp(-^a()] 

x-a,  ' 

which  is  equal  to  the  number  of  larvae  that  recruited  during  the  last  a,  days  and  survived 
(Aksncs  and  Ohman  1996).  If  it  is  also  assumed  that  the  mortality  rate  is  constant  for  a 
period  equal  to  the  duration  of  two  consecutive  stages  i  and  i+l,  then  the  number  of 
nauplius  larvae  in  stage  i+\  at  day  x  is  : 

x 

"m  =  Pm  Jexp[-0(x -t)]dt  =  ^[l-exp (-0tf,+1)]  (2.3) 

The  number  of  cyprid  larvae,  for  which  the  stage  duration  ac  is  assumed  infinite 
by  the  analysis,  can  be  expressed  as  : 

* 

nc  =  Pc  jexp[-0c (*-/)]<//  =  Y  (2.4) 

—oo  ( 

Finally,  the  rate  of  recruitment  to  a  stage  ;+l  is  the  product  of  recruitment  to  the 
previous  stage  i  and  the  stage-specific  survival  : 

Pm  =Aexp(-fy*<)  (2-5) 

An  equation  that  relates  the  relative  abundance  of  two  consecutive  stages  with 
their  duration  and  mortality  rate  can  be  obtained  by  combining  Eqs.  2.2,  2.3,  and  2.5,  and 
by  setting  Q,  -  6 : 


3_=  exp((ta,)-1 
”m  1-exp  (~0aM) 

Similarly,  an  equation  that  relates  the  number  of  nauplius  VI  and  cyprid  larvae 
can  be  obtained  by  combining  Eqs.  2.2,  2.4,  and  2.5. 
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(2.7) 


Equation  2.6  was  solved  iteratively  to  obtain  estimates  of  the  mortality  rate  at  the 
transition  between  naupliar  stages  NII-NV.  The  estimated  mortality  of  NVI  larvae  was 
obtained  directly  from  Equation  2.7.  Estimates  of  the  number  of  individuals  per  stage  for 
each  day  with  plankton  observations  (Days  1-7)  were  obtained  by  pooling  larval  counts 
recorded  at  the  three  stations.  An  average  stage  duration  of  3  days  was  used  for 
Chthamalus  spp.  (Miller  et  al.  1989),  whereas  durations  of  1,  2,  2,  3,  and  3  days  were 
used  for  B.  glandula  Nil,  NIII,  NIV,  NV,  and  NVI,  respectively  (Brown  and 
Roughgarden  1985). 
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2.4.  Results 


2.4.1.  Spatial-temporal  distribution  of  larval  stages 

Patterns  of  spatial  distribution  for  different  larval  stages  were  consistent  for  both 
species.  In  general,  earlier  naupliar  stages  (Nil  and  NIII)  were  more  abundant  at  the 
inshore  stations,  whereas  later  stages  (NV-NVI)  increased  in  abundance  towards  the 
offshore  station  (Fig.  2.2,  2.3).  The  statistical  analysis  of  stage-specific  distributions 
indicated  no  significant  between-station  differences  for  all  but  one  larval  stage  in  Balanus 
glandula  (Nil,  Table  2.1).  Conversely,  the  analysis  of  Chthamalus  spp.  distributions 
indicated  significant  spatial  differences  in  all  but  one  case  (Nil,  Table  2.1).  The 
concentration  of  B.  glandula  Nil  larvae  at  the  offshore  station  was  significantly  lower 
than  at  the  inshore  stations  (Table  2.1).  On  the  other  hand,  concentrations  of  Chthamalus 
spp.  NIV,  NV,  and  NVI  larvae  were  significantly  higher  at  the  offshore  station.  Although 
not  statistically  significant  for  B .  glandula ,  cyprid  concentrations  were  higher  at  the 
inshore  stations  (Fig.  2.2f,  2.3f)  and  resembled  the  distributions  of  Nil  rather  than,  for 
instance,  NVI. 

High  concentrations  of  cyprids  relative  to  NVI  larvae  at  the  inshore  station  suggest 
that  cyprids  were  advected  from  an  adjacent  source,  or  that  their  distribution  was  dictated 
by  a  different  set  of  physical  and  behavioral  factors.  These  high  concentrations  of  cyprids 
at  the  inshore  station  could  also  be  the  result  of  accumulation  over  an  undetermined 
number  of  days.  Patchiness  in  larval  distribution  appeared  to  be  species  specific.  For 
example,  on  Day  1  the  concentration  of  B.  glandula  Nil  decreased  steadily  from  14.5  to 
1.5  indiv  m'3  between  the  inshore  and  offshore  station  (Fig.  2.2a),  whereas  Chthamalus 
spp.  Nil  peaked  at  Station  2  with  108  indiv  m’  .  This  concentration  was  two  orders  of 
magnitude  higher  than  observed  at  the  other  two  stations  (Fig.  2.3a). 
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2.4.2.  Temporal  variability  in  stage  distribution 


Daily  stage  distributions  for  B.  glandula  and  Chthamalus  spp.  (Fig.  2.4a-g,  2.4h~n) 
were  obtained  by  pooling  larval  counts  from  the  three  sampling  stations.  All  stages  of 
larval  development  in  B.  glandula  were  observed  on  all  but  one  day  of  the  survey  (Fig. 
2.4f).  Stage  distributions  of  B.  glandula  were  dominated  by  Nil  and  NIV  during  the  first 
three  days  of  the  survey  (Fig.  2.4a-c).  The  distributions  observed  on  subsequent  days 
were  dominated  by  Nil  (Fig.  2.4d-g),  suggesting  a  constant  input  of  newly  hatched  larvae, 
either  released  by  local  adults  of  advected  into  the  study  area  from  adjacent  populations. 
Stage  distributions  of  Chthamalus  spp.  larvae  were  dominated  by  Nil  on  five  out  of  seven 
days  of  observations  (Fig.  2.4h-n),  but  most  notably  on  days  5-7  (Fig.  2.4k-n).  This 
suggests  a  continuous  input  of  larvae  released  locally  or  at  a  site  from  which  they  can  be 
advected  within  the  time  it  takes  for  the  NII-NIII  transition  to  occur  (-3  days,  Miller  et  al. 
1989).  A  high  abundance  of  cyprid  larvae  on  days  3-5  and  7,  relative  to  the  abundance  of 
NVI  larvae  (Fig.  2.4j-l,n),  suggests  that  cyprids  are  either  transported  mto  the  study  area 
from  an  external  source  or  accumulated  over  time. 


2. 4. 3.  Settlement  at  Dike  Rock 

Larvae  of  both  species  settled  at  Dike  Rock  throughout  the  survey  (Fig.  2.5), 
suggesting  that  cyprids  found  in  the  plankton  samples  were  competent  to  settle.  In 
general,  and  consistent  with  their  higher  abundance  in  the  plankton  (Fig.  2.3f), 
Chthamalus  spp.  cyprids  settled  at  higher  rates  than  those  of  B.  glandula  (Fig.  2.5a). 
While  the  timing  of  settlement  in  B.  glandula  was  not  correlated  with  changes  in  cyprid 
abundance  in  the  nearshore  (Fig.  2.5b),  daily  changes  in  the  rate  of  settlement  of 
Chthamalus  spp.  larvae  were  positively  correlated  with  daily  changes  in  their 
concentration  at  the  inshore  station  (Fig.  2.5c). 
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2. 4 . 4.  En vironmental  variables 


Wind  data  for  March  2003  indicated  that  the  survey  was  conducted  during  a 
relatively  calm  period  (Fig.  2.6a,  inset).  Cross-shore  wind  speeds  ranged  from  -2.5  to  3.2 
m  s"!  (positive  is  onshore,  111°  east  of  true  north),  whereas  along-shore  speeds  ranged 
between  -2.0  and  2.8  m  s'1  (positive  is  21°  east  of  true  north).  Fluctuations  in  cross-shore 
wind  velocities,  up  to  3  m  s"1  in  amplitude  (Fig.  2.6a),  followed  a  diurnal  cycle  that  was 
consistent  with  the  daily  sea  breeze. 

Horizontal  currents  measured  during  the  survey  were  dominated  by  along-shore 
flows  (Fig.  2.6b).  Along-shore  currents  were  more  energetic  (up  to  17  cm  s'1)  and 
uniformly  distributed  through  the  water  column  than  cross-shore  currents  (<  10  cm  s'1), 
where  most  of  the  energy  was  concentrated  near  the  surface  (Fig.  2.6c).  Variability  in  the 
along-shore  flow  was  dominated  by  the  barotropic  tide,  with  semidiurnal  changes  in 
current  direction  throughout  the  water  column  (Fig.  2.6b).  The  vertical  structure  of  cross- 
shore  flows,  on  the  other  hand,  was  consistent  with  the  structure  of  mode  one  internal 
motions,  often  showing  two  layers  of  variable  thickness  that  were  flowing  in  opposite 
directions  (Fig.  2.6c). 

Temperature  measurements  showed  a  stratified  water  column,  with  differences 
between  the  surface  and  the  bottom  that  reached  ca.  5°C  on  Day  7  (Fig.  2.6d). 
Semidiurnal  fluctuations  in  temperature  distribution,  associated  with  the  tidal  cycle, 
occurred  throughout  the  survey.  A  propagating  internal  tidal  bore  warm  front  was 
observed  in  the  afternoon  of  Day  3  (Fig.  2.6d).  The  occurrence  of  this  bore  coincided 
with  an  increase  in  the  velocity  of  onshore  currents  throughout  the  water  column, 
especially  near  the  bottom  (Fig.  2.6c),  as  well  as  a  change  in  the  vertical  structure  of 
along-shore  velocities  (Fig.  2.6b).  Although  a  high  tide  made  it  impossible  to  record 
settlement  immediately  after  such  feature  occurred  (Day  4),  maximum  Chthamalus  spp. 
settlement  was  observed  on  the  following  day  (Fig.  2.5a). 
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2.4.5 .  Progressive  Vector  Diagrams 


In  no  instance  did  the  PVD  produce  trajectories  that  suggest  an  offshore  dispersal 
of  larvae  (Fig.  2.7).  According  to  the  predicted  trajectories,  larvae  would  remain  within 
ca.  2  km  of  their  release  point  in  all  but  one  case.  Only  when  larvae  were  assumed  to 
occur  near  the  bottom  at  all  times  did  the  PVD  predict  a  longer,  northward  trajectory  that 
intersected  with  the  shoreline  -8  km  north  of  the  starting  point  (Fig.  2.7b).  Trajectories 
predicted  when  PVD  computations  were  started  at  different  phases  of  the  tidal  cycle 
(symbols  in  Fig.  2.7)  suggested  that  the  tidal  phase  during  which  particles  are  released 
should  not  affect  the  expected  range  of  dispersal  distances. 


2 . 4. 6.  Mortality  estimates 

Highly  variable  daily  stage  distributions  (Fig.  2.4)  precluded  the  estimation  of 
mortality  for  all  transitions  between  naupliar  stages.  Larval  counts  were  pooled  for  the 
contiguous  stages  NII-NIII  and  NIV-NV,  so  that  reported  mortality  estimates  correspond 
to  the  transitions  NII+NIII  ->  NIV+NV,  NIV+NV  ->  NVI,  and  NVI  C  (Table  2.2). 
Average  mortality  estimates  ranged  between  0.298  d'1  and  0.396  d’1  for  B.  glandula ,  and 
between  0.176  d"1  and  0.309  d  1  for  Chthamalus  spp..  There  was  no  trend  in  the  average 
values  or  in  the  variability  of  mortality  estimates  with  stage  transitions  (Fig.  2.8). 
Average  mortality  rates  computed  across  stage  transitions  and  days  were  0.329  d’1 
(SE=0.045  d'1)  for  if  glandula  larvae,  and  0.232  d”1  (SE=0.033  d'1)  for  Chthamalus  spp. 
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2.5.  Discussion 


I  have  shown  stage-specific  patterns  of  spatial  distribution  that  point  to  a  limited 
dispersal  of  intertidal  barnacle  larvae  in  the  nearshore,  and  estimated  mortality  rates  that 
are  substantially  higher  than  previously  assumed.  However,  large  daily  fluctuations  in  the 
observed  stage  distributions  suggest  that  the  spatial  coverage  and  resolution  of  the  survey 
were  insufficient  to  compensate  for  nearshore  advection  and  small-scale  patchiness.  Thus, 
the  basic  assumption  of  the  mortality  estimation  method  may  have  not  been  met.  Despite 
the  apparent  effect  of  advection  on  the  observed  stage  structures,  spatial  distributions  of 
late  larval  stages  suggested  that  larvae  may  complete  their  development  within  a  short 
distance  from  shore.  The  generality  of  these  results  should  be  tested  at  sites  with  different 
coastal  configuration  and  bathymetry,  as  two  submarine  canyons  flanking  the  study  area 
may  have  played  a  role  in  the  nearshore  retention  of  larvae.  Horizontal  flows  and  their 
variability  at  scales  of  kilometers,  together  with  stage-specific  patterns  of  vertical 
distribution,  must  be  better  described  before  any  conclusions  can  be  drawn  as  to  dispersal 
distances  expected  for  these  larvae. 


2.5.7.  Spatial  distribution  and  settlement  of  larvae 

Observed  temporal  changes  in  stage  distributions  suggested  a  continuous  input  of 
early  nauplii  (Nil)  into  nearshore  waters,  especially  after  Day  4  (Fig.  2.4).  Consistently 
higher  concentrations  of  Nil  larvae  at  Stations  1  and  2  could  indicate  a  continuous 
production  of  larvae  at  the  study  area,  superimposed  with  continuous  input  of  larvae 
produced  at  neighboring  populations.  Food  availability  and  water  temperature  affect  the 
rate  of  development  in  barnacle  larvae  (Barnes  and  Barnes  1958;  Scheltema  and  Williams 
1982).  Assuming  that  food  concentration  is  sufficiently  high,  it  should  take  2-3  days  for 
Nil  larvae  to  molt  into  NIII  at  the  average  temperature  of  ca.  15°C  recorded  during  this 
study  (Brown  and  Roughgarden  1985;  Miller  et  al.  1989).  During  such  time,  a  passive 
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particle  in  the  area  could  have  been  transported  2.4  -  3.6  km  northward,  based  on  a  depth- 
averaged  mean  velocity  of  1.4  cm  s_l  computed  from  alongshore  current  measurements. 
Thus,  intertidal  populations  located  within  a  few  kilometers  of  each  other  could  be 
connected  through  the  dispersal  of  early  larval  stages,  assuming  that  their  behavior 
resembles  that  of  a  passive  particle  more  closely  than  that  of  later  stages. 

Cyprid  concentrations  at  the  inshore  station  were  higher  than  those  of  NVI  larvae, 
suggesting  that  some  of  the  observed  cyprids  did  not  settle  immediately  after  molting 
from  a  NVI  larva  and  were  accumulated  over  an  undetermined  number  of  days.  Nauplii 
and  cyprids  are  likely  to  exhibit  different  behaviors,  especially  in  terms  of  their  response 
to  transient  hydrodynamic  features  that  could  transport  them  onshore  (Pineda  1999; 
Helfrich  and  Pineda  2003).  Cyprids  collected  at  the  inshore  station  may  have  also  been 
transported  from  an  adjacent  population,  or  from  an  offshore  source  located  beyond 
Station  3.  Current  velocities  recorded  during  this  study,  and  the  change  in  relative 
importance  of  advective  versus  diffusive  forces  as  a  function  of  distance  to  shore  (see 
Largier  2003),  suggest  that  advection  by  currents  is  more  likely  to  affect  the  along-shore 
distribution  of  larvae.  Thus,  higher  concentrations  of  cyprid  larvae  at  the  inshore  station 
are  probably  the  result  of  an  interaction  between  stage-specific  behavior,  duration  of  the 
cyprid  stage,  and  the  variability  of  cross-shore  flows.  Vertical  distributions  of  cyprids  and 
their  responses  to  transient  hydrodynamic  features  must  be  assessed  in  order  to  further 
investigate  this  question. 

Cross-shore  flows  of  opposite  directions,  dominated  by  the  internal  tide,  were 
often  observed  in  the  current  meter  data,  whereas  along-shore  flows  were  predominately 
uniform  throughout  the  water  column.  Under  such  conditions,  vertically  migrating  larvae 
may  be  able  to  control  their  distance  to  shore,  but  unable  to  control  their  position  on  the 
along-shore  axis,  which  would  be  strongly  affected  by  the  energetic  flows  associated  with 
the  barotropic  tide.  Net  alongshore  transport  may  be  restricted  due  to  the  oscillatory 
nature  of  tidal  flows,  although  it  has  been  shown  that  oscillating  along-shore  flows  can 
generate  complex  spatial-temporal  patterns  of  larval  distribution  (Richards  et  al.  1995). 
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An  increase  in  the  abundance  of  later  stages  at  the  offshore  station  towards  the  end  of  the 
survey  could  be  ascribed  to  an  accumulation  of  larvae  produced  at  a  number  of 
neighboring  populations,  rather  than  to  the  retention  of  larvae  released  only  at  the  study 
area.  No  conclusion  can  be  reached,  however,  without  knowledge  on  the  ontogenetic 
patterns  in  vertical  distribution  and  migration  behavior  of  these  larvae. 

Another  factor  that  could  contribute  to  the  observed  patterns  of  cross-shore  larval 
distribution  is  the  effect  of  wind  forcing  on  the  transport  of  surface  materials.  As 
previously  observed  at  a  similar  site  in  the  region  (Tapia  et  al.  2004),  cross-shore  winds 
blew  onshore  during  the  day  and  slightly  reversed  their  direction  at  night.  Although  at 
this  point  it  is  not  known  whether  late  larval  stages  occur  close  enough  to  the  surface  to 
be  affected  by  onshore  winds,  larval  settlement  observed  at  the  appropriate  frequencies 
could  provide  a  means  to  test  for  this  potential  association.  The  frequency  of  settlement 
observations  in  this  study  did  not  allow  us  to  test  for  a  correlation  with  diurnal  changes  in 
wind  forcing.  Had  settlement  been  monitored  at  a  semi-diurnal  frequency  (e.g.  at  dawn 
and  dusk),  significantly  higher  numbers  of  settlers  should  have  been  observed  at  dusk  if 
cross-shore  winds  have  any  effect  on  the  distribution  of  competent  larvae. 


2.5.2.  Assessing  larval  dispersal  in  nearshore  waters 

As  larvae  develop  in  the  plankton,  they  acquire  behaviors  that  could  either 
increase  dispersal  or  facilitate  retention.  For  example,  crustacean  larvae  use  diurnal 
vertical  migrations  or  other  vertical  swimming  behaviors  to  exploit  vertical  shear  in 
horizontal  currents  to  move  in  and  out  of  estuaries  (Forward  1988)  or  bays  (DiBacco  et 
al.  2001),  or  to  remain  close  to  shore  (Sponaugle  et  al.  2002,  p.  349).  Mortality  and 
dilution,  on  the  other  hand,  cause  a  decrease  in  the  abundance  of  larvae  in  nearshore 
waters  to  an  extent  that  is  currently  unknown  for  most  species  (but  see  Ellien  et  al.  2004). 
Estimating  these  rates  in  the  field  requires  tracking  aggregations  of  larvae  with  a  spatial 
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and  temporal  resolution  that  is  rarely  achieved  (e.g.  Natunewicz  et  al.  2001;  Arnold  et  al. 
2005). 


Tagging  larvae  to  estimate  mortality  and  dispersal  through  a  mark-recapture 
approach  has  numerous  disadvantages  and  logistical  limitations  (see  Levin  1990), 
especially  for  larvae  that  are  as  small  and  abundant  as  the  larvae  of  intertidal  barnacles. 
Natural  or  environmentally-induced  tags  have  provided  a  means  to  assess  larval  dispersal 
over  temporal  scales  of  days  to  months  (e.g.  Swearer  et  al.  1999;  DiBacco  and  Levin 
2000;  Becker  et  al.  2005).  However,  the  applicability  of  this  approach  to  tracking 
barnacle  larvae  is  very  limited,  given  the  duration  of  their  planktonic  life  and  the  loss  of 
hard  structures  during  ecdysis.  Nauplius  I  larvae,  for  instance,  molt  into  Nil  within  hours 
of  hatching  (Brown  and  Roughgarden  1985;  Miller  et  al.  1989).  Therefore,  in  order  to 
assess  the  dispersal  of  barnacle  larvae  over  temporal  scales  of  days  to  weeks  we  must 
resort  to  a  high-frequency  characterization  of  nearshore  circulation  patterns  and  plankton 
distribution. 

Although  more  spatial  coverage  is  required  to  characterize  horizontal  flows  and 
their  variability  in  the  study  area,  it  was  assumed  that  horizontal  current  velocities 
measured  at  the  study  site  were  representative  of  nearshore  conditions.  This  assumption 
is  critical  when  interpreting  results  of  the  PVD  analysis,  which  pointed  to  a  limited 
dispersal  (1-10  km)  of  larval  barnacles  in  nearshore  waters  off  Southern  California. 
Before  these  results  can  be  extrapolated  to  other  sites  in  the  region,  the  spatial  variability 
of  nearshore  flows  must  be  assessed.  It  is  also  necessary  to  investigate  the  role  of  two 
submarine  canyons  flanking  the  sampling  area  (Fig.  2.1)  on  the  retention  ancLor 
aggregation  of  larvae.  Furthermore,  hypotheses  regarding  vertical  migration  behavior, 
ontogenetic  changes  in  the  range  of  vertical  distribution,  and  swimming  responses  to 
transient  hydrodynamic  features,  remain  to  be  tested. 
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2.5.3 .  Mortality  estimates 


Using  vertical  life  tables  to  estimate  mortality  at  every  transition  between  larval 
stages  yielded  a  large  number  of  estimates  that  were  negative,  zero,  or  greater  than  one 
(results  not  shown).  Such  values  are  clearly  outside  the  range  of  values  expected  for  a 
closed  system,  or  for  a  stage  distribution  that  is  representative  of  the  larval  population. 
Although  it  would  not  be  correct  to  use  these  values,  for  instance,  to  parameterize  a 
population  model,  it  is  possible  to  use  them  to  infer  how  open  the  area  was  in  terms  of 
larval  dispersal  over  the  duration  of  our  study. 

On  four  out  of  seven  days,  mortality  estimates  for  the  NII-NIII  transition  in 
Balanus  glandula  were  greater  than  one,  suggesting  that  an  input  of  Nil  and/or  a  loss  of 
NIII  larvae  had  occurred.  Estimates  for  the  NIII-NIV  transition  in  B.  glandula  were  either 
negative  or  zero  on  four  out  of  seven  days,  suggesting  a  loss  of  NIII,  or  an  input  of  NIV 
larvae.  Negative  or  zero  estimates  for  the  transitions  NII-NIII  (2  out  of  7)  and  NIII-NIV 
(5  out  of  7)  in  Chthamalus  spp.  also  suggested  a  loss  of  NIII,  and  probably  of  Nil,  during 
the  first  two  days.  This  is  probably  related  to  ontogenetic  differences  in  swimming 
abilities  and  vertical  distribution.  Perhaps  NII-NIII  larvae  are  less  able  to  determine  their 
position  within  the  water  column,  becoming  easily  entrained  in  along-shore  flows. 

The  estimated  rates  of  larval  mortality  presented  here  were  clearly  affected  by 
insufficient  spatial  coverage  and  resolution  of  our  sampling.  Although  the  vertical  life 
table  method  does  not  require  tracking  an  individual  larval  aggregation,  instantaneous 
stage  structures  utilized  in  the  estimation  of  mortality  are  assumed  to  represent  the 
population’s  stage  composition  (Aksnes  and  Ohman  1996).  Daily  stage  distributions 
observed  in  this  study  were  highly  variable,  and  suggested  that  the  plankton  survey 
lacked  the  spatial  coverage  needed  to  compensate  for  small-scale  patchiness  and 
advection  of  early  naupliar  stages.  Furthermore,  the  assumption  of  constant  and  equal 
development  time  for  all  larvae  in  a  given  stage  (Aksnes  and  Ohman  1996)  could  be 
unrealistic  for  the  larvae  of  B.  glandula  and  Chthamalus  spp.,  which  have  shown  an 
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increase  in  the  variability  of  stage  duration  between  early  and  late  naupliar  stages  (Brown 
and  Roughgarden  1985;  Miller  et  al.  1989). 

The  results  suggest  that  average  mortality  rates  of  barnacle  nauplii  in  nearshore 
waters  fluctuate  around  20-40%  per  day  (Fig.  2.8).  These  values  are  substantially  higher 
than  previously  assumed  for  B.  glanclula  and  Chthamalus  spp.  (Connolly  and 
Roughgarden  1998,  p.  323),  but  within  the  range  of  mortality  estimates  found  in  the 
literature  for  other  benthic  invertebrate  larvae  (Rumrill  1990;  Morgan  1995;  Lamare  and 
Barker  1999).  A  nominal  mortality  rate  of  0.05  d'1  has  been  used  repeatedly  in  modeling 
studies  that  describe  the  distribution  and  population  dynamies  of  B.  glanclula  and 
Chthamalus  spp.  from  Northern  California  and  Oregon  (see  Gaylord  and  Gaines  2000,  p. 
774).  This  estimate  was  reportedly  obtained  from  stage-specific  counts  of  Semibalanus 
balanoides  larvae,  collected  on  3  occasions  over  a  period  of  20  days  from  a  pier  in 
Millport,  Scotland  (Pyefinch  1949).  The  characteristics  of  my  sampling  design  and  target 
species  make  the  data  presented  in  this  contribution  more  likely  to  provide  a  realistic 
estimate  for  the  mortality  of  B.  glandula  and  Chthamalus  spp.  larvae  in  nearshore  waters. 
However,  similar  surveys  must  be  conducted  at  other  sites  in  the  region  to  contrast  and 
validate  these  estimates. 

A  four-  to  six-fold  increase  in  mortality  rates  would  have  a  substantial  effect  on 
the  number  of  larvae  completing  their  planktonic  development.  When  constant  mortality 
rates  of  20  -  40%  d'1  are  used  together  with  a  simple  exponential  decay  function  to 
project  the  numbers  of  larvae  in  a  closed  system  after  2.5  weeks  (average  development 
time  for  barnacle  larvae  at  ~15  °C),  the  numbers  of  larvae  expected  to  complete  their 
development  are  15  -  450  times  smaller  than  the  number  expected  with  a  5%  d'1 
mortality.  Depending  on  larval  duration,  such  increases  in  mortality  rates  could  change 
the  relevant  larval  transport  mechanisms,  shorten  the  mean  expected  travel  distances  (e.g. 
Ellien  et  al.  2004),  and  ultimately  affect  the  scales  of  connectivity  among  populations. 

Results  presented  in  this  paper  have  implications  for  the  current  view  on 
population  openness  and  connectivity  in  coastal  marine  invertebrates.  If  the  high 
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mortality  rates  and  potential  for  limited  larval  dispersal  inferred  from  these  data  are 
characteristic  of  intertidal  populations  in  the  region,  the  distances  at  which  populations 
are  connected  by  larval  exchange  could  be  no  longer  than  1-10  km.  If  populations  with  a 
similar  life  cycle  are  indeed  disconnected  at  such  scales,  future  efforts  to  model  the 
relationship  between  local  population  dynamics  and  hydrodynamics  must  shift  from  an 
emphasis  on  mesoscale  processes  to  a  better  description  of  nearshore  processes  and 
spatial  changes  in  dispersal  and  self-recruitment. 
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2.8.  Tables 
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Table  2.1 


Results  of  two-way  ANOVAs  used  to  test  for  spatial  (between-station)  differences  in  the 
concentration  of  each  larval  stage.  Different  days  of  the  study  were  treated  as  repeat 
observations  (df=  2,6  for  all  tests).  Significant  between-station  differences  are  indicated 
in  bold.  The  results  of  Tukey  pair-wise  comparisons,  performed  for  each  test  with  a 
significant  between-station  effect,  are  shown  on  the  third  column  for  each  species  (  : 
P<0.05;  ”  :  P<0.01). 


Stage 

Balatuis  glandula 

Chthamalus  spp. 

F  2.6 

P 

Pairwise  differences 

F  2.6 

p 

Pairwise  differences 

Nil 

10.973 

0.002 

Stl  >  St3*,  St2>St3” 

3.120 

0.089 

NIII 

3.168 

0.082 

6.508 

0.015 

St2  >  St  f 

NIV 

1.254 

0.323 

25.373 

0.000 

st3  >  st  r,  st2>st  r 

NV 

1.872 

0.200 

32.301 

0.000 

St3  >  St  1*\  Sl2  >  St  1** 

NVI 

3.533 

0.065 

28.256 

0.000 

st3  >  st  r*,  st2  >  st  r 

Cyprids 

2.451 

0.132 

10.326 

0.004 

St  1  >  St  3” 
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Table  2.2 


Estimates  of  mortality  obtained  for  each  day  and  stage  transition.  Empty  cells  correspond 
to  days/transitions  for  which  the  estimated  mortality  rate  was  zero  or  negative. 


MORTALITY  ESTIMATES  (d^l 


Transition 

Day  I 

Day  2 

Day  3 

Day  4 

Day  5 

Day  6 

Day  7 

Balanus  glandula 

II+1II  ->  IV+V 

0.279 

0.043 

0.128 

0.421 

0.288 

0.322 

0.603 

IV+V  VI 

0.393 

0.282 

0.366 

0.036 

0.311 

0.108 

- 

VI  ->  c 

0.123 

0.629 

0.558 

0.210 

0.462 

- 

0.693 

Chthamalus  spp. 

II+1I1  -»  IV+V 

0.239 

- 

0.083 

0.119 

0.243 

0.320 

IV+V  ->  VI 

0.441 

0.224 

0.300 

0.475 

0.315 

0.149 

0.262 

VI  ^  c 

0.082 

0.501 

0.135 

0.015 

0.086 

0.327 

0.089 
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2.9.  Figures 
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32°52' 


117°15' 


32°51' 


1 17°16* 


Figure  2.1 

Study  area  in  La  Jolla,  Southern  California.  The  numbered  crosses  indicate  the 
sampling  stations.  The  letters  A  and  T  show  the  location  of  the  Doppler  current 
meter  (32°51 ,833’N,  1 17°15.786’W)  and  temperature  mooring  (32°51.827’N, 
1 17°  15. 894’ W).  Contour  lines  indicate  depth  in  meters. 
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Temporal  variability  in  the  stage-specific  spatial  distributions  of  Balanus  glandula  larvae. 


/ — 7 — T — t 

/ 

lr  ; - i - (— 

/  /OTjO 

/  '  V  #  "X 

\  \  *.  ■  . 

v  \  -HFfJ 

\  \  V  /  ¥ 

VA  v 

\  \  \ 

\  \  V 1 

\  \  \  \  i 

\  \ 

) - 7 — 

-1 - 1 - 7 - 7 

\  \ 

v  A 

A 

\ 

VT - S ^ V 

f  I  }  i 

f 

/ 

'  1 

!  [Ij 

/ 

?8  8  ?  ° 

t—  r- 

£LU  Jed  0BAJB~J 


o  eo  co  v  <v  o 

£U1  JBd  QBAJB~l 


r*-> 

ci 

a> 

u 


P 


50 


Temporal  variability  in  the  stage-specific  spatial  distributions  of  Chthamalas  spp.  larvae. 


Relative  frequency 


0.75 

0.50 

0.25 

0.75 

0.50 

0.25 

0.75 

0.50 

0.25 

0.75 

0.50 

0.25 

0.75 

0.50 

0.25 

0.75 

0.50 

0.25 

0.75 

0.50 

0.25 


B 


I  ir 


Day  1 
N=  358 


i E 


1 


Day  2 
N=  297 


IL 


Day  3 
A/=219 


Day  4 
AM  37 


“LL 


Day  5 
N-225 


Day  6 
AM  08 


Day  7 
A/=88 


K 


IV  V  VI  C  II 

Larval  stage 


Day  1 
A/=6  83 


Day  2 
A/=  379 


Day  3 
A/-93 


Day  4 
A/— 399 


Day  5 
AM  62 


IV  V  VI  c 


Figure  2.4 

Stage  distributions  of  Balanns  glandula  (A-G)  and  Chthamalus  spp.  (H-N)  larvae 
for  each  day  of  the  survey  (Day  1=6  March  2003).  Relative  stage  frequencies 
were  computed  from  total  counts  of  larvae  (N)  at  the  three  sampling  stations. 
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Day  1  Day  2  Day  3  Day  4  Day  5  Day  6  Day  7  Day  8 


16 

B 

12 

c 

• 

8 

o 

4 

• 

O 

o 

0 

o 

-4 

• 

• 

o 

-8 

-12 

• 

Pearson's  r=  0.955 
P=  0.011 

-L - 1 - L 

-16 

-2.0  -1.5  -1.0  -0.5  0.0  0.5  1.0  1.5  2.0  -8  -6  -4  -2  0  2  4  6  8 

D  Cyprid  abundance  (#  m  3)  D  Cyprid  abundance  (#  m"3) 


Figure  2.5 

Larval  settlement  at  Dike  Rock  (A),  and  the  correlation  between  daily  changes  in 
settlement  and  daily  changes  in  nearshore  cyprid  concentration  for  B.  glandula  (B) 
and  Chthamalus  spp.  (C). 
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Figure  2.6 

Time  series  of  environmental  variables  measured  at  the  study  area.  Along-shore  and 
cross-shore  winds  (A)  were  computed  from  hourly  measurements  gathered  at  the  Scripps 
Institution  of  Oceanography  pier  for  the  month  of  March  2003  (inset).  Box  in  inset 
corresponds  to  A.  The  series  of  along-shore  and  cross-shore  horizontal  currents  (B  and  C, 
both  in  cm  s‘l),  and  temperature  distribution  (D),  were  plotted  using  low-pass  filtered 
series  of  10-point  (i.e.  10  min)  running  means. 
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Along-shore  (cm/s)  Cross-shore  (cm/s)  Temp  (°C) 
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Figure  2.7 

Progressive  Vector  Diagrams  (PVD)  computed  from  current  meter  data  collected 
during  this  study.  Results  are  presented  in  three  panels  for  clarity.  All  trajectories 
(lines)  started  from  Station  1  (square),  and  were  computed  for  six  cases  of  larval 
vertical  distribution  (see  Materials  and  Methods).  Numbers  at  the  end  of  each 
trajectory  correspond  to  the  approximate  time  (in  hours)  after  which  a  particle 
would  have  collided  with  the  shore.  Symbols  along  the  shore  correspond  to  the 
endpoints  of  PVDs  computed  at  different  times  during  a  tidal  cycle.  Circles  and 
diamonds  correspond  to  the  thin  and  thick  lines,  respectively.  The  Scripps 
Institution  of  Oceanography  pier  is  indicated  by  P. 
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Figure  2.8 

Average  estimates  of  larval  mortality  rate  (±  SE-SD Nn)  for  three  stage 
transitions  in  Balanus  glandula  and  Chthamalus  spp.  Negative  and  zero 
values  (see  Table  2.2)  were  not  included  in  the  computations. 
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3.  CHAPTER  3 

STAGE-SPECIFIC  VERTICAL  DISTRIBUTIONS  OF  BARNACLE 

LARVAE  IN  THE  NEARSHORE 

3.1.  Abstract 

The  vertical  distribution  of  barnacle  larvae  was  monitored  at  a  nearshore  station 
off  Southern  California  over  a  48-hour  period  in  June  2003.  Counts  of  nauplii  and  cyprids 
were  tested  for  the  existence  of  between-stage  differences,  as  well  as  diurnal  changes  in 
the  vertical  distribution  of  these  stages.  Plankton  samples  were  collected  every  2  hours 
from  3  depth  intervals  (0.5-9,  11-19,  and  21-29  m).  Environmental  information  was 
collected  from  CTD  casts  conducted  immediately  after  each  plankton  station.  Collected 
larvae  were  enumerated  and  identified  to  species  and  stage  of  development.  Nauplii  of 
the  intertidal  barnacle  Chthamalus  spp.  were  numerically  dominant  (81.7%)  and  present 
at  all  stages  of  larval  development,  yet  cyprids  of  this  species  were  absent  from  our 
samples.  Most  collected  cyprids  (84.2%)  were  identified  as  the  subtidal  barnacle  Balanus 
nubilus ,  which  only  accounted  for  5.6%  of  total  nauplii.  Distinct  patterns  of  distribution 
were  found  when  species-specific  distributions  of  nauplii  and  cyprids  were  compared. 
Naupliar  stages  of  all  species  occurred  near  the  surface,  whereas  cyprids  were 
concentrated  within  the  bottom  half  of  the  water  column  at  all  times.  These  vertical 
patterns,  which  did  not  change  between  day  and  night  or  across  species,  are  considered  in 
relation  to  previously  observed  differences  in  the  horizontal  distribution  of  nauplii  and 
cyprids.  I  review  the  factors  that  are  likely  to  explain  this  vertical  separation  and  their 
implications  for  the  dispersal  of  larval  stages  in  the  nearshore. 
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3.2.  Introduction 


The  life  cycle  of  many  benthic  invertebrates  includes  a  dispersive  larval  stage 
(Thorson  1950).  Sedentary  adults  produce  larvae  that  spend  from  hours  to  months  in  the 
plankton  before  they  become  competent  to  settle  on  adult  habitats  (Levin  and  Bridges 
1995).  During  that  time  they  can  disperse  over  potentially  long  distances  (Scheltema 
1968,  1986),  and  larvae  that  survive  to  settle  often  do  so  in  habitats  that  are  far  from  their 
parental  populations.  Larval  dispersal  is  thus  the  mechanism  by  which  many  disjunct 
populations  of  benthic  invertebrates  are  demographically  and  genetically  connected.  The 
spatial  scale  of  such  exchange  of  individuals  is  relevant  for  understanding  recruitment 
variability  (Sponaugle  et  al.  2002)  and  for  the  design  of  marine  reserves  (Gaines  et  al. 
2003).  Dispersal  of  planktonic  larvae  also  allows  sedentary  species  to  re-colonize  and 
invade  coastal  habitats.  Therefore,  the  characterization  of  factors  that  control  the 
horizontal  distribution  and  dispersal  of  larvae  is  key  to  understanding  species  geographic 
ranges  and  rates  of  invasion  in  coastal  environments. 

Swimming  abilities  of  most  invertebrate  larvae  are  limited  relative  to  the 
horizontal  velocities  of  coastal  currents.  Swimming  speeds  range  from  less  than  1  cm  s'1 
in  ciliated  larvae  to  a  maximum  of  2-8  cm  s’1  in  the  larvae  of  some  crustaceans  (Chia  et 
al.  1984),  whereas  typical  velocities  of  alongshore  currents  may  fluctuate  between  .10  and 
30  cm  s'1  (Largier  2003  and  references  therein).  However,  the  majority  of  larvae  do  not 
behave  as  passive  particles,  and  dispersal  cannot  be  predicted  simply  from  mean  flows 
without  information  on  their  behavior  and  vertical  distribution.  Vertical  shear  in 
horizontal  currents  can  generate  large  differences  in  the  dispersal  patterns  of  larvae 
occurring  at  different  depths  in  the  water  column  (Rothlisberg  et  al.  1995;  Hannah  et  al. 
1997;  Hill  1998).  For  instance,  it  has  been  suggested  that  a  number  of  coastal  and 
estuarine  larvae  can  control  their  horizontal  displacements  by  moving  up  and  down  the 
water  column  at  various  times  of  the  tidal  or  diurnal  cycle  (e.g.  Bousfield  1955;  DiBacco 
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et  al.  2001).  Consequently,  knowledge  of  the  vertical  distribution  of  invertebrate  larvae  is 
essential  to  understand  patterns  of  dispersal  and  horizontal  distribution. 

Intertidal  and  subtidal  barnacles  are  conspicuous  members  of  coastal  communities 
around  the  world.  Their  high  abundance,  wide  geographic  distribution,  and  accessibility 
make  them  ideal  organisms  to  study  the  interaction  between  nearshore  physical  processes 
and  larval  dispersal.  Extensive  information  is  available  on  the  reproductive  cycle, 
ecology  and  larval  development  of  a  number  of  intertidal  species  (e.g.  Barnes  and  Barnes 
1956;  Crisp  1962;  Hines  1979;  Crisp  and  Bourget  1985;  Barnes  1989;  Crisp  et  al.  1991; 
Barnes  1999;  Wethey  2002).  Most  adult  barnacles  are  cross-fertilizing  hermaphrodites 
that  maintain  fertilized  eggs  in  the  mantle  cavity  for  several  weeks  and  release  up  to  tens 
of  thousands  of  nauplius  larvae  into  the  plankton  (e.g.  Hines  1976).  Most  nauplii  must 
feed  on  microalgae  and  undergo  six  molts  before  attaining  a  cyprid  stage,  which  is  the 
non-feeding  larva  that  eventually  settles  onto  hard  substrata  and  metamorphoses  into  a 
juvenile  barnacle. 

Recent  observations  in  Southern  California  show  substantial  differences  in  the 
horizontal  distribution  of  cyprids  and  the  naupliar  stages  that  precede  them  (see  Chapter 
2),  suggesting  that  nauplii  and  cyprids  are  affected  by  different  sets  of  advective  forces  in 
the  nearshore.  Differences  in  the  vertical  distribution  of  nauplii  and  cyprids  may  have 
contributed  to  the  observed  horizontal  segregation  of  barnacle  larval  stages.  Here  I 
document  temporal  patterns  in  the  vertical  distribution  of  barnacle  larvae  observed  in 
nearshore  waters  off  Southern  California  over  a  48-hour  period.  The  objective  of  this 
study  was  to  test  for  the  existence  of  (1)  between-stage  differences  in  vertical  distribution 
and  (2)  diurnal  changes  in  the  vertical  distribution  of  these  stages.  The  results  indicate 
that  nauplii  and  cyprids  of  the  subtidal  Balanus  nubilus  and  the  intertidal  gooseneck 
barnacle  Pollicipes  polymerus  occur  at  different  depths,  and  that  their  vertical 
distributions  do  not  change  between  day  and  night. 
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3.3.  Materials  and  Methods 


Sampling  was  conducted  in  June  2003  onboard  the  R/V  Robert  Gordon  Sproul  of 
the  Scripps  Institution  of  Oceanography  (SIO).  The  ship  was  anchored  for  48  hours, 
starting  on  8  June  at  23:00  Pacific  Daylight  Time  (PDT),  at  a  coastal  station  located  at  a 
depth  of  30  m  and  approximately  2  km  off  Del  Mar  (32°55.9I’N,  117°16.79’W)  in 
Southern  California,  USA  (Fig.  3.1).  Plankton  samples  were  collected  every  2  hours  from 
3  depth  intervals  (0.5-9,  11-19,  and  21-29  meters)  using  a  semi-vortex  AC  pump  (Ebara 
DVS,  Ebara  Corp.).  To  insure  an  even  sampling  of  water,  the  pump  was  hauled  up  and 
down  within  each  depth  interval  at  ca.  1  m  s"1  for  approximately  20  min,  until  6  m3  of 
water  were  filtered  through  a  1 10  p.m  mesh  net.  Volume  of  water  filtered  was  measured 
with  a  digital  flow  meter  attached  to  the  pump  hose.  Plankton  samples  were  preserved 
onboard  using  buffered  4%  formaldehyde,  and  sorted  under  a  dissecting  microscope  in 
the  laboratory. 

A  CTD-Rosette  (Seabird  SBE9)  equipped  with  a  fluorometer  (Seatech)  was 
lowered  down  to  29  m  immediately  after  a  plankton  station  was  completed.  Additional 
environmental  data  was  obtained  from  the  SIO  Coastal  Data  Information  Program 
website  (http://cdip.ucsd.edu).  Hourly  data  recorded  at  the  SIO  pier  (32°52.00’N, 
1 17°15.42’W,  ca.  8  km  SE  of  the  ship’s  position)  were  used  to  compute  alongshore  and 
cross-shore  components  of  wind.  Winds  were  rotated  in  order  to  align  them  with  the 
shoreline  orientation  in  the  area  (350°  relative  to  true  North).  Data  on  wave  significant 
height  and  direction  recorded  every  30  min  at  a  nearby  buoy  (32°57.72’N,  1 17°16.38’W, 
Fig.  3.1)  were  obtained  from  the  same  database. 

All  barnacle  larvae  in  the  samples  were  enumerated.  Nauplii  and  cyprids  were 
identified  to  the  lowest  possible  taxonomic  level.  The  species  identification  and  staging 
of  nauplius  larvae  were  based  on  morphological  criteria  given  by  Lewis  (1975),  Standing 
(1980),  Brown  and  Roughgarden  (1985),  Miller  et  al.  (1989),  and  Miller  and 
Roughgarden  (1994).  Chthamalid  larvae  encountered  in  our  samples  were  probably  those 
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of  Chtharnalus  fissus,  a  species  that  occurs  between  San  Francisco  and  Baja  California 
(Newman  and  Abbott  1980)  and  that  numerically  dominates  the  high  to  upper  middle 
intertidal  in  La  Jolla.  However,  I  refer  to  the  larvae  in  our  samples  as  Chtharnalus  spp. 
because  the  geographic  distribution  of  C.  fissus  occasionally  overlaps  with  that  of  C.  dalli, 
a  northern  species  with  a  southern  limit  that  fluctuates  between  Point  Conception  and  San 
Diego  (Newman  and  Abbott  1980;  Wares  2001).  Adult  C.  fissus  and  C.  dalli  can  be 
identified  upon  dissection  and  microscopic  examination,  but  their  larvae  arc 
morphologically  identical  (Miller  et  al.  1989). 

Counts  of  nauplii  and  cyprids  for  each  species  were  used  to  test  for  between-stage 
differences  in  vertical  distribution,  as  well  as  day/night  differences  within  stages.  1  used  a 
statistical  test  that  allows  for  comparisons  of  vertical  plankton  distributions  in  the 
presence  of  patchiness  and  when  multiple  replicates  are  available  (Beet  et  al.  2003). 
Samples  taken  at  different  times  were  pooled  into  two  sets  of  observations,  daytime  and 
nighttime,  and  considered  as  replicates.  Three  plankton  profiles  collected  at  transition 
times  (i.e.  dawn  or  dusk)  were  not  used  in  the  analysis. 

I  tested  the  null  hypothesis  that  means  of  larval  abundance  at  each  depth  under 
one  condition  (i.e.  daytime)  are  the  same  as  the  means  under  another  condition  (i.e. 
nighttime),  using  the  test  statistic 


i=l  /=!  Mi 


T-TT) 


(3.1) 


where  Y{J  is  the  average  count  of  n  replicates  for  condition  i  and  depth  j,  and  jutj  and  c  are 

maximum  likelihood  estimates  of  the  mean  (□/,)  and  dispersion  coefficient  (cf  of  a 
negative  binomial  distribution,  which  is  recommended  to  describe  count  data  with  a 
variance  to  mean  ratio  greater  than  one  (Hayek  and  Buzas  1997).  The  log-likelihood 
function  and  details  on  the  procedure  used  to  obtain  the  maximum  likelihood  estimates 
are  given  by  Beet  et  al.  (2003). 
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To  correct  for  potential  effects  of  serial  autocorrelation  within  groups  (i.e.  day  and 
night),  statistical  comparisons  of  vertical  distributions  were  performed  using  subsets  of 
plankton  profiles  selected  at  random  from  the  pooled  datasets.  This  procedure  was 
repeated  ten  times  for  each  comparison,  and  ranges  of  values  obtained  for  the  test  statistic 
B  and  its  corresponding  probability  under  the  null  hypothesis  (i.e.  p-values)  are  reported 
below.  The  same  sub-sampling  scheme  was  used  when  the  vertical  distributions  of 
nauplii  and  cyprids  were  compared. 

Finally,  temporal  changes  in  the  vertical  distribution  of  nauplii  and  cyprids  were 
compared  with  those  observed  for  temperature  and  chlorophyll-a  in  the  water  column 
using  the  Mean  Depth  of  Distribution  (MDD)  of  each  larval  group.  MDDs  were 
calculated  for  each  plankton  profile  j  as  the  weighted  average 

(3.2) 

where  z,  is  the  mean  depth  of  interval  /,  r\y{  is  the  number  of  nauplii  or  cyprids  collected 
within  that  interval,  and  Nj  is  the  total  number  of  nauplii  or  cyprids  found  at  time  j . 
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3.4.  Results 


Total  larval  counts  showed  large  and  persistent  differences  in  the  abundance  of 
nauplii  and  cyprids  throughout  the  sampling  period  (Fig.  3.2).  Nauplii  were  more 
abundant  near  the  surface  (Fig.  3.2a),  with  89%  of  total  nauplii  collected  between  0.5-9  m. 
In  contrast,  most  cyprids  were  found  away  from  the  surface  (Fig.  3.2b),  with  34%  and 
59%  of  the  total  counts  collected  within  the  mid-depth  and  near-bottom  intervals, 
respectively. 

Chthamalus  spp.  represented  78%  of  the  total  nauplii  found  at  any  given  depth 
interval  (Fig.  3.3).  Nauplii  of  the  intertidal  gooseneck  barnacle  Pollicipes  polymerus  and 
the  subtidal  species  Balanus  mibilus  occurred  with  lower  abundance  and  were  mostly 
found  near  the  surface  (Fig.  3.3).  A  contrasting  species  composition  was  observed  in 
cyprids,  with  a  numerical  dominance  of  B.  mibilus  (84.2%  of  total  cyprids)  and  the 
absence  of  Chthamalus  spp.  (Fig.  3.3).  All  naupliar  stages  of  Chthamalus  spp.  were 
observed  at  the  near-surface  depth  interval,  with  a  shift  in  stage  composition  over  the  48 
hours  of  sampling  (Fig.  3.4).  The  stage  structures  of  P.  polymerus  and  B .  mibilus  nauplii, 
on  the  other  hand,  were  dominated  by  early  stages  (i.e.  NII-NIV)  throughout  the 
observation  period  (Fig.  3.4). 

Mean  fractions  of  nauplii  and  cyprids  computed  from  species-specific  counts 
showed  no  day/night  differences  within  each  group,  but  suggested  a  consistent  difference 
in  the  distributions  of  nauplii  and  cyprids  (Fig.  3.5).  This  was  confirmed  by  the  statistical 
comparisons  of  vertical  profiles,  which  showed  that  the  distributions  of  nauplii  or  cyprids 
do  not  significantly  change  between  day  and  night,  but  are  significantly  different  from 
each  other  at  all  times  (Table  3.1).  This  vertical  separation  between  nauplii  and  cyprids 
was  also  apparent  when  total  counts  were  used  to  compute  mean  depths  of  distribution 
(MDD).  At  all  times  during  the  observation  period,  the  majority  of  cyprids  were  found  at 
a  greater  depth  than  nauplii  (Fig.  3.6).  There  was  no  obvious  correspondence  between  the 
observed  temporal  changes  in  MDD  of  nauplii  or  cyprids  and  fluctuations  in  cross-shore 
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or  alongshore  winds  (Fig.  3.6c),  wave  action  or  tidal  height  (data  not  shown). 
Temperature  and  chlorophyll  contours  from  CTD  measurements  showed  that  the  majority 
of  cyprids  remained  below  the  surface  mixed  layer  at  all  times  (Fig.  3.6a),  and  that  the 
MDD  of  naupliar  stages  did  not  overlap  with  depths  at  which  maximum  chlorophyll 
concentrations  were  detected  (Fig.  3.6b). 
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3.5.  Discussion 


Significant  differences  were  found  in  the  vertical  distributions  of  nauplii  and 
cyprids  of  Balanus  nubilus  and  Pollicipes  polymerus.  No  day/night  differences  were 
observed  in  the  distribution  of  sampled  larval  stages  for  B.  nubilus ,  Chthamalus  spp.,  and 
P.  polymerus.  A  majority  of  the  collected  cypnds  (84.2%)  were  those  of  Balanus  nubilus , 
a  subtidal  barnacle  that  occurs  from  the  low  intertidal  to  depths  of  90  m  (Newman  and 
Abbott  1980).  This  species  accounted  for  only  5.6%  of  total  nauplii.  In  contrast,  nauplii 
of  the  intertidal  barnacle  Chthamalus  spp.  were  numerically  dominant  (81.7%)  and 
present  at  all  stages  of  development,  yet  cyprids  of  this  species  were  not  found  in  the 
samples.  This  absence  of  Chthamalus  spp.  cyprids,  despite  the  abundance  of  their  nauplii, 
is  an  intriguing  result  that  I  discuss  below. 


3.5.7.  Larval  distribution  in  relation  to  environmental  variables 

Mean  Depths  of  Distribution  (MDD)  computed  for  total  nauplii  did  not  appear  to 
track  the  vertical  distribution  of  chlorophyll  concentration.  Although  the  observed 
chlorophyll  values  were  rather  low  for  coastal  waters  (e.g.  Barnett  and  Jahn  1987; 
Phillips  2005),  nauplii  were  expected  to  occur  at  depths  with  the  highest  pigment 
concentration  if  their  constant  vertical  distribution  is  to  be  interpreted  as  an  adaptation  for 
feeding.  Transient  increases  in  near-surface  and  mid-depth  chlorophyll  fluorescence 
associated  with  high-frequency  (>  2  cycles  IT1)  internal  waves  have  been  observed  during 
spring  and  summer  in  this  area  (Lennert-Cody  and  Franks  2002).  Given  the  sampling 
frequency  and  short  duration  of  this  study,  such  features  cannot  be  ruled  out  as  a  potential 
source  of  food  for  nauplii  in  near-surface  depths. 

Forcing  of  surface  layers  by  wind  may  have  a  potential  effect  on  the  horizontal 
distribution  of  larvae  (e.g.  Bertness  et  al.  1996;  Garland  et  al.  2002).  The  cross-shore 
component  of  wind  at  our  study  site  was  either  positive  or  just  below  zero,  and  a  signal 
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corresponding  to  the  onshore  sea  breeze  could  be  detected  during  the  daytime  hours  (Fig. 
3.6c).  This  persistent  onshore  breeze  could  have  an  influence  on  the  horizontal 
distribution  of  near-surface  material  (Tapia  et  al.  2004).  However,  its  effect  on  the 
position  of  larvae  with  respect  to  shore  can  only  be  speculated  at  this  point. 


3.5.2 .  Stage-specific  differences  in  distribution 

Differences  in  the  vertical  distributions  of  barnacle  nauplii  and  cyprids  appeared 
consistent  with  observations  previously  reported  for  coastal  and  estuarine  waters.  Over  4 
years  of  sampling  at  a  site  ca.  50  km  north  of  our  sampling  station,  Barnett  and  Jahn 
(1987)  found  that  total  cyprids  (species  not  identified)  were  consistently  more  abundant 
near  the  bottom  and  within  4  km  from  shore.  Pineda  (1991)  found  unidentified  cyprids 
that  were  more  abundant  near  the  bottom  than  at  the  surface  in  shallow  waters  (4-16  m) 
off  La  Jolla,  Southern  California.  Bousfield  (1955)  observed  that  early  nauplii  of  Balanus 
improvisus  and  Balanus  crenatus  in  the  Miramichi  Estuary  occurred  closer  to  the  surface, 
whereas  nauplius  VI  and  cyprids  were  concentrated  near  the  bottom.  At  a  shallow  bay  in 
New  Brunswick,  Canada,  cyprids  of  Semibalamis  balanoides  were  most  abundant  within 
the  bottom  half  of  the  water  column  and  showed  no  apparent  day/night  difference  in  their 
vertical  distribution  (Miron  et  al.  1999,  Fig.  3).  A  contrasting  distribution  was  found  for 
cyprids  of  the  intertidal  barnacle  Verruca  stroemia  at  tens  of  kilometers  off  the  west  coast 
of  Brittany  (Le  Fevre  and  Bourget  1991).  Cyprids  occurred  predominantly  in  the  upper 
meter  of  the  water  column,  whereas  nauplii  had  a  more  uniform  distribution.  This 
suggests  that  differences  in  the  vertical  distribution  of  nauplii  and  cyprids  are  probably 
species  specific,  and  that  vertical  differences  observed  in  this  study  for  B.  nubilus  and  P. 
polymerus  cannot  be  extrapolated  to  Chthamalus  spp.  until  both  nauplii  and  cyprids  are 
observed  over  similar  temporal  scales. 

The  absence  of  Chthamalus  spp.  cyprids,  despite  high  concentrations  of  nauplii, 
could  be  due  to  differences  in  the  behavior  and  swimming  abilities  of  these  stages.  In 
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March  2003,  Tapia  and  Pineda  (submitted,  see  Chapter  2)  observed  consistently  higher 
concentrations  (up  to  22  per  nr)  of  Chthamalus  spp.  cyprids  within  300  m  of  shore,  at  ca. 
10  m  depth,  over  a  7-day  period  in  La  Jolla,  approximately  10  km  SE  of  our  sampling 
station.  Concentrations  declined  to  almost  zero  at  a  station  located  1.1  km  offshore  (ca. 
40  m  depth).  Thus,  it  is  possible  that  Chthamalus  spp.  cyprids  were  not  found  in  this 
study  because  at  a  station  located  ca.  2  km  offshore  their  concentrations  were  too  low  to 
be  detected.  Differences  in  behavior  and  perhaps  swimming  abilities  between 
Chthamalus  spp.  cyprids  and  nauplii  might  explain  these  low  concentrations,  as  cyprids 
could  be  better  able  to  move  shorewards  and  to  aggregate  near  adult  habitats.  Cyprids 
might  accomplish  this  by  responding  to  cross-shore  moving  hydrodynamic  features 
(Pineda  1999;  Helfrich  and  Pineda  2003). 

The  above  argument  must  be  qualified  by  the  fact  that  different  sampling  gear 
was  utilized  in  the  present  study  and  that  conducted  by  Tapia  and  Pineda  (see  Chapter  2). 
Plankton  samples  taken  by  Tapia  and  Pineda  corresponded  to  vertical  tows  of  a  75  cm 
diameter  net  with  1 10  pm  mesh  size.  Average  volumes  of  water  filtered  ranged  from  8.25 
m  at  the  nearshore  station  to  18.15  m  at  the  station  located  1.1  km  from  shore.  Samples 
taken  for  the  present  study  were  collected  by  filtering  6  m3  of  water  with  a  pump  that  was 
lowered  to  a  maximum  depth  of  29  m,  at  a  site  with  an  average  depth  of  30  m.  It  could  be 
argued  that  Chthamalus  spp.  cyprids  were  not  collected  because  they  occurred  within  the 
space  between  maximum  pumping  depth  and  the  bottom.  However,  Pineda  (1999) 
collected  cyprids  of  Chthamalus  spp.  and  Pollicipes  polymerus  using  the  same  type  of 
pump  to  filter  4  m3  of  water  at  maximum  depths  of  4  m  below  the  surface,  at  least  2  m 
above  the  bottom.  Thus,  it  is  unlikely  that  the  absence  of  Chthamalus  spp.  cyprids  was 
due  to  low  sampling  efficiency.  Still,  confirming  their  absence  from  our  sampling  station 
does  not  prove  that  their  horizontal  distribution  resembled  that  observed  by  Tapia  and 
Pineda  in  La  Jolla  (see  Chapter  2),  and  a  future  study  that  includes  spatial  replication  is 
needed  to  verify  these  patterns. 
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3.5.  3.  Implications  for  larval  horizontal  distribution 


The  occurrence,  forcing  factors,  and  ecological  implications  of  vertical 
movements  in  the  plankton  have  been  the  subject  of  extensive  research  for  the  past  70 
years  (e.g.  Hardy  1936;  Forward  1988).  Diurnal  vertical  migrations,  in  which  planktonic 
organisms  move  between  near-bottom  and  near-surface  waters  at  dawn  and  dusk,  have 
been  interpreted  as  a  predator  avoidance  behavior  (e.g.  Zaret  and  Suffem  1976;  Ohman 
and  Cohen  1983),  as  well  as  a  mechanism  by  which  many  planktonic  taxa  control  their 
position  with  respect  to  shore  (e.g.  Shanks  1986;  Rothlisberg  et  al.  1995;  Poulin  et  al. 
2002).  Numerous  examples  are  available  in  the  literature  where  the  vertical  motions  of 
certain  taxa,  mostly  coastal,  are  such  that  vertical  differences  in  horizontal  flows  allow 
them  to  remain  in,  or  to  move  out  of  a  given  area  (e.g.  Bousfield  1955;  Pineda  1999; 
DiBacco  et  al.  2001).  Ontogenetic  changes  in  vertical  distribution  could  result  from 
changes  in  swimming  and  feeding  behavior.  It  has  been  suggested  that  earlier  stages  may 
spend  more  time  in  near-surface  waters  due  to  a  phototactic  response  enhanced  by 
starvation  (e.g.  Singarajah  et  al.  1967),  whereas  later  stages  may  spend  more  time  in 
deeper  waters  to  avoid  visual  predators  (e.g.  Hays  1995).  Transient  shifts  in  the  vertical 
distribution  of  barnacle  cyprids  can  result  from  a  behavioral  response  to  transient 
hydrodynamic  features  (Pineda  1999).  It  has  been  shown  for  copepods  that  individual 
variability  in  energy  reserves,  measured  as  lipid  storage,  can  explain  intra-specific 
differences  in  the  amplitude  of  vertical  migrations  (Hays  et  al.  2001). 

Why  do  cyprids  remain  close  to  the  bottom,  whereas  nauplii  occur  in  the  upper 
half  of  the  water  column?  The  studies  cited  above  suggest  that  energetic  advantages  and 
behavioral  differences  between  cyprids  and  nauplii  could  be  invoked  to  answer  this 
question.  First,  cyprids  are  non-feeding  larvae,  hence  feeding  cannot  be  considered  as  a 
mechanism  forcing  their  vertical  distribution.  Second,  lower  temperatures  in  bottom 
waters  might  reduce  the  consumption  rate  of  lipid  reserves,  thus  allowing  cyprids  to 
remain  competent  for  longer  periods  of  time.  Laboratory  and  field  observations  indicate 
that  low  temperatures,  or  even  freezing,  can  lengthen  the  competency  period  for  some 
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species  (Pechenik  et  al.  1993;  Pineda  et  al.  2005).  Near-bottom  waters  at  the  sampling 
station  were  ca.  7°C  colder  than  surface  waters  throughout  the  48-h  survey,  and  similar 
vertical  differences  have  been  recorded  at  La  Jolla  through  the  spring  and  summer  (e.g. 
Pineda  and  Lopez  2002,  Fig.  4).  A  near-bottom  habitat  could  thus  increase  the 
competency  time  of  cyprids,  and  therefore  increase  their  chance  to  find  suitable  substrate 
and  settle.  Also,  bottom  waters  may  provide  cyprids  with  better  conditions  to  remain  in 
the  nearshore.  Wind-driven  cross-shelf  transport  of  bottom  waters  (Garland  et  al.  2002) 
and  onshore  moving  features  such  as  internal  tidal  bores  (Pineda  1999)  could  aid  these 
larvae  in  reaching  a  suitable  substrate  to  settle.  Observations  presented  here  on  the 
vertical  distribution  of  cyprids,  and  their  apparent  concentration  in  nearshore  waters 
(Chapter  2)  raise  the  question  of  whether  cyprids  of  intertidal  barnacles  found  tens  of 
kilometers  offshore  (e.g.  Roughgarden  et  al.  1988;  Le  Fevre  and  Bourget  1991)  are  likely 
to  return  to  adult  habitats  and  settle  successfully  during  their  competency  period.  Perhaps 
larvae  found  at  such  distances  from  shore  have  been  lost  from  the  pool  of  potential 
settlers,  and  should  not  be  considered  as  viable  propagules  in  models  that  integrate 
mesoscale  physical  processes  and  the  dynamics  of  benthic  populations  (e.g.  Alexander 
and  Roughgarden  1996). 

In  summary,  I  have  presented  distinct  patterns  of  vertical  distribution  for  different 
stages  of  barnacle  larvae  in  nearshore  waters.  Nauplii  occur  near  the  surface,  whereas  the 
non-feeding  cyprids  are  concentrated  in  the  bottom  half  of  the  water  column.  These 
vertical  patterns,  which  did  not  appear  to  change  between  day  and  night  or  across  species, 
suggest  that  previously  observed  differences  in  the  horizontal  distribution  of  larvae  may 
be  caused  by  ontogenetic  changes  in  vertical  distribution,  and  probably  behavior,  that 
allow  cyprids  to  respond  to  transient  hydrodynamic  features  and  to  move  shorewards 
more  effectively  than  nauplii. 
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3.8.  Tables 


74 


Table  3.1 


Comparisons  of  vertical  distributions  for  species-specific  and  total  counts  of 
barnacle  larvae.  The  test  statistic  B  (Beet  et  al.  2003)  and  its  corresponding  p-value 
were  computed  ten  times  for  each  comparison.  The  range  of  observed  values  is 
reported.  Asterisks  indicate  comparisons  were  significant  differences  were  found 
(a=0.05). 


Comparison 

B 

P 

B  alarms  nub  Hus 

NaupJii  day  vs  Nauplii  night 

Cyprids  day  vs  Cyprids  night 

Nauplii  day  vs  Cyprids  day  * 

Nauplii  night  vs  Cyprids  night  * 

0.0000  -2.9288 

0.0976-  1.6367 

9.6415  -  15.666 

15.589  -21.426 

0.2312  -  1.0000 

0.4412  -0.9524 

0.0004  -  0.0081 

0.0000  -  0.0004 

Pollicipes  polymer  us 

Nauplii  day  vs  Nauplii  night 

Nauplii  day  vs  Cyprids  day  * 

0.4864  -  2.4200 

14.630  -  19.584 

0.2982  -  0.7841 

0.0000  -  0.0007 

Chthamalus  spp. 

Nauplii  day  vs  Nauplii  night 

0.3087  -3.1196 

0.2102  -  0.8570 

Total  counts 

Nauplii  day  vs  Nauplii  night 

Cyprids  day  vs  Cyprids  night 

0.8009  -4.3161 

0.0002  -  1.2561 

0.1156-  0.6700 

0.5336  -  0.9999 
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3.9.  Figures 


76 


Fraction  of  total  count 


N AUPLM 


0  6  12  18  24  30  36  42  48 


C  Y  P  R  I  D  S 


0  6  12  18  24  30  36  42  48 


o 

Ln 

I 

CD 

3 


<£> 

3 


Time  (hours) 


Time  (hours) 


Megabalanus  californicus 
E— j.M  Balanus  nubilus 
i.  -  -  i  Pollicipes  polymerus 
Chthamalus  spp. 


Figure  3.3 

Species  composition  of  total  counts  for  nauplii  (left  panels)  and  cyprids  (right 
panels)  at  three  depth  intervals.  The  four  most  abundant  species  are  shown. 
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Temporal  variability  in  the  stage  composition  of  Chthamalus  spp.,  Pollicipes  polymerus,  and  Balanus  nubilus 
nauplii  at  the  three  depth  intervals  sampled  in  this  study. 


Depth  (m) 


All  larvae  Chthamalus  spp.  B.  nubilus  P.  polymerus 


0.0  0.5  1.0  0.0  0.5  1.0  0.0  0.5  1.0  0.0  0.5  1.0 


0.0  0.5  1.0  0.0  0.5  1.0  0.0  0.5  1.0  0.0  0.5  1.0 


Mean  fraction  of  total  count 


Figure  3.5 

Mean  vertical  distributions  of  barnacle  nauplii  and  cyprids  during  the  day  (white 
bars)  and  night  (black  bars)  at  a  coastal  station  in  Southern  California  (32°55.91’N, 
1 17°16.79’W).  Fractions  of  total  counts  were  computed  separately  for  nauplii  and 
cyprids  in  each  plankton  profile,  and  averaged  for  each  depth  interval  across 
daytime  and  nighttime  samples.  Error  bars  are  standard  errors. 
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Figure  3.6 

Mean  Depths  of  Distribution  (MDD)  of  barnacle  larvae  and  environmental 
variability.  MDDs  are  superimposed  on  temperature  (A)  and  chlorophyll  (B) 
contours  from  CTD  data  collected  after  each  plankton  station.  Hourly  winds  (C) 
were  measured  at  a  site  ca.  8  km  SE  of  the  sampling  station.  Horizontal  black  and 
white  bars  indicate  nighttime  and  daytime.  Dashed  vertical  lines  correspond  to 
sunrise  and  sunset  time. 
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4.  CHAPTER  4 

SPATIAL  PATTERNS  IN  LIFE  HISTORY  TRAITS  OF  INTERTIDAL 
BARNACLES:  IMPLICATIONS  FOR  GEOGRAPHIC  RANGES 

4.1.  Abstract 

Regional  changes  in  demography  of  the  intertidal  barnacle  Balanus  glandula  were 
studied  over  400  km  of  shoreline  in  Southern  California  and  Baja  California.  Significant 
spatial  differences  were  found  in  the  size  distribution,  growth,  and  survival  rates  of  adult 
individuals  along  the  study  region.  Size-specific  fertility  did  not  significantly  change  over 
the  same  geographic  range.  Field  observations  were  used  to  parameterize  a  stage- 
structured  matrix  population  model  for  sites  at  the  extremes  and  at  the  middle  of  the 
region.  Estimates  of  growth  rate  and  net  reproductive  rate  were  greater  than  one 
throughout  the  region,  although  values  calculated  for  the  southernmost  site,  Punta  Baja, 
were  substantially  smaller  and  closer  to  a  demographic  threshold  that  separates 
persistence  from  local  extinction.  Sensitivity  calculations  indicated  that  larval  survival 
had  the  largest  effect  on  population  growth  and  net  reproductive  rate.  Robustness  of  the 
three  local  models  to  environmental  stochasticity  was  assessed  by  simulating  different 
regimes  of  interannual  variability  in  the  survival  of  larvae  and  juveniles.  Simulation 
results  showed  large  between-site  differences,  with  the  model  parameterized  for  Punta 
Baja  being  the  least  robust  to  environmental  stochasticity  and  more  prone  to  local 
extinction.  I  conclude  that  a  spatial  gradient  in  vital  rates  and  temporal  variability  in  the 
environment  would  be  sufficient  conditions  to  generate  a  range  endpoint  for  Balanus 
glandula  near  Punta  Baja. 
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4.2.  Introduction 


Biogeography  is  concerned  with  the  geographic  distribution  of  species  and  the 
factors  that  determine  their  range  endpoints  over  long  periods  of  time.  The  endpoints  of  a 
species’  range  often  correlate  with  physical  and  biological  factors  that  prevent  or  limit 
dispersal  of  individuals  (for  a  review  see  Gaston  2003,  Chapter  2).  Mountains,  lakes,  or 
an  ocean  can  be  effective  barriers  for  the  dispersal  of  terrestrial  species.  In  the  marine 
environment,  barriers  for  dispersal  are  more  difficult  to  characterize.  In  the  absence  of 
physical,  topographic  barriers,  the  exchange  of  individuals  between  two  areas  is 
constrained  by  hydrodynamic  features  and  discontinuities  in  physical-chemical 
conditions  that  often  change  in  their  intensity  and  spatial  extent.  Thus,  the  geographic 
limits  of  distribution  for  marine  species  are  intrinsically  more  flexible  than  those  of 
terrestrial  species. 

Many  coastal  species,  particularly  benthic  invertebrates,  have  complex  life  cycles 
(Thorson  1950).  Sedentary  adults  produce  planktonic  larvae  that  can  disperse  over 
potentially  long  distances  during  their  development  (Schcltema  1986).  Surviving  larvae 
settle  as  adults  on  habitats  that  are  often  far  from  their  parental  populations.  Larval  stages 
are  thus  the  “vector”  by  which  many  benthic  species  can  colonize  new  habitats,  and  the 
potential  limits  of  a  species  range  are  determined  by  physical  and  biological  processes 
that  affect  the  horizontal  distribution  of  larvae.  The  factors  that  ultimately  determine  the 
geographic  range  of  a  species,  however,  are  those  that  affect  post-settlement  survival  and 
the  subsequent  recruitment  of  settlers  into  reproductive  adults  (see  Gosselin  and  Qian 
1997).  Regional  differences  in  these  factors  must  affect  the  spatial  variability  of  a 
population’s  demography  (i.e.  its  structure  and  vital  rates)  and,  by  extension,  the  position 
and  variability  of  its  range  endpoints. 

Range  endpoints  have  been  traditionally  attributed  to  discontinuities  in  the 
environment  and  to  physiological  limitations  of  a  species  (e.g.  Hutchins  1947;  Cox  and 
Moore  2000;  Gaston  2003).  Regional  differences  in  survival,  growth,  and  fertility  of 


84 


benthic  invertebrates  have  been  correlated  with  changes  in  environmental  factors  such  as 
the  availability  of  food  and  temperature  (e.g.  Sanford  and  Menge  2001;  Phillips  2005), 
availability  of  substrate  (e.g.  Wethey  1984),  physiological  stress  and  inter-specific 
competition  (Wethey  2002),  as  well  as  regional  differences  in  larval  supply  (e.g. 
Connolly  and  Roughgarden  1998).  Recent  analyses  of  metapopulation  models,  however, 
have  shown  that  range  endpoints  could  arise  from  smooth  spatial  gradients  in 
demographic  rates,  in  the  absence  of  environmental  discontinuities  (Lennon  et  al.  1997; 
Holt  and  Keitt  2000).  Other  modeling  studies  have  suggested  that  dynamic  physical 
boundaries  such  as  coastal  currents  could  set  and  maintain  limits  in  the  distribution  of 
coastal  species  (Gaylord  and  Gaines  2000)  in  the  absence  of  demographic  gradients  or 
physiological  limitations. 

In  this  paper  I  combine  field  observations  and  a  stage-structured  matrix 
population  model  to  submit  that  the  southern  range  endpoint  of  Balanus  glandula ,  an 
intertidal  barnacle  from  the  west  coast  of  North  America,  could  be  explained  by  spatial 
changes  in  demographic  rates  that  make  local  populations  less  robust  to  environmental 
stochasticity.  Matrix  models  have  been  extensively  applied  to  terrestrial  and  marine 
populations  (Caswell  et  al.  1997).  Age-structured  and  size-structured  models  have  been 
previously  formulated  for  sedentary  benthic  species  with  complex  life  cycles 
(Roughgarden  et  al.  1985;  Pascual  and  Caswell  1991;  Hyder  et  al.  2001;  Svensson  et  al. 
2004).  I  use  a  simple  two-stage  model  to  explore  the  effects  of  environmental 
stochasticity  and  spatial  changes  in  vital  rates  on  the  persistence  of  local  populations. 
Parameters  for  this  model  were  estimated  from  field  observations  on  size  structure, 
growth,  survival,  and  reproductive  activity  of  adult  B.  glandula  over  1.5  years  in  a  400- 
kilometer  region  of  Southern  California  and  Baja  California.  I  conclude  that  a  spatial 
gradient  in  vital  rates  and  temporal  variability  in  the  environment  would  be  sufficient 
conditions  to  generate  a  range  endpoint  for  B.  glandula  near  the  southern  limit  of  our 
study  region. 
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4.3.  Materials  and  Methods 


4.3.  L  Species ,  study  sites  and  sampling  procedure 

The  acorn  barnacle  Balanus  glandula  Darwin,  1854  is  a  widespread  and  abundant 
intertidal  species  that  occurs  in  the  middle  to  high  intertidal  zone  along  the  west  coast  of 
North  America,  from  the  Aleutian  Islands  to  Punta  Baja,  Baja  California  (Newman  and 
Abbott  1980  and  author’s  personal  observations).  Recent  reports  indicate  its  occurrence 
as  an  invasive  species  of  intertidal  habitats  in  Japanese  (Kado  2003)  and  Argentinean 
shores  (Orensanz  et  al.  2002).  The  reproduction  of  B.  glandula  is  seasonal  in  Southern 
California  and  Baja  California.  Larvae  are  released  by  adults  from  early  winter  to  late 
spring  (Barnes  and  Barnes  1956;  Hines  1979;  Newman  and  Abbott  1980),  and  settlement 
takes  place  between  late  winter  and  early  summer  (Hines  1979;  Pineda  and  Lopez  2002). 

The  study  region  encompassed  -400  kilometers  of  shoreline,  from  Dana  Point 
(33°27.78’N,  117°42.92’W)  in  Southern  California,  USA,  to  Punta  Baja  (29°57.25’N, 
115°48.64,W)  in  Baja  California,  Mexico  (Fig.  4.1).  Seven  sites  were  visited 
approximately  monthly  over  16-month  period  (Table  4.1).  On  each  visit,  ca.  100  adult 
Balanus  glandula  were  collected  haphazardly  from  an  area  that  spanned  -  100  m  of 
coastline  and  the  entire  vertical  range  of  distribution.  Individuals  were  preserved  in  95% 
ethanol  for  further  analysis  in  the  laboratory. 

Pictures  of  four  5x5  cm  quadrats  were  taken  on  each  visit  using  a  4  megapixel 
digital  camera  (Nikon  Coolpix  4500,  Nikon  Inc.).  Low  magnification  pictures  of  whole 
quadrats  were  taken,  as  well  as  high-magnification  shots  of  smaller  portions  that  resol  ved 
the  smallest  settlers  (-300  pm).  Thermistors  (Tidbit  loggers,  Onset  Computer)  that 
recorded  temperature  at  hourly  intervals  were  deployed  in  the  lower  intertidal  at  each  site. 
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4.3.2.  Measurements  and  data  analysis 


Up  to  100  individuals  from  each  site  and  sampling  month  were  measured  and 
inspected  for  reproductive  status.  They  were  classified  into  three  different  stages: 
immature  (no  apparent  egg  masses  or  gonad  development),  mature  with  no  egg  masses 
(yolk  accumulation,  no  distinct  egg  mass),  and  mature  with  egg  masses  (egg  masses 
distinct,  with  individual  eggs  clearly  visible).  Egg  masses  that  appeared  intact  were 
removed  from  gravid  individuals  and  preserved  in  95%  ethanol  for  subsequent  fertility 
analysis  (see  below). 

The  following  measurements  were  taken  to  the  nearest  0.01  mm  from  each 
individual:  basal  length  (longest  axis  of  base),  basal  width  (perpendicular  to  longest  axis), 
opercular  length,  opercular  width  (perpendicular  to  opercular  length),  and  height.  We 
investigated  the  existence  of  spatial  trends  in  morphometry  with  a  Principal  Component 
Analysis  (PCA)  of  measurements  collected  at  all  sites  and  sampling  times  (N  =  6,957 
individuals).  In  a  PCA,  new  variables  (i.e.  components)  are  computed  as  linear 
combinations  of  the  variables  originally  measured  (Morrison  1976).  The  sets  of 
coefficients  that  defined  those  linear  combinations  in  our  case  corresponded  to  the 
eigenvectors  of  the  correlation  matrix  for  all  five  sets  of  measurements.  We  used  the 
principal  component  that  explained  the  largest  percentage  of  total  variance  (i.e. 
eigenvector  corresponding  to  the  largest  eigenvalue)  to  summarize  size  information  along 
the  study  region,  and  to  conduct  between-site  comparisons  of  mean  size  for  each  season. 
Mean  sizes  within  each  season  were  compared  with  a  one-way  Analysis  of  Variance 
(ANOVA)  that  used  sites  as  treatments.  Homogeneity  of  variances  was  tested  prior  to  this 
analysis  using  the  Hartley  test  of  largest  to  smallest  variance  ratio  (Sokal  and  Rohlf  1995). 


4.3.3.  Estimation  of  growth  and  survival  rates 

Photographs  collected  over  16  months  were  used  to  estimate  a  growth  function  for 
B.  glandula  at  the  extremes  of  the  study  region,  Dana  Point  and  Punta  Baja.  Individuals 
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that  had  been  photographed  repeatedly  at  those  sites  for  1-14  months  were  tracked  and 
measured.  Lack  of  individuals  that  survived  for  sufficiently  long  periods  after  settlement 
precluded  the  estimation  of  growth  parameters  at  other  sites.  Assignments  of  age  to 
recently  settled  individuals  were  based  on  a  set  of  daily  measurements  of  settlers  in  La 
Jolla,  Southern  California.  Individuals  with  opercular  lengths  of  440  -  530  pm  were 
assigned  age  0  months  with  an  error  assumed  to  be  no  greater  than  a  few  days. 

Size-at-age  data  were  used  to  estimate  the  parameters  of  a  von  Bertalanffy  growth 
function: 


Lt=Lm( (4.1) 

where  t  is  time,  Loo  is  the  asymptotic  maximum  size,  k  is  a  parameter  proportional  to  the 
rate  at  which  L{  approaches  Loo,  and  to  corresponds  to  the  age  expected  for  an  individual 
of  size  zero.  Estimates  of  L^  k  and  to  were  obtained  for  each  site  with  a  non-linear 
ordinary  least  squares  regression.  The  two  sets  of  parameters  were  compared  using  a 
randomization  algorithm  devised  by  Weinberg  and  Helser  (1996),  which  is  briefly 
outlined  below. 

Under  the  null  hypothesis  of  no  difference  in  growth  parameters  between  sites, 
age-size  pairs  are  all  drawn  from  the  same  population.  If  that  is  the  case,  a  growth  curve 
fitted  to  the  combined  dataset  of  size  Nq  -  N i  +  yV2,  where  1  and  2  are  site  indices,  should 
not  differ  from  curves  fitted  separately  for  each  site  using  N\  and  N2  observations.  Thus, 
the  sum  of  squared  residuals  obtained  from  a  curved  fitted  to  the  combined  dataset  (SSRq) 
should  not  be  greater  than  the  sums  of  squared  residuals  obtained  from  fitting  two 
separate  curves  (SSR\+SSR2),  hence  the  test  statistic 

D0  =  SSRe  -  ( SSR ,  +  SSR2 )  (4.2) 
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should  be  small.  A  large  value  of  D0  indicates  that  one  curve  fitted  to  the  combined 
dataset  cannot  explain  part  of  the  variability  that  two  separate  site-specific  fits  do  explain, 
hence  there  is  a  small  chance  that  the  null  hypothesis  is  true. 

Prior  to  the  calculation  of  D0,  a  von  Bertalanffy  growth  curve  was  fitted  to  the 
combined  set  of  size-age  observations  gathered  from  Dana  Point  (N\=28)  and  Punta  Baja 
(7^2=57).  The  significance  of  Do  was  evaluated  by  comparing  it  with  a  probability 
distribution  obtained  from  the  repeated  computation  of 


(4.3) 


where  the  asterisks  indicate  sums  of  squared  residuals  from  growth  curves  fitted  to 
groups  of  Nj  and  N?  observations  that  are  randomly  allocated  to  each  group  from  the 
pooled  dataset.  This  random  allocation  of  data  pairs  and  computation  of  Dx  was  repeated 
Nt=  10,000  times,  and  a p-value  was  computed  as  P(DX  >  Do)  =  (x+l)/(AVH),  where  x  is 
the  number  of  iterations  in  which  Dx  >  Do  (Edgington  1986). 

Survival  probabilities  were  estimated  separately  for  adults  and  juveniles  from 
frequency  distributions  of  minimum  longevity,  which  was  measured  as  the  number  of 
months  an  individual  was  photographed  alive.  Lack  of  sufficient  data  precluded  the 
estimation  of  survival  rates  at  all  sites.  Only  Dana  Point  and  Punta  Baja,  together  with 
pooled  observations  from  La  Jolla  and  San  Antonio  (hereby  referred  to  as  LJ-SA), 
provided  sufficient  data  to  estimate  survival  rates  of  adults  and  juveniles.  The  assignment 
of  juvenile  or  adult  status  to  each  individual  in  the  photographs  was  done  by  comparing 
their  size  with  the  average  minimum  opercular  length  of  gravid  individuals  found  in  the 
field  samples  (2.03  mm,  SD  =  0.33  mm).  Because  photographs  were  taken  monthly, 
many  juveniles  were  only  photographed  once.  Minimum  ages  for  these  specimens  were 
estimated  from  their  opercular  length  and  the  above  mentioned  growth  curves. 

The  probability  of  surviving  until  time  /  was  computed  as 
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(4.4) 


/>(/</)  =  ! 


where  m  is  the  number  of  individuals  with  longevity  i  and  N  is  the  total  number  of 
individuals  tracked  in  the  photographs.  Temporal  changes  in  survival  probability  were 
modeled  for  t  -  1-12  months  using  a  negative  exponential  function,  which  was  fitted  to 
probabilities  computed  as  in  Eq.  4.4  by  least-squares  regression. 


4.3.4.  Fertility  analysis 

A  subset  of  gravid  individuals  collected  during  the  2002-2003  winter  was  used  to 
conduct  an  analysis  of  fertility  at  size,  i.e.  number  of  eggs  per  brood  as  a  function  of 
individual  size.  Individuals  spanning  the  range  of  sizes  found  at  each  site  were  chosen 
among  those  with  intact  egg  masses.  A  volumetric  method  was  used  to  estimate  the  total 
number  of  eggs  per  brood.  The  egg  masses  extracted  from  an  individual  barnacle  were 
suspended  in  water  and  separated  mechanically  using  a  combination  of  vortexing  and 
manual  brushing,  after  which  they  were  placed  in  a  suspension  of  80  ml  of  water.  Ten  1- 
ml  aliquots  were  examined  under  the  microscope,  and  the  numbers  of  eggs  were  recorded, 
averaged,  and  used  to  estimate  total  number  of  eggs.  Aliquots  were  taken  with  an 
automatic  pipette  in  sets  of  five  at  a  time,  so  as  to  minimize  the  effect  of  sampling 
without  replacement.  Extra  sets  of  aliquots  were  taken  whenever  the  coefficient  of 
variation  of  egg  counts  was  greater  than  20%.  Comparisons  of  volumetric  estimates  and 
the  total  number  of  eggs  counted  in  a  subset  of  individuals  indicated  that  the  method  was 
accurate  within  6-15%. 

A  Principal  Component  Analysis  (PCA)  was  performed  on  the  morphometric  data 
of  individuals  used  for  egg  counting,  in  order  to  find  the  best  predictor  variable  for  the 
analysis  of  fertility  at  size,  as  well  as  to  compare  the  size  structure  of  gravid  individuals 
between  sites.  The  first  principal  component  (PCI)  was  used  as  predictor  variable  in  a 
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reduced  major  axis  regression  of  log-transformed  fertility  versus  size  data.  This  type  of 
model  II  regression  is  recommended  for  cases  in  which  both  the  predictor  and  response 
variable  are  observed  with  error  (Rayner  1985).  Prior  to  the  analysis,  an  upper  cutoff  was 
chosen  for  the  predictor  variable,  so  that  comparisons  would  be  made  among  regression 
lines  fitted  within  similar  size  ranges.  A  threshold  of  20  mm  was  chosen  for  PCI,  which 
resulted  in  a  total  of  11  observations  removed  from  the  original  dataset  (2,  1,  1,  and  7 
from  DP,  LJ,  SA  and  ST,  respectively). 


4.3.5.  Stage-structured  matrix  population  model 

Observations  were  modeled  with  a  seasonal  two-stage  matrix  model  for  the 
population  of  Balanus  glandula  (Fig.  4.2).  The  model  relates  the  abundance  of  juvenile 
and  adult  individuals  with  transition  probabilities  at  two  different  times  of  the  year:  a 
settlement  season  centered  on  March,  and  a  reproductive  season  centered  on  December. 
The  choice  of  these  periods  was  based  on  the  available  literature  (e.g.  Hines  1979)  as 
well  as  on  the  author’s  field  observations.  The  separation  of  a  yearly  cycle  into  two 
different  transition  matrices  with  different  projection  intervals  is  warranted  by  the 
difference  in  time  scales  at  which  such  transitions  take  place  (Caswell  2001).  On  average, 
the  transition  between  settlement  and  reproduction  takes  place  over  a  period  of  nine 
months,  whereas  the  transition  between  larval  production  and  settlement  as  juveniles 
occurs  within  3  months  (Fig.  4.2),  For  the  first  transition,  the  number  of  adults  (^)  that 
will  be  alive  to  reproduce  in  December  is  a  function  of  the  survival  probabilities  for 
individuals  that  were  juveniles  (a)  and  adults  ( b )  in  March: 


n2  =  (a  b) 


(  »  \ 


\n2j 


(4.5) 


In  the  second  transition  (December  to  March),  the  numbers  of  juvenile  and  adult 
individuals  are  given  by 
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n 


(4.6) 


where  d  is  the  probability  of  adult  survival  from  December  to  March,  and  the 
reproductive  term  fs  (Fig.  4.2)  is  a  product  of  the  per  capita  fertility  of  adults  (/)  and  the 
survival  of  larvae  in  the  plankton  (5). 

Both  transitions  can  be  combined  to  obtain  one  matrix  equation  that  projects  the 
number  of  individuals  from  one  yearly  census  to  the  next.  However,  different  equations 
are  obtained  depending  on  the  time  of  year  for  which  projections  arc  made.  For  example, 
if  a  census  is  conducted  every  December,  there  will  be  no  juveniles  in  the  population,  and 
the  abundance  of  adults  will  be  given  by  the  juveniles  that  survived  after  settlement  plus 
adults  that  survived  from  the  previous  year.  The  equation  that  projects  the  population 
from  year  to  year  is  found  by  substituting  Eq.  (4.6)  into  Eq.  (4.5): 


i2(/  +1)  =  (afs  +  bd)n2{t) 


(4.7) 


n 


where  afs  +  bd  =  X  is  the  population  growth  rate  (i.e.  rate  at  which  the  population  will 
increase  from  year  to  year). 

For  a  March-March  projection  interval,  both  juveniles  and  adults  are  encountered 
on  each  census,  and  a  projection  matrix  is  found  by  substituting  Eq.  (4.5)  into  Eq.  (4.6): 


(4.8) 


where  afs  amounts  to  a  recruitment  probability  (i.e.  production  and  survival  of  larvae 
times  post-settlement  survival),  ad  is  the  probability  of  surviving  the  first  year  of  benthic 
life,  bfs  is  the  reproductive  output  of  surviving  adults,  and  bd  is  the  year-round 
probability  of  adult  survival  (Fig.  4.2). 
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If  the  projection  matrix  in  Eq.  (4.8)  is  separated  into  a  transition  matrix  (T)  and  a 
reproduction  matrix  (F), 


T  +  F  = 

'  0 

0  " 

+ 

afs 

bfs 

ad 

bd 

_  0 

0 

(4.9) 


formulas  to  compute  life  expectancies  and  net  reproductive  rates  can  be  obtained  from 
the  fundamental  matrix  (N) 


N  =  (l  -  T)~' 


(4.10) 


where  I  is  the  identity  matrix.  Life  expectancies  (i.e.  number  of  years  an  individual  is 
expected  to  live)  for  juveniles  and  adults  in  the  population  were  calculated  as  the  sum  of 
the  first  and  second  column  of  this  matrix,  respectively.  Finally,  the  largest  eigenvalue  of 
the  matrix  that  results  from  pre-multiplying  the  fundamental  matrix  (Eq.  4.10)  by  the 
reproductive  matrix  (F  in  Eq.  4.9)  corresponds  to  the  Net  Reproductive  Rate  ( Rq ) 


=■ 


afs 

1  —  bd 


(4.11) 


which  is  the  mean  number  of  individuals  a  juvenile  will  produce  by  the  end  of  its  life, 
hence  the  rate  by  which  the  population  will  increase  from  one  generation  to  the  next 
(Caswell  2001). 

I  studied  the  asymptotic  and  transient  behavior  of  the  model  when  parameterized 
for  three  different  localities  along  the  study  region:  Punta  Baja,  La  Jolla  -  San  Antonio, 
and  Dana  Point.  The  original  goal  was  to  estimate  model  parameters  for  at  least  three 
sites,  at  the  extremes  and  middle  of  the  study  region.  Insufficient  data  precluded  the 
satisfactory  parameterization  of  the  model  for  any  single  locality  other  than  Punta  Baja 
and  Dana  Point.  Observations  taken  at  La  Jolla  and  San  Antonio  were  pooled  in  order  to 
obtain  a  third  case  study.  For  each  site,  population  growth  rate  (k=afs+bd)  and  net 
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reproductive  rate  (Eq.  4.1 1)  were  computed,  as  well  as  their  sensitivity  to  changes  in  the 
model  parameters.  Sensitivities  were  calculated  as  partial  derivatives  of  the  respective 
rate  with  respect  to  each  parameter.  Life  expectancies  were  calculated  for  juveniles  and 
adults  at  each  site. 

To  analyze  the  effect  that  interactions  between  parameters  could  have  on 
population  persistence,  1  first  defined  persistence  as  a  population  state  where  the  growth 
rate  X.  >  1 .  A  condition  for  persistence  in  terms  of  the  model  parameters  can  be  found  by 
substituting  the  definition  of  A.  (Eq.  4.7)  into  this  inequality  and  solving  for  the 
reproductive  term  fs: 


Ji>b“  (4.12) 

a 

where  the  ratio  on  the  right  hand  side  is  proportional  to  mortality  in  the  benthic  part  of 
the  life  cycle.  Thus,  the  population  will  persist  if  the  reproductive  term  fs  -  which  would 
be  proportional  to  larval  supply  -  balances  the  combined  effect  of  juvenile  and  adult 
mortality.  Dividing  both  sides  of  Eq.  4.12  by  (1  -bd)/a,  the  left  hand  side  of  the  inequality 
is  equivalent  to  the  Net  Reproductive  Rate  (Eq.  4.1 1)  and  the  condition  for  persistence  is 
reduced  io  Ro>  1-  Thus,  Ro  could  also  be  interpreted  as  the  ratio  by  which  larval  supply 
exceeds  mortality  in  a  population,  or  as  how  far  a  population  is  from  crossing  the  line 
between  persistence  and  local  extinction.  Assuming  that  per  capita  fertility  (/)  is  constant 
at  any  given  site  -  owing  to  a  characteristic  size  distribution  of  reproductive  individuals 
and  a  space  invariant  fertility  function  -  I  calculated  the  minimum  per  capita  fertility 
required  to  attain  persistence  given  a  combination  of  survival  probabilities. 

To  explore  the  effects  of  environmental  stochasticity  on  the  persistence  of  each 
population,  I  simulated  random  interannual  fluctuations  in  the  parameters  to  which 
growth  and  net  reproductive  rates  were  most  sensitive.  Two  scenarios  were  used.  In  one 
of  them,  larval  survival  varied  from  year  to  year  while  all  other  parameters  were  kept 
constant.  In  the  second  scenario,  the  stochastic  parameter  was  juvenile  survival.  Values 
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for  larval  survival  (5)  and  juvenile  survival  (a)  were  drawn  at  random  from  uniform 
probability  distributions.  Parameters  for  these  distributions  were  chosen  such  that  the 
range  of  possible  values  would  encompass,  and  go  beyond,  those  estimated  from  our  field 
observations.  For  example,  ranges  used  for  the  lower  and  upper  limits  of  the  larval 
survival  distribution  were  smjn  =  .00005  -  .0025  and  ^max  =  .0025  -  .01,  which  assuming  a 
constant  mortality  rate  and  a  21 -day  larval  period  would  correspond  to  instantaneous 
mortality  rates  of  ca.  0.5  -  0.3  and  ca.  0.3  -  0.2  d’1,  respectively.  Recent  estimates  of 
instantaneous  mortality  rates  for  B.  glandula  larvae  ranged  from  0.2  to  0.4  d'1  (see 
Chapter  2) 

For  each  site  and  simulation,  the  transition  matrix  and  initial  conditions  «i=l 
juveniles  and  n2=\  adults  were  used  to  project  population  size  over  7=40,000  time  steps. 
On  each  time  step,  a  value  for  the  respective  stochastic  survival  parameter  was  taken  at 
random  from  its  probability  distribution.  The  log  of  the  mean  stochastic  growth  rate  ( Xs ) 
was  estimated  as 


r-i 


lo§4  =7Z,0g 
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(4.13) 


where  N{  is  the  total  number  of  individuals  in  the  population  at  time  t.  Population  size  N(t) 
is  lognormally  distributed  when  the  population  is  projected  over  a  sufficiently  long  period, 
so  that  log  Xs  is  normally  distributed  with  mean  t  log  Xs  and  variance  tcf  (Caswell  2001). 
Survival  values  in  the  simulations  were  independent  and  identically  distributed,  hence  it 
was  possible  to  estimate  c?  (i.e.  rate  of  increase  of  the  variance  of  log  TV)  as 

a2  =  2(log /i  -  log As)  (4.14) 


where  log  p  is  the  growth  rate  of  mean  population  size,  computed  as  the  log  of  the 
dominant  eigenvalue  of  the  average  projection  matrix.  This  estimate  was  used  to  compute 
the  probability  of  quasi-extinction,  defined  as  the  probability  that  population  size  declines 
below  a  fraction  6  of  the  initial  population  size,  and  calculated  as 
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if  log/ls  <  0 


(4.15) 
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if  \ogXs  >  0 


4.3.6.  Estimation  of  model  parameters 

Mean  per  capita  fertility  if)  was  estimated  from  the  mean  size  of  gravid 
individuals  at  each  site  and  the  corresponding  fertility  regression.  Confidence  limits  for / 
were  obtained  from  bootstrap  distributions,  which  were  computed  by  resampling  with 
replacement  from  the  size  and  fertility  data  available  for  each  site  (see  Manly  1997).  A 
reduced  major  axis  regression  of  fertility  (see  above)  was  performed  on  each  of  Ni(  = 
10,000  iterations  per  site,  using  bootstrap  datasets  sampled  from  the  original  sets  of 
fertility-at-size  observations.  The  same  number  of  size  datasets  was  generated  by 
bootstrapping  from  size  observations  at  each  site,  and  used  to  compute  the  mean  size  of 
gravid  individuals.  Size  and  fertility  observations  for  La  Jolla  and  San  Antonio  were 
pooled  pnor  to  this  procedure. 

The  mean  estimate  for  larval  survival  (5)  was  based  on  recent  estimates  of  the 
instantaneous  mortality  rate  of  larval  B.  glandula  in  nearshore  waters  off  La  Jolla, 
Southern  California  (see  Chapter  2).  I  assumed  that  the  number  of  surviving  larvae 
decayed  exponentially  at  a  constant  rate  of  0.2  d"1  over  a  21 -day  period.  The  choice  of  21 
days  for  larval  duration  was  based  on  published  observations  of  development  time 
(Brown  and  Roughgarden  1985)  at  temperatures  similar  to  those  recorded  at  the  intertidal 
sites  during  the  reproductive  season.  Mean  estimates  for  the  survival  probabilities  of 
juveniles  and  adults  ( a ,  b  and  d)  were  obtained  from  negative  exponential  curves  fitted  to 
survival  data  from  each  site  (see  above).  All  of  these  mean  estimates  (5,  a ,  b,  cl)  had 
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associated  sampling  errors  that  I  attempted  to  include  in  the  calculations  of  population 
growth  rate,  net  reproductive  rate,  and  life  expectancies  with  a  Monte  Carlo  procedure. 
Rather  than  using  only  one  parameter  value  on  each  rate  calculation,  vectors  of  values 
were  generated  by  random  sampling  from  their  respective  probability  distributions. 
Values  of  larval  survival  (s)  were  obtained  from  a  uniform  distribution  with  limits 
defined  by  the  fraction  of  larvae  that  would  survive  after  14  -  28  days  with  a  constant 
mortality  rate  of  0.2  d  1  (smin  =  0.0037,  smQX  =  0.061).  Juvenile  survival  (a)  was  sampled 
from  a  uniform  distribution  that  ranged  between  the  minimum  and  maximum  estimates 
obtained  from  field  measurements.  Adult  survival  after  3  and  9  months  ( d  and  b)  was 
calculated  using  the  exponential  curves  fitted  to  survival  data  from  each  site:  y  =  e''x, 
where  y  is  survival,  x  is  time  (in  months)  and  the  parameter  r  is  a  normally  distributed 
random  variable  with  mean  and  standard  deviation  given  by  its  least-squares  estimator 
and  standard  error.  On  each  iteration  of  the  Monte  Carlo  procedure,  a  different  vector  of 
parameter  values  was  used  to  compute  growth  rate,  net  reproductive  rate,  and  juvenile 
and  adult  life  expectancies.  I  used  10,000  repetitions  to  obtain  probability  distributions 
from  which  95%  confidence  limits  were  obtained  using  the  bias-corrected  percentile 
method  (Manly  1  997). 
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4.4.  Results 


4.4. 1.  Size  structure  and  timing  of  reproductive  activity 

A  Principal  Component  Analysis  of  morphometric  data  for  all  sites  and  measured 
individuals  (Table  4.2)  produced  a  first  principal  component  that  explained  75.4%  of  the 
total  variability  in  size  structure.  A  new  variable  (PCI)  was  produced  with  this  set  of 
weights,  and  is  used  for  the  remainder  of  this  paper  to  describe  the  size  structure  of 
Balanus  glandula.  The  weights  used  to  compute  PCI  from  five  originally  measured 
variables  were  very  similar  in  magnitude  relative  to  PC2  and  PC3  (Table  4.2),  hence  the 
new  variable  approximately  corresponds  to  the  summation  of  all  five  measurements,  and 
could  be  visualized  as  the  outer  line  or  external  perimeter  of  each  individual. 

The  frequency  distributions  of  PC  I  (Fig.  4.3)  were  similar  in  shape  for  some  pairs 
of  neighboring  sites  (e.g.  San  Quintin  and  Punta  Baja  in  winter  2003-2004  and  summer 
2003),  with  a  weak  trend  of  decreasing  mean  size  from  Dana  Point  to  Punta  Baja. 
Significant  between-site  differences  were  found  when  means  were  compared  for  each 
season  using  a  one-factor  ANOVA  (Z^oos =  62.03  for  winter  2002-03,  ^6,2004  -  91 .68  for 
summer  2003,  /r6,i927==  40.85  for  winter  2003-04,  F6,n6i  =  41.60  for  gravid  individuals; 
all  P<0.001).  Bonferroni-adjusted  pairwise  comparisons  showed  consistent  spatial 
patterns  (Table  4.3).  For  example,  differences  in  mean  sizes  recorded  at  the  extremes  of 
the  study  region  were  always  significant  (e.g.  DP-LJ  versus  SQ-PB).  San  Quintin  and 
Punta  Baja  were  significantly  different  from  the  rest,  but  not  from  each  other,  on  two  out 
of  three  seasons  and  when  gravid  individuals  were  compared.  Finally,  differences  in 
mean  size  were  always  significant  between  La  Jolla  and  sites  located  south  of  it  (Table 
4.3).  In  addition  to  spatial  differences  in  their  means,  size  distributions  recorded  at  the 
extremes  of  the  study  region  differed  significantly  in  their  shapes.  A  comparison  of 
standardized  residuals  from  a  two-way  contingency  table  analysis  of  PCI  distributions 
showed  that  large  individuals  were  more  abundant,  and  smaller  ones  more  scarce,  than 
expected  at  Dana  Point  and  La  Jolla  under  the  hypothesis  of  independence  (Fig.  4.4).  The 
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opposite  was  true  for  San  Quintin  and  Punta  Baja:  small  and  large  individuals  were  more 
and  less  numerous,  respectively,  than  expected  under  the  hypothesis  of  independence 
between  site  and  size  frequency  (Fig.  4.4). 

Temporal  changes  in  the  proportion  of  brooding  individuals  occurred 
synchronously  across  the  study  region  at  the  monthly  sampling  interval  (Fig.  4.5).  A 
small  and  narrow  peak  in  brooding  activity,  centered  on  March  2003,  was  observed 
during  the  first  winter.  Brooding  increased  earlier  in  the  reproductive  season  and 
remained  high  for  a  longer  time  during  the  second  winter  (Fig.  4.5),  with  maximum 
percentages  observed  between  January  and  March  2004.  It  is  possible  that  the  small  peak 
recorded  during  the  first  winter  of  our  study  corresponded  to  a  late-season  pulse  of 
reproduction.  High  percentages  of  brooding  individuals  were  found  at  San  Quintin  and 
Santo  Tomas  in  September  2002  (author’s  personal  observations),  suggesting  that 
reproduction  started  earlier  on  that  season.  The  synchronous  increase  in  brooding  activity 
observed  between  November  2003  and  January  2004  coincided  with  a  downward  trend  in 
intertidal  temperatures  (Fig.  4.5).  Differences  in  temperature  of  up  to  7  °C  between  Dana 
Point  and  Punta  Baja  were  observed  from  spring  to  late  summer. 


4.4.2 .  Growth  and  survival  rates 

Growth  parameters  estimated  from  size-at-age  observations  at  each  site  were 
significantly  different  (D0  =  3.872,  P  =  0.0001)  and  indicated  that  individuals  at  Punta 
Baja  grew  more  slowly  than  at  Dana  Point  (ApB=0.11  <  &dp=0.20,  Fig.  4.6).  Estimated 
asymptotic  maximum  sizes  ( L «>)  suggest  that  barnacles  at  Punta  Baja  could  grow  to  larger 
sizes  than  those  at  Dana  Point  (Toopb=4.20  >  Zoodp^.OS,  Fig.  4.6).  However,  the 
comparatively  low  estimate  of  adult  survival  obtained  for  Punta  Baja  (Fig.  4.7)  suggests 
that  these  individuals  are  unlikely  to  survive  the  time  necessary  to  attain  larger  sizes.  This 
is  consistent  with  observed  differences  in  size  distributions  at  the  extremes  of  the  region 
(Fig.  4.4).  A  rapid  decline  in  the  survival  of  juveniles  was  observed  at  all  sites  during  the 


99 


first  four  months  after  settlement  (Fig.  4.7).  This  decrease  was  less  pronounced  at  Punta 
Baja,  where  the  probability  of  a  juvenile  to  survive  for  nine  months  after  settlement  was 
greater  than  at  LJ-SA  and  Dana  Point  (Fig.  4.7). 


4.4.3.  Fertility  at  size  and  mean  per  capita  fertility 

A  PCA  of  morphometric  data  for  gravid  barnacles  utilized  in  the  fertility  analysis 
produced  a  first  principal  component  that  explained  81.7%  of  the  total  variability  in  size 
structure  (Table  2).  As  with  the  PCA  performed  on  the  full  dataset,  these  weights  were 
very  similar  in  magnitude,  so  the  interpretation  of  the  variable  PCI  as  body  perimeter 
remained  valid.  The  second  principal  component  explained  an  additional  7.8%  of  the 
total  variability  (Table  2).  Scores  computed  with  these  two  components  for  all  individuals 
used  in  the  fertility  analysis  showed  no  apparent  separation  between  sites  when  plotted  on 
a  PC1-PC2  plane  (results  not  shown),  suggesting  that  differences  in  morphometry  could 
be  ruled  out  as  the  cause  of  any  between-site  differences  in  size  specific  fertility. 

With  the  exception  of  Santo  Tomas,  fertility  regressions  yielded  similar 
relationships  between  number  of  eggs  and  individual  size  throughout  the  study  region 
(Fig.  4.8).  The  slope  of  the  fertility  function  at  Santo  Tomas  was  larger  than  those 
estimated  for  neighboring  sites.  However,  there  was  partial  overlap  among  bootstrap 
confidence  intervals  when  a  correction  was  made  to  allow  for  pairwise  comparisons  with 
a  group  significance  of  a  =  0.05  (Table  4.4).  Similar  estimates  of  mean  per  capita  fertility 
were  obtained  for  Dana  Point  and  LJ-SA  (ca.  6,000  eggs),  whereas  an  average  individual 
at  Punta  Baja  was  expected  to  produce  less  than  half  as  many  eggs  (Table  4.5).  Bootstrap 
distributions  of  per  capita  fertility  were  symmetrical  around  mean  estimates  for  each  site 
(Fig.  4.9)  and  very  similar  for  Dana  Point  and  LJ-SA,  with  ca.  70%  overlap  in  their  95% 
confidence  intervals  (Table  4.5).  There  was  no  overlap  between  these  two  distributions 
and  that  obtained  for  Punta  Baja  (Fig.  4.9,  Table  4.5). 
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4.4.4.  Matrix  population  model 


Results  from  the  analyses  detailed  in  previous  sections  were  used  to  parameterize 
the  projection  matrix  in  Eq.  (4.8)  for  Dana  Point,  LJ-SA,  and  Punta  Baja.  Median 
estimates  of  population  growth  rate  ()i)  were  greater  than  one  at  all  three  sites  (Table  4.5), 
with  values  obtained  for  the  northern  sites  being  three  times  higher  than  at  Punta  Baja 
(Table  4.5).  Probability  distributions  of  X  were  very  similar  in  shape  and  width  at  Dana 
Point  and  LJ-SA,  with  wide  confidence  intervals  that  almost  completely  overlapped 
(Table  4.5).  Greater  differences  were  found  in  the  probability  distributions  of  Net 
Reproductive  Rate  (Ro).  Median  estimates  were  all  greater  than  one,  but  increased  by  a 
factor  of  two  from  Dana  Point  to  LJ-SA,  and  decreased  by  a  factor  of  thirteen  between 
LJ-SA  and  Punta  Baja  (Table  4.5).  In  terms  of  the  condition  for  persistence  previously 
defined  (Eq.  4.12),  Punta  Baja  was  thirteen  times  closer  than  LJ-SA  to  crossing  the  line 
between  persistence  and  local  extinction.  A  very  narrow  and  positively  skewed 
distribution  at  Punta  Baja  (Fig.  4.10b)  yielded  a  2%  probability  of  finding  values  of  Rq  < 

1 .  The  difference  in  R0  between  Dana  Point  and  LJ-SA  was  due  to  differences  in 
estimated  adult  survival  (e.g.  Table  4.5),  which  also  generated  substantial  differences  in 
the  probability  distributions  of  juvenile  and  adult  life  expectancies  at  these  sites  (Fig. 
4.10c,d).  Juveniles  and  adults  at  LJ-SA  were  expected  to  live  0.23  and  2.53  years  longer 
than  at  Dana  Point.  Distributions  of  juvenile  and  adult  expectancies  at  Punta  Baja  were 
extremely  narrow  and  close  to  one  (Fig.  4.10c,d  and  Table  4.5),  indicating  that  barnacles 
at  this  site  are  not  likely  to  be  in  the  population  for  much  longer  than  one  yearly  census. 

Sensitivity  calculations  indicated  that  growth  rates  and  net  reproductive  rates  were 
disproportionately  more  sensitive  to  changes  in  the  survival  of  early  stages  of  the  life 
cycle.  Larval  survival  (s)  had  the  largest  effect  at  all  three  sites,  although  the  importance 
of  juvenile  survival  (a)  appeared  to  increase  at  LJ-SA  (Table  4.6).  Based  on  these 
findings  and  the  assumption  that  per  capita  fertility  (/)  is  constant  for  a  given  site,  I  used 
the  inequality  in  Eq.  4.12  to  compute  minimum  fertilities  required  to  attain  persistence 
under  different  combinations  of  larval  survival  and  ‘benthic  mortality  ratios’,  i.e.  (1  -bd)!a. 
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The  results  showed  that  northern  sites,  given  their  per  capita  fertility,  could  persist  with 
larval  survival  an  order  of  magnitude  lower  than  values  that  would  drive  Punta  Baja  to 
local  extinction  (Fig.  4.1 1).  The  horizontal  arrows  on  Figure  4.1 1  point  to  survival  values 
of  0.0029,  0.0006  and  0.0005  for  Punta  Baja,  LJ-SA  and  Dana  Point,  respectively. 
Assuming  an  exponential  decay  in  the  number  of  surviving  larvae  and  an  average 
development  time  of  21  days,  these  values  would  correspond  to  instantaneous  mortality 
rates  of  0.28  d’1,  0.35  d"1,  and  0.36  d'1.  Mortality  rates  that  are  comparable  to  these  values 
have  been  recently  estimated  for  the  larvae  of  Balanus  glandula  in  nearshore  waters  off 
La  Jolla  (see  Chapter  2). 

Sensitivity  calculations  suggested  that  the  persistence  of  a  local  population  could  be 
greatly  affected  by  temporal  fluctuations  in  the  survival  probabilities  of  larvae  (5)  and 
juveniles  ( a ).  The  sensitivity  of  Rq  to  these  parameters  was  a  function  of  per  capita 
fertility  and  adult  survival  (Table  4.6),  which  showed  the  largest  differences  between 
Punta  Baja  and  the  two  northern  sites  (Table  4.5).  The  effects  that  interannual  variability 
in  these  parameters  would  have  on  population  persistence  are  summarized  in  Figures  4.12 
and  4.13.  Each  panel  in  these  figures  shows  the  results  of  900  simulations,  corresponding 
to  900  different  probability  distributions  from  which  survival  values  were  drawn  at 
random.  Contours  on  Figure  4.12  correspond  to  the  stochastic  growth  rate  (log  As),  which 
was  always  positive  for  Dana  Point  and  LJ-SA,  and  negative  on  ca.  30%  of  the 
simulations  for  Punta  Baja  (Fig.  4.12).  A  similar  pattern  was  observed  for  quasi- 
extinction  probabilities  (Fig.  4.13).  Dana  Point  and  LJ-SA  had  a  zero  probability  of 
quasi-extinction  throughout  the  parameter  space  used  in  the  simulations.  Quasi-extinction 
had  a  probability  of  one  in  30%  of  the  simulations  conducted  for  Punta  Baja  (Fig.  4.13). 
Consistent  with  Figure  4.11,  persistence  was  not  possible  at  Punta  Baja  when  larval 
survival  ranged  between  values  smaller  than  0.0029  (Fig.  4.13).  Persistence  at  Punta  Baja 
was  possible  only  when  years  with  low  survival  (of  either  larvae  or  juveniles)  were  as 
likely  to  occur  as  years  with  high  survival  (Fig.  4.13).  In  other  words,  this  population 
could  withstand  very  low  recruitment  years  only  if  years  with  very  high  recruitment  were 
equally  likely. 
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4.5.  Discussion 


I  have  found  spatial  differences  in  the  size  distribution,  growth,  and  survival  rates 
of  B alarms  glandula  within  a  400-km  range  in  Southern  and  Baja  California.  The  results 
suggest  that  a  range  endpoint  for  this  species  could  result  from  spatial  changes  in 
demography  in  a  temporally  variable  environment.  Previous  studies  based  on  spatially 
explicit  models  have  suggested  that  range  endpoints  in  coastal  populations  could  result 
from  the  interaction  between  larval  advection  and  substrate  availability  (Roughgarden  ct 
al.  1988),  or  from  advective  transport  only  (Gaylord  and  Gaines  2000).  On  the  other  hand, 
a  recent  field-based  experimental  study  shows  that  interactions  between  physiological 
stress  and  inter-specific  competition  could  limit  the  distribution  of  intertidal  species 
(Wethey  2002).  Here  1  used  field  observations  and  a  stage-structured  matrix  model  to 
show  that  a  spatial  gradient  in  vital  rates  and  environmental  stochasticity  are  sufficient 
conditions  to  generate  a  range  endpoint  for  Balarms  glandula  near  Punta  Baja. 


4.5.1.  Population  structure  and  vital  rates 

There  were  consistent  spatial  patterns  in  the  size  structure  of  Balanus  glandula. 
Both  the  mean  and  shape  of  size  distributions  changed  significantly  from  Dana  Point  to 
Punta  Baja  in  all  three  seasons.  Mean  sizes  were  significantly  smaller  at  Punta  Baja 
because  small  and  large  individuals  were,  respectively,  more  and  less  abundant  than 
expected  under  the  hypothesis  of  spatial  homogeneity  in  size  distributions.  Such 
differences  may  be  caused  by  differences  in  the  growth  rate  and  survival  of  adults  at 
Punta  Baja  and  Dana  Point. 

Spatial  changes  in  the  growth  rate  of  adult  Balanus  glandula  have  been  ascribed 
to  differences  in  temperature  (Phillips  2005),  food  availability  (Duggins  et  al.  1989),  or 
both  (Sanford  and  Menge  2001).  In  general,  warmer  temperatures  and  higher  food 
availability  should  increase  growth  rates,  although  temperature  appears  to  be  a  more 
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important  factor.  Differences  in  B.  glandula  growth  rates  recently  found  around  Point 
Conception  (Phillips  2005)  occurred  in  the  absence  of  patterns  in  chlorophyll-# 
concentration  (used  as  a  proxy  for  food  availability),  and  attributed  to  spatial  differences 
in  temperature.  These  results  are  consistent  with  the  differences  in  growth  found  between 
Punta  Baja  and  Dana  Point.  Growth  at  Punta  Baja  was  significantly  slower  than  at  Dana 
Point,  where  intertidal  water  temperatures  were  up  to  7°C  warmer  during  the  summer  and 
early  autumn.  Although  in  situ  chlorophyll  concentrations  were  not  measured,  shipboard 
measurements  conducted  in  the  spring  of  2002,  2003  and  2004  off  Punta  Baja  and  La 
Jolla  (Pineda  et  ah,  unpublished  results)  showed  higher  chlorophyll-#  concentrations  and 
colder  temperatures  in  nearshore  waters  off  Punta  Baja.  Satellite-derived  chlorophyll 
maps  for  the  region  (Espinosa-Carreon  et  al.  2004)  suggest  that  such  spatial  differences  in 
pigment  concentration  are  consistent  over  inter-annual  scales.  Thus,  differences  in  growth 
rate  between  the  extremes  of  the  sampling  region  could  be  caused  by  spatial  differences 
in  temperature  rather  than  food  availability.  Differences  in  adult  survival,  on  the  other 
hand,  could  be  linked  to  spatial  changes  in  a  number  of  physical  and  biological  factors, 
including  predation  pressure,  substrate  availability  and  substrate  quality.  Given  the  low 
sampling  frequency  used  in  this  study,  at  this  point  I  could  only  speculate  as  to  the  role  of 
these  factors  in  determining  survival.  Experimental  work  at  local  scales  and  higher 
sampling  frequencies  is  needed  to  gather  a  more  complete  picture  of  the  interaction 
between  physical  and  biological  factors  that  determine  the  observed  differences  in  vital 
rates. 


Although  size  distributions  were  significantly  different  between  the  extremes  of 
the  study  region,  no  spatial  differences  in  morphometry  (i.e.  relative  magnitudes  of 
different  body  dimensions)  were  detected.  This  finding  supported  the  use  of  one  global 
variable  (PCI)  to  describe  and  compare  the  size  structure  of  B.  glandula  along  the  study 
region.  Furthermore,  it  allowed  the  direct  comparison  of  fertility  regressions  from 
different  sites.  With  the  exception  of  Santo  Tomas,  between-site  differences  in  the  slope 
of  fertil  ity  regressions  were  small,  and  suggested  that  changes  in  reproductive  output  over 
regional  scales  are  determined  primarily  by  changes  in  the  size  structure  and  density  of 
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reproductive  adults.  Results  from  other  comparisons  of  crustacean  fertility  conducted  at 
similar  latitudes  and  spatial  scales  are  consistent  with  these  results.  For  instance,  Dugan 
et  al.  (1991)  found  no  significant  differences  in  the  slopes  of  their  fertility  regressions  for 
sand  crabs  Emerita  analoga  from  six  sites  along  400  kilometers  of  the  Southern 
California  coast. 


4.5.2 .  Stage  structured  matrix  model 

The  analysis  of  a  stage-structured  matrix  model  indicated  that  an  imbalance 
between  the  reproductive  term  fs  and  the  mortality  of  adults  and  juveniles  (Eq.  4.12)  is  a 
necessary  condition  for  population  persistence.  Field  observations  verified  that  such 
condition  holds  at  all  three  sites  for  which  the  model  was  parameterized.  However,  the 
factor  by  which  population  renewal  was  larger  than  the  benthic  mortality  ratio  - 
equivalent  to  the  net  reproductive  rate  (Table  4.5)  -  decreased  more  than  ten  times 
between  La  Jolla  -  San  Antonio  and  Punta  Baja,  a  site  that  is  near  the  southern  limit  of 
distribution  of  Balcinus  glandula.  This  latitudinal  change  in  net  reproductive  rate  suggests 
that,  as  barnacle  populations  get  closer  to  Punta  Baja,  they  get  closer  to  crossing  a 
threshold  between  persistence  and  local  extinction. 

Sensitivity  calculations  showed  that  larval  survival  and  juvenile  survival  had  the 
largest  proportional  effects  on  growth  rate  (X)  and  net  reproductive  rate  ( Ro ),  and  by 
extension  on  population  persistence.  The  simulations  confirmed  these  results  and 
suggested  that  inter-annual  fluctuations  in  the  survival  of  larvae  and  juveniles  could 
generate  between-site  differences  in  susceptibility  to  quasi-extinction,  defined  here  as  the 
event  of  population  size  falling  below  a  fraction  #=0.1  of  the  initial  population  size.  Dana 
Point  and  LJ-SA,  the  two  northern  sites,  appeared  robust  to  large  changes  in  the  survival 
of  larvae  and  juveniles.  These  sites  yielded  positive  estimates  of  the  log  mean  stochastic 
growth  rate,  and  zero  probability  of  quasi-extinction,  for  the  entire  range  of  parameter 
values  used  in  the  simulations.  Results  were  dramatically  different  for  Punta  Baja,  which 
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was  susceptible  to  quasi-extinction  under  ca.  30%  of  the  simulated  regimes  of 
stochasticity  in  larval  survival  and  juvenile  survival. 

The  stage-structured  matrix  model  was  formulated  in  a  non-spatial  context,  where 
larval  numbers  decrease  due  to  natural  mortality  (e.g.  predation  and  starvation)  as  well  as 
by  advection  and  diffusion.  Thus,  stochasticity  in  the  larval  survival  parameter  (s)  must 
reflect  the  combined  variability  of  biological  and  physical  processes  that  affect  mortality 
and  advection-diffusion  of  larvae  in  the  nearshore.  It  follows  that  spatial  differences  in 
physical  forcing  should  generate  spatial  differences  in  larval  supply.  In  the  context  of  the 
model,  this  amounts  to  saying  that  different  local  populations  may  experience  degrees  of 
stochasticity  found  at  different  regions  of  the  parameter  space  shown  on  Figures  4.12  and 
4.13  (left  panels),  which  could  increase  between-site  differences  in  susceptibility  to  local 
extinction. 

Physical  observations  conducted  in  recent  years  show  more  intense  advectivc 
forcing  at  Punta  Baja  than  at  La  Jolla.  Coastal  upwelling  is  more  prevalent  at  Punta  Baja, 
and  can  be  detected  as  net  offshore  flows  down  to  depths  of  ca.  40  m  (Lopez  and  Pineda 
2004).  Such  conditions  might  increase  the  offshore  transport  of  larvae  (e.g.  Roughgardcn 
et  al.  1988;  but  see  Garland  and  Zimmer  2002)  and  the  interannual  variability  in  larval 
supply  at  Punta  Baja.  Low  per  capita  fertilities  at  San  Quintin  and  Punta  Baja  could 
further  contribute  to  a  reduction  of  larval  abundance  in  the  nearshore. 


4.5.3.  Caveats  and  limitations  of  this  study 

The  spatial  coverage  that  allowed  to  show  geographic  patterns  in  the  demography 
of  Balanus  glandula  was  attained  at  the  expense  of  sampling  frequency.  Therefore,  it  was 
not  possible  to  simultaneously  investigate  processes  that  could  generate  some  of  these 
patterns.  Settlement  rate,  i.e.  number  of  larvae  settling  per  unit  time  per  unit  area  of 
available  substrate,  has  been  traditionally  used  as  a  proxy  for  larval  supply  (e.g.  Gaines  et 
al.  1985;  Miron  et  al.  1995),  although  small  scale  hydrodynamics,  behavioral  factors  and 
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space  availability  may  decouple  these  two  processes  (Pineda  2000).  Information  on 
settlement  rates  at  the  sampling  sites  could  have  provided  valuable  insight  on  inter-site 
differences  in  larval  supply  and  susceptibility  to  extinction.  Sampling  intervals,  however, 
would  have  had  to  be  much  shorter  (e.g.  1  day),  as  inter-site  differences  due  to  larval 
supply  can  be  obscured  by  differences  due  to  post-settlement  mortality. 

Another  limitation  imposed  by  a  low  sampling  frequency  was  the  lack  of  temporal 
resolution  in  the  observation  of  reproductive  timing.  Although  spatial  differences  in  the 
timing  or  duration  of  the  brooding  season  were  not  apparent  in  the  data,  the  frequency  of 
field  observations  was  probably  insufficient  to  resolve  differences  of  days  to  weeks  in  the 
timing  of  reproduction  across  the  study  region.  Laboratory  and  field  observations  indicate 
that  warm  water  may  shorten  the  brooding  season  of  Balanus  glandula  and  that,  to  some 
extent,  cold  water  would  have  the  opposite  effect  (Hines  1976,  1979).  Temperature 
differences  between  Dana  Point  and  Punta  Baja  fluctuated  around  7°C  in  spring-summer 
and  approached  zero  in  late  fall  to  early  winter.  The  effect  of  such  spatial  differences  in 
temperature  regime  on  the  timing  and  duration  of  a  brooding  season,  and  potentially  on 
the  number  of  broods  produced  per  individual,  remains  to  be  investigated.  Spatial 
differences  in  the  number  of  broods  produced  per  individual  in  one  season  could  change 
the  latitudinal  gradient  observed  in  mean  per  capita  fertility  and  perhaps  affect  the 
model's  predictions. 

Finally,  the  simplicity  of  this  matrix  model  made  it  possible  to  obtain  simple 
linear  expressions  for  growth  rate  ( X )  and  net  reproductive  rate  ( Ro ),  which  made  it 
possible  to  calculate  sensitivities  of  these  rates  to  individual  parameters.  However, 
leaving  density-dependent  effects  out  of  the  model  might  have  been  an  over 
simplification  that  could  potentially  affect  this  work’s  conclusions.  It  has  been  shown  that 
sharp  boundaries  in  spatial  distribution  could  be  generated  by  smooth  gradients  in  intra- 
specific  interactions  (Wilson  et  al.  1996),  many  of  which  could  be  attributed  to  density 
dependence.  For  example,  survival  rates  of  juveniles  and  adults  are  likely  to  depend  on 
local  density  (Svensson  et  al.  2004).  Per  capita  fertility  could  decrease  at  low  densities 
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due  to  an  Alice  effect  (Wethey  1984),  as  adult  barnacles  must  be  within  a  minimum 
distance  of  their  neighbors  in  order  for  fertilization  to  take  place  (Barnes  and  Crisp  1956). 
Settlement,  whether  gregarious  (Knight-Jones  and  Stevenson  1950;  Kent  et  al.  2003)  or 
proportional  to  free  space  (Gaines  and  Roughgarden  1985;  Hyder  et  al.  2001,  but  see 
Pineda  and  Caswell  1997),  could  also  be  a  function  of  local  adult  density. 

Thus,  the  possibility  that  a  range  endpoint  for  Balanus  glandula  at  Punta  Baja  is 
the  manifestation  of  one  or  more  density  dependent  effects  on  local  demography  cannot 
be  ruled  out.  I  can  only  argue  that,  in  the  framework  of  these  field  observations  and 
model  results,  spatial  differences  in  adult  survival  and  mean  per  capita  fertility  alone  can 
largely  affect  the  susceptibility  of  a  population  to  local  extinction,  and  possibly  the 
geographic  limits  of  a  species’  distribution. 
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Table  4.1 


Sampling  effort  for  each  site  and  month,  showing  the  numbers  of 
individuals  sampled,  measured  and  inspected  for  reproductive  condition. 
Sampling  sites  from  north  to  south  (see  Fig.  1),  are  Dana  Point  (DP),  La 
Jolla  (LJ),  San  Antonio  (SA),  Santo  Tomas  (ST),  Erendira  (ER),  San 
Quintin  (SQ),  and  Punta  Baja  (PB). 


Sampling 

Date 

S 

1  T 

E 

DP 

LJ 

SA 

ST 

ER 

SQ 

PB 

Winter  2002  -  2003 

November  2002 

117 

102 

96 

93 

94 

100 

92 

December  2002 

100 

97 

100 

- 

94 

100 

68 

January  2003 

99 

98 

100 

96 

93 

100 

- 

February  2003 

85 

54 

97 

53 

95 

101 

40 

March  2003 

100 

98 

91 

105 

75 

92 

87 

Summer  2003 

April  2003 

98 

79 

75 

43 

57 

37 

40 

May  2003 

92 

65 

99 

- 

- 

89 

100 

July  2003 

62 

86 

60 

76 

72 

47 

64 

September  2003 

100 

96 

99 

99 

94 

94 

88 

Winter  2003  -  2004 

November  03 

100 

69 

96 

73 

96 

100 

60 

December  03 

99 

95 

95 

99 

96 

100 

44 

January  04 

69 

56 

- 

69 

- 

- 

66 

March  04 

78 

64 

- 

79 

59 

64 

108 

Total 

1199 

1059 

1008 

885 

925 

1024 

857 
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Table  4.2 


Principal  Component  Analysis  performed  on  morphometric  data  collected  for 
Balanus  glandula  at  all  localities.  Columns  show  the  weights  assigned  to 
originally  measured  variables.  Weights  were  computed  as  eigenvectors 
corresponding  to  the  three  largest  eigenvalues  of  the  correlation  matrix  for  each 
dataset.  Only  the  components  that  explained  more  than  5%  of  the  total  variance 
are  reported. 


Variable 

Global  data 
(N  ~  6,957) 

(N- 

Fertility  subset 
=  202  gravid  indiv.) 

PCI 

PC2 

PC3 

PCI 

PC2 

PC3 

Basal  length 

0.459 

0.020 

-0.491 

0.456 

0.279 

-0.363 

Basal  width 

0.452 

0.011 

-0.588 

0.448 

0.247 

-0.576 

Opercular  length 

0.458 

-0.426 

0.377 

0.460 

-0.448 

0.114 

Opercular  width 

0.466 

-0.332 

0.376 

0.447 

-0.596 

0.159 

Height 

0.398 

0.841 

0.359 

0.425 

0.552 

0.706 

%  Total  var. 

75.408 

9.816 

8.357 

81.732 

7.786 

5.718 
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Table  4.3 


Pairwise  comparisons  of  mean  size  (PCI)  between  sites.  Black  bars  connect  sites 
with  no  significant  differences  in  their  mean  size  (i.e.  Bonferroni-adjusted  P  > 
0.05).  Means  are  sorted  in  decreasing  order  from  left  to  right  on  each  set  of 
comparisons.  Site  labels  correspond  to  Dana  Point  (DP),  La  Jolla  (LJ),  San 
Antonio  (SA),  Santo  Tomas  (ST),  Erendira  (ER),  San  Quintin  (SQ),  and  Punta 
Baja  (PB). 


Season  / 
Category 


Winter 

2002-2003 


SITE 

Decreasing  means - > 


LJ  DP  ER  ST  SA  PB  SQ 


LJ  DP  SA  ST  ER  SQ  PB 

Summer 


Winter 

2003-2004 


LJ 


SQ 


LJ  DP  ER  ST  SA  PB  SQ 

Gravid 

individuals 
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Table  4.4 


Mean  slope  estimates  from  the  reduced  major  axis  regression  of  fertility  at 
size.  Confidence  intervals  were  obtained  by  bootstrapping.  The  99.6% 
confidence  limits  ( a  =  0.004)  were  used  for  pairwise  comparisons  with  a 
group  significance  of  a=  0.05. 


Site 

Mean 

estimate 

95%  Conf.  Limits 

99.6%  Conf.  Limits 

Lower 

Upper 

Lower 

Upper 

Dana  Point 

2.949 

2.404 

3.549 

2.159 

3.888 

La  Jolla 

3.327 

2.855 

3.788 

2.529 

4.1 15 

San  Antonio 

2.833 

2.223 

3.299 

1.910 

3.623 

Santo  Tomas 

4.388 

3.857 

4.962 

3.548 

5.398 

San  Quinti'n 

3.640 

2.729 

5.261 

1.618 

7.159 

Punta  Baja 

3.095 

2.621 

3.552 

2.368 

3.801 
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Table  4.5 


Parameterization  and  asymptotic  analysis  of  projection  matrices  for  Punta  Baja  (PB),  La 
Jolla  -  San  Antonio  (LJ-SA)  and  Dana  Point  (DP).  Numbers  in  brackets  correspond  to 
95%  confidence  limits  computed  from  bootstrap  (#)  and  Monte  Carlo  (+)  distributions 
(see  Materials  and  Methods). 


Parameter 

Symbol/Formula 

PB 

SITE 

LJ-SA 

DP 

Per  capita  fertility 

/ 

2436 

6625 

6232 

[2032  -  2998]* 

[5844  -7530]* 

[5513  -  7025]* 

Larval  survival 

s* 

0.015 

0.015 

0.015 

Juvenile  survival 

a 

0.135 

0.044 

0.110 

Mar-Dec  Adult  survival 

b 

0.154 

0.860 

0.736 

Dec-Mar  Adult  survival 

d 

0.535 

0.951 

0.903 

Growth  rate 

afs  +  bd 

6.5 

18.1 

17.2 

(M 

[1.1  -  17.1]* 

[3.5  -  46.5]+ 

[3.2-44.2]* 

Net  Reproductive  Rate 

afs/(  1  -bd) 

7.0 

95.2 

49.0 

(*o) 

[1.1  -  18.6]+ 

[14.8  -  255.1]* 

[7.7  -  130.6]* 

Juvenile  life  expectancy 

\+[ad/(\-bd )] 

1.05 

1.47 

1.24 

(JLE) 

[1.03  —  1  08]* 

[1.24  1.72]* 

[1.12  -  1.37]* 

Adult  life  expectancy 

1/(1  -bd) 

1.09 

5.50 

2.97 

(ALE) 

[1.06  -  1.13]* 

[4.93  -  6.28]* 

[2.69  -  3.35]* 

*  :  Based  on  a  21 -day  larval  duration  and  a  constant  mortality  rate  of  0.2  d'1 
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Table  4.6 


Sensitivity  (5)  of  the  population  growth  rate  (A.)  and  net  reproductive  rate  (7?0)  to 
each  of  the  model  parameters.  Parameter  estimates  obtained  for  each  site  (Table 
5)  were  used  in  the  computations. 


Model 

parameter 

Sensitivity 

formula 

Punta  Baja 

La  Jolla 

San  Antonio 

Dana  Point 

5 

% 

5 

% 

s 

% 

Population  growth  rate  (X) 

/ 

as 

0.002 

0.00 

0.00! 

0.00 

0.002 

0.00 

s 

af 

328.86 

89.83 

291.50 

74.23 

685.52 

87.82 

a 

fi 

36.54 

9.98 

99.38 

25.31 

93.48 

11.97 

b 

d 

0.54 

0.15 

0.95 

0.24 

0.90 

0.12 

d 

b 

0.15 

0.04 

0.86 

0.22 

0.74 

0.09 

Net  reproductive  rate  (Rn) 

/ 

as/(\-bd) 

0.002 

0.00 

0.004 

0.00 

0.005 

0.00 

s 

afl(  1  - bd ) 

358.39 

89.10 

1600.42 

67.11 

2043.94 

82.67 

a 

fs/(\-bd) 

39.82 

9.90 

545.60 

22.88 

278.72 

1 1.27 

b 

afsdi(l-bd)2 

3.13 

0.78 

125.34 

5.26 

82.55 

3.34 

d 

afsbl(l-bd)2 

0.90 

0.22 

113.35 

4.75 

67.28 

2.72 
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4.9.  Figures 
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Figure  4.1 

Map  of  the  Southern  California  and  Baja  California  region  where  this 
study  was  conducted.  Labels  correspond  to  intertidal  sampling  sites. 
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Figure  4.2 

Life  cycle  diagram  of  the  stage-structured  matrix  model  for  the  intertidal 
bamacle  Balanus  glandula.  Arrows  correspond  to  transitions  in  the  model. 
Survival  probabilities  are  indicated  by  a,  b,  and  d.  The  product  fs  is  a 
reproductive  term,  proportional  to  per  capita  fertility  (/)  and  the  survival 
probability  for  larvae  (5). 
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Figure  4.3 

Size  distribution  of  Balanus  glandula  per  season,  site  and  reproductive  state  in 
Southern  and  Baja  California.  Size  is  the  first  principal  component  of  a  PCA 
performed  on  the  global  set  of  size  measurements  (Table  2).  Vertical  solid  and 
dashed  lines  are  the  mean  and  median  size  for  each  distribution.  Black  bars 
correspond  to  gravid  individuals.  Two-letter  codes  for  site  names  (see  Figure  1) 
correspond  to  Dana  Point  (DP),  La  Jolla  (LJ),  San  Antonio  (SA),  Santo  Tomas 
(ST),  Erendira  (ER),  San  Quintin  (SQ),  and  Punta  Baja  (PB). 
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Figure  4.4 

Distribution  of  residuals  from  a  two-way  contingency  table  analysis  of  size 
frequency  distributions  at  the  extremes  of  the  sampling  region.  Standardized 
residuals  were  computed  as  where  Ol}  corresponds  to  the  frequency 

of  size  interval  j  observed  at  site  i,  and  Ey  is  the  frequency  expected  for  that  site 
and  size  interval  under  the  hypothesis  of  independence. 
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Figure  4.5 


Time  series  of  the  proportion  of  gravid  individuals  at  each  site  and  sampling  date. 
Gray  lines  correspond  to  water  temperatures  measured  in  the  lower  intertidal  at 
each  site.  Data  for  La  Jolla  were  obtained  from  the  Network  of  Environmental 
Observations  of  the  Coastal  Ocean  database  (http://www.es.ucsc.edu/~neoco/). 
Dashed  vertical  lines  separate  the  seasons  used  in  the  analysis  of  size  data. 
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Figure  4.6 

Growth  parameters  for  Balcinus  glandula  at  Dana  Point  and  Punta  Baja.  Lines 
are  von  Bertalanffy  growth  curves  fitted  to  age-size  observations.  Equations 
with  the  corresponding  estimates  of  von  Bertalanffy  parameters  are  shown  on 
the  lower  right  comer  of  each  panel. 
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Figure  4.7 


Survival  curves  estimated  for  juveniles  (gray  lines  and  symbols)  and  adults 
(black  lines  and  symbols)  from  Dana  Point  (circles),  La  Jolla  -  San  Antonio 
(squares),  and  Punta  Baja  (triangles). 


127 


log  Fertility  (#  eggs) 


4.5 

DP 

©o  / 

4.5 

LJ 

4.0 

®>L<r 

-  V*?  • 

4.0 

•  • 

- 

3.5 

3.5 

9^  9 

© 

© 

3.0 

• 

3.0 

2.5 

logN  =  2.95*/ogPC  +  0.39 

2.5 

logN  =  3.33*logPC  -  0.04 

0.9  1.0  1.1  1.2  1.3  1.4  0.9  1.0  1.1  1.2  1.3  1.4 


4.5 

4.5 

4.0 

SA 

o 

4.0 

ST  *  v$ 

•  •  yv*  ° 

3.5 

° 

3.5 

co/  • 

o  So  • 

s'  o 

©  0/  9  o 

%  X  0° 

3.0 

3.0 

O 

o  / 

2.5 

2.5 

• 

logN  =  2.Q3*logPC  +  0.49 

logN  =  4.39* logPC  -  1.41 

0.9  1.0  1.1  1.2  1.3  1.4  0.9  1.0  1.1  1.2  1.3  1.4 


4.5 

4.5 

SQ 

PB 

4.0 

4.0 

0 

•°  ^ 

° 

3.5 

•o  *y 

3.5 

o°  C >y™°9o 

y'  °  • 

oo  o  *  # 

3.0 

° /  ° 

3.0 

os'  *  ft  O 

2.5 

2.5 

• 

O 

logN  =  3.64*/ogPC  -  0.33 

logN  =  3  09 VogPC  +  0.07 

0.9  1.0  1.1  1.2  1.3  1.4  0.9  1.0  1.1  1.2  1.3  1.4 


log  Size  (PCI) 


Figure  4.8 

Regressions  of  fertility  versus  size  in  adult  Balanus  glandula  from  Dana  Point 
(DP),  La  Jolla  (LJ),  San  Antonio  (SA),  Santo  Tomas  (ST),  San  Quintin  (SQ), 
and  Punta  Baja  (PB).  Parameter  estimates  for  the  linear  equations  shown  on 
each  panel  were  obtained  by  reduced  major  axis  regression  (see  Materials  and 
Methods).  Size  data  used  in  the  regression  correspond  to  the  first  principal 
component  of  a  PCA  performed  using  morphometric  measurements  on  the  202 
individuals  analyzed  for  fertility  (see  Table  4.2). 
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Figure  4.9 

Bootstrap  distributions  (10,000  repetitions  each)  of  per  capita  fertility  computed 
for  Dana  Point  (DP),  La  Jolla  -  San  Antonio  (LJ-SA),  and  Punta  Baja  (PB). 
Vertical  lines  correspond  to  mean  estimates  obtained  for  each  site  from  their 
mean  size  of  gravid  individuals  (Fig.  4.3)  and  size  specific  fertility  (Fig.  4.8). 
Confidence  limits  computed  with  the  bias-corrected  percentile  method  (Manly 
1997)  are  given  on  Table  4.5. 
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Figure  4.10 

Probability  distributions  of  population  growth  rate  (A),  net  reproductive  rate  (B),  and 
juvenile  and  adult  life  expectancies  (C,  D)  computed  for  Dana  Point  (DP),  La  Jolla  - 
San  Antonio  (LJ-SA),  and  Punta  Baja  (PB).  Vertical  lines  correspond  to  the  median 
of  each  distribution.  The  probability  of  finding  a  net  reproductive  rate  smaller  than 
one  (B)  were  calculated  as  the  number  of  cases  with  Ro<  1  divided  by  the  total  number 
of  repetitions.  Confidence  limits  computed  with  the  bias-corrected  percentile  method 
(Manly  1997)  are  given  on  Table  4.5. 
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Figure  4.11 

Fertility  thresholds  for  Balanus  glandula  at  Dana  Point  (DP),  La  Jolla  -  San 
Antonio  (LJ-SA),  and  Punta  Baja  (PB).  Contours  correspond  to  the 
minimum  per  capita  fertility  required  to  attain  population  persistence  given 
the  benthic  mortality  ratio  and  larval  survival  shown  on  the  x  and  y  axes. 
Gray  arrows  originate  from  the  x-axis  value  computed  for  each  site  (see 
Table  4.5),  intersect  their  corresponding  contour  of  mean  per  capita  fertility, 
and  point  to  the  minimum  larval  survival  (y-axis)  that  would  allow  them  to 
persist  (A,  >  1). 
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Figure  4.12 

Stochastic  growth  rate  log  Xs  computed  for  Dana  Point  (DP),  La  Jolla  -  San 
Antonio  (LJ-SA)  and  Punta  Baja  (PB)  under  different  regimes  of 
environmental  stochasticity.  Axes  correspond  to  lower  and  upper  limits  of  the 
uniform  probability  distributions  used  to  simulate  interannual  fluctuations  in 
larval  survival  (left  panels)  and  juvenile  survival  (right  panels).  The  thick 
contour  indicates  combinations  of  parameters  for  which  log  Xs=  0. 
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Figure  4.13 

Quasi-extinction  probabilities  Pg(0=O.l)  computed  for  Dana  Point  (DP),  La 
Jolla  -  San  Antonio  (LJ-SA)  and  Punta  Baja  (PB)  under  different  regimes  of 
environmental  stochasticity.  Axes  correspond  to  lower  and  upper  limits  of  the 
uniform  probability  distributions  used  to  simulate  interannual  fluctuations  in 
larval  survival  (left  panels)  and  juvenile  survival  (right  panels).  Regions  of  the 
parameter  space  for  which  Pq(Q)=l  are  shown  in  black. 


133 


5.  CHAPTER  5 


GENERAL  CONCLUSIONS 


The  results  of  this  dissertation  emphasize  the  importance  of  a  small-scale,  high- 
frequency  approach  to  study  nearshore  phenomena,  especially  in  relation  with  the 
dispersal  of  invertebrate  larvae.  The  importance  of  choosing  right  temporal  and  spatial 
scales  cannot  be  overstated.  In  recent  years  it  has  become  clear  that  often  there  is  not  a 
unique  right  scale  of  observation,  and  that  small-scale  physical  and  biological  processes 
in  the  nearshore  may  be  modulated  by  large-scale  phenomena  (e.g.  Pineda  and  Lopez 
2002).  The  question  of  whether  a  population  functions  as  an  open  or  closed  system  is  also 
dependent  on  the  spatial  and  temporal  scales  of  observation.  For  example,  populations 
separated  by  100  kilometers  may  not  exchange  individuals  at  demographically  relevant 
temporal  scales  (i.e.  years),  but  still  be  genetically  connected  by  larval  dispersal 
occurring  at  scales  of  decades  to  thousands  of  years. 

Observed  patterns  of  nearshore  larval  distribution  (Chapter  2)  agree  with  recent 
suggestions  that  long  distance  dispersal  of  larvae  in  the  nearshore  may  not  be  as 
pervasive  as  traditionally  thought  (e.g.  Jones  et  al.  1999;  Swearer  et  al.  1999).  Consistent 
small-scale  (0.1-1  km)  patterns  in  the  stage-specific  horizontal  distribution  of  Balanus 
glandula  and  Chthamalus  spp.  larvae  suggest  that  larvae  may  be  retained  within  short 
distances  to  shore  over  the  duration  of  their  planktonic  life.  Horizontal  segregation  of 
nauplii  and  cyprids  in  both  species,  with  the  latter  being  substantially  more  abundant  at 
the  inshore  station,  suggests  that  ontogenetic  changes  in  behavior  may  impact  the 
expected  horizontal  distribution  and  dispersal  distances  for  the  larvae  of  coastal 
invertebrates. 

The  mortality  estimates  obtained  for  barnacle  nauplii  in  nearshore  waters 
fluctuated  around  20-40%  per  day  (Chapter  2).  These  are  substantially  higher  than  a  5% 
d"1  previously  assumed  for  B.  glandula  and  Chthamalus  spp.  (Connolly  and  Roughgarden 
1998,  p.  323),  but  within  the  range  of  mortality  estimates  found  in  the  literature  for  other 
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benthic  invertebrate  larvae  (Rumrill  1990;  Morgan  1995;  Lamare  and  Barker  1999). 
Although  more  studies  are  needed  to  verity  these  estimates,  a  four-  to  six-fold  increase  in 
mortality  should  reduce  the  number  of  larvae  that  complete  their  planktonic  development 
and,  by  extension,  shorten  mean  travel  distances  and  scales  of  connectivity  among 
populations.  The  shape  of  a  dispersal  kernel  -  the  distribution  of  dispersal  distances  from 
a  point  source  -  depends  on  larval  duration  relative  to  local  hydrodynamics  (e.g.  Siegel  et 
al.  2003).  For  example,  dispersal  kernels  with  long  tails  may  be  observed  when  strong 
residual  currents  disperse  larvae  with  a  short  development  time.  Increasing  larval 
mortality  could  truncate  the  tails  of  such  distributions  and  cause  a  reduction  in  the  mean 
distance  of  dispersal.  A  recent  modeling  study  (Ellien  et  al.  2004)  has  shown  that  losses 
due  to  mortality  can  be  more  important  than  those  due  to  hydrodynamics  when  larval 
duration  is  short  relative  to  the  intensity  of  residual  currents.  Increasing  mortality  rates 
from  5%  to  9%  d'1  greatly  reduced  effective  dispersal  distances  for  larvae  of  the 
polychaete  Peclinaria  koreni.  Thus,  depending  on  the  species  and  local  hydrodynamics, 
changes  in  larval  mortality  rates  may  shift  our  view  as  to  whether  a  population  is  open  or 
closed. 


Results  presented  in  Chapter  2  suggest  that  coastal  invertebrate  populations 
separated  by  distances  as  short  as  tens  of  kilometers  might  be  demographically 
disconnected.  Such  disconnection  would  reduce  the  spatial  scales  of  coherence  in 
population  dynamics,  affect  the  robustness  of  local  populations  to  environmental 
perturbations,  and  ultimately  impact  the  persistence  of  a  species  within  its  geographic 
limits  of  distribution. 

The  observation  of  distinct  and  consistent  patterns  of  vertical  distribution  in  the 
nauplii  and  cyprids  of  two  barnacle  species  at  a  nearshore  site  (Chapter  3)  suggested  that 
horizontal  segregation  of  larval  stages  might  result  from  the  interaction  of  ontogenetic 
changes  in  behavior  and  vertical  shear  of  horizontal  currents  in  the  nearshore.  The 
absence  of  Chthamalus  spp.  cyprids  from  the  sampling  station,  located  ca.  2  km  offshore 
and  approximately  10  km  NW  of  La  Jolla  -  where  high  concentrations  of  Chthamalus 
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spp.  cyprids  have  been  observed  (Pineda  1999  and  Chapter  2)  -  is  intriguing  because 
Chthamalus  spp.  nauplii  were  numerically  dominant  and  present  at  all  stages  of 
development.  This  absence  did  not  allow  testing  whether  Chthamalus  spp.  nauplii  and 
cyprids  occur  at  different  depths.  Thus,  a  connection  between  the  horizontal  segregation 
observed  between  nauplii  and  cyprids  of  Chthamalus  spp.  at  La  Jolla  and  differences  in 
the  vertical  distribution  of  these  stages  must  be  tested.  More  field  observations  are  needed 
in  a  future  attempt  to  verify  the  biological-physical  mechanisms  that  generate  the 
observed  patterns  of  larval  distribution.  Future  efforts  to  model  the  relationship  between 
local  population  dynamics  and  hydrodynamics  must  shift  from  an  emphasis  on  mesoscale 
processes  to  a  better  description  of  nearshore  processes  and  spatial  changes  in  dispersal 
and  self-recruitment.  Also,  the  variability  in  nearshore  hydrodynamics  must  be  better 
resolved,  perhaps  with  high-frequency  observations  of  horizontal  flows  collected 
simultaneously  at  several  points  over  a  larger  area.  Estimates  of  dispersal  trajectories  and 
distances  would  benefit  from  an  assessment  of  the  spatial  variability  in  mean  horizontal 
flows  in  the  nearshore. 

Finally,  the  analysis  of  adult  demography  over  scales  of  hundreds  of  kilometers 
yielded  valuable  insights  on  the  interaction  between  spatial  changes  in  vital  rates, 
robustness  to  environmental  stochasticity,  susceptibility  to  local  extinction,  and 
geographic  ranges  of  distribution.  The  combination  of  extensive  field  observations  on 
population  structure,  growth,  and  reproductive  rates  with  a  simple  mathematical  model 
proved  to  be  a  fruitful  approach  to  the  assessment  of  factors  that  might  determine  the 
geographic  distribution  of  coastal  benthic  populations.  Future  studies  should  include  the 
actual  geographic  limit  of  distribution  of  this  species,  and  monitor  the  regional  variability 
in  recruitment  and  larval  supply,  as  well  as  the  variability  in  physical-chemical  factors 
that  are  likely  to  affect  rates  of  growth  and  reproduction  and  ultimately  generate  spatial 
patterns  in  population  demography. 
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Abstract 

We  investigated  the  effect  of  the  daily  sea  breeze  on  the  surface  flow  at  Bahia  Salsipuedes,  Baja  California,  Mexico. 
Drifter-tracking  experiments  were  conducted  over  a  2-week  period  in  July,  1999.  Surface  drifters  were  deployed  daily 
and  their  trajectories  tracked  from  a  small  boat  for  3-7  h.  Wind  speed  and  direction,  as  well  as  the  vertical  distribution 
of  temperature  and  current  velocities  were  measured  throughout  the  observation  period.  Our  results  indicate  that  strong 
winds  associated  with  the  daily  sea  breeze  at  the  study  area  can  produce  onshore  surface  currents  with  velocities  of  up  to 
10  cm  s”1  (net  transport  of  up  to  2  km  in  a  diurnal  cycle).  There  was  a  hyperbolic  relationship  between  onshore  wind 
stress  and  net  onshore  transport  of  the  drifters.  A  decrease  in  stratification  and  more  consistent  pattern  of  onshore 
transport  were  associated  with  sustained  winds  >6ms_l  towards  the  end  of  the  study.  Progressive  vector  diagrams 
computed  from  near-surface  currents  failed  to  predict  the  trajectory  of  drifters,  especially  on  days  where  wand  forcing 
was  more  consistent  in  terms  of  direction  and  magnitude.  We  discuss  the  implications  of  these  results  regarding  the 
collection  of  settlement  data  and  their  interpretation  in  terms  of  larval  supply  to  coastal  populations. 

©  2004  Elsevier  Ltd.  All  rights  reserved. 

Keywords:  Sea  breeze;  Onshore  transport;  Drifters;  Neuston;  Baja  California 


1.  Introduction 

The  dynamics  of  coastal  populations  is  highly 
dependent  on  the  mechanisms  and  processes 
affecting  the  spatial  distribution  and  the  fate  of 
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larvae  prior  to  settlement.  Research  efforts  during 
the  past  decades  have  been  focused  on  the  physical 
processes  that  are  likely  to  affect  the  distribution 
of  larvae  relative  to  adult  populations  (e.g. 
Johnson,  1939;  Scheltema,  1968;  Butman,  1987; 
Pineda,  1994b;  Epifanio  and  Garvine,  2001),  and 
on  the  importance  of  substrate  availability  (e.g. 
Bertness  et  ah,  1992;  Pineda,  1994a).  Such  a 
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‘supply  side’  approach  (see  Underwood  and  Fair- 
weather,  1989)  has  yielded  valuable  insights  on  the 
relevance  of  physical  processes  as  predictors  of  the 
local  and  regional  dynamics  of  adult  populations. 
In  order  to  characterize  the  variability  of  larval 
supply  to  coastal  populations,  however,  it  is 
necessary  to  assess  patterns  of  larval  distribution 
at  spatial  and  temporal  scales  that  are  relevant  to  a 
species’  life  cycle  and  ecology.  Plankton  patchiness 
and  logistic  constraints  preclude  the  use  of  such 
sampling  strategy  in  most  cases,  thus  forcing  the 
use  of  proxies  for  larval  supply  that  are  more 
feasible  to  monitor  over  extended  periods  of  time 
(e.g.  settlement — see  Pineda,  2000). 

Settlement  has  been  the  proxy  of  choice  in  a 
number  of  studies  focused  on  the  supply  of  larvae 
to  adult  populations  of  coastal  invertebrates,  such 
as  crabs  (Wing  et  al.,  1995b;  Reyns  and  Sponaugle, 
1999),  sea  urchins  (Ebert  et  al.,  1994;  Wing  et  al., 
1995a),  and  barnacles  (Gaines  et  al.,  1985;  Pineda, 
1994a).  Artificial  and/or  natural  surfaces  are 
provided  for  the  larvae  to  settle  on,  and  monitored 
at  intervals  of  days  to  weeks.  Settlement  time  series 
data  are  then  analyzed  in  conjunction  with 
environmental  data  series  to  identify  temporal 
correlations  between  larval  supply  and  physical 
variability  (e.g.  Farrel  et  ah,  1991;  Wing  et  ah, 
1995a;  Pineda  and  Lopez,  2002).  A  number  of 
published  studies  have  documented  correlations 
between  larval  settlement  and  physical  processes 
occurring  at  a  range  of  temporal  scales:  from  intra- 
seasonal  changes  in  wind  forcing  (Bertness  et  ah, 
1996),  upwelling  intensity  (Wing  et  ah,  1995a),  and 
the  stratification  of  the  water  column  (Pineda  and 
Lopez,  2002),  to  interannual  changes  in  river 
discharge  (Gaines  and  Bertness,  1992).  The  im¬ 
portance  of  higher-frequency  physical  variability 
for  settlement,  however,  has  not  been  extensively 
studied  and  cannot  be  resolved  at  the  sampling 
frequencies  typically  used  in  settlement  surveys. 
Therefore,  inferences  made  so  far  about  larval 
transport  are  probably  biased  towards  lower- 
frequency  mechanisms  if  based  only  on  settlement 
data  (Pineda,  2000). 

Few  high-frequency  mechanisms  with  a  poten¬ 
tial  effect  on  the  onshore  transport  of  larvae  have 
been  documented  so  far.  Internal  tidal  bores,  a 
seasonal  and  somewhat  predictable  phenomenon 


that  occurs  in  coastal  environments  at  scales  of 
hours  to  days,  are  an  example  of  high-frequency 
forcing  likely  to  affect  the  onshore  transport  of 
both  planktonic  and  neustonic  larvae  (Pineda, 
1994b,  1999). 

The  diurnal  sea  breeze,  which  is  also  seasonal 
and  locally  predictable  in  its  occurrence,  represents 
another  highly  episodic  phenomenon  that  could 
potentially  affect  larval  settlement  (Kaplan  et  al., 
2003;  but  see  Epifanio  and  Garvine,  2001,  p.  54). 
The  sea  breeze  is  characterized  by  a  daily  regime  of 
variability,  generated  by  oscillations  in  the  land- 
sea  temperature  gradient,  which  in  some  regions 
can  explain  fluctuations  of  up  to  7 ms-1  in 
onshore  wind  speed  (Simpson,  1994).  Such  daily 
surges  in  onshore  winds  may  affect  the  cross-shore 
distribution  of  material  suspended  in  the  upper 
centimeters  of  the  water  column  (i.e.  neuston), 
including  the  larvae  of  a  number  of  coastal 
invertebrates.  Thus,  by  modulating  the  onshore 
flow  of  surface  waters,  the  sea  breeze  could  have 
an  effect  on  the  timing  and  variability  of  recruit¬ 
ment  in  coastal  populations  with  neustonic  larvae. 

In  this  contribution  we  investigate  the  effect  of 
the  sea  breeze  on  cross-shore  surface  currents  by 
means  of  drifter-tracking  experiments,  which  were 
conducted  over  a  period  of  2  weeks  at  Bahia 
Salsipuedes  in  northern  Baja  California.  Our 
results  demonstrate  the  potential  for  the  diurnal 
sea  breeze  to  substantially  affect  onshore  transport 
of  neustonic  material.  We  discuss  the  implications 
of  these  results  regarding  the  collection  of  settle¬ 
ment  data  and  their  interpretation  in  terms  of 
larval  supply  to  coastal  populations. 

2.  Materials  and  methods 

The  experiments  were  conducted  on  13-22  July, 
1999  (Table  1),  at  Bahia  Salsipuedes  (31°  56'  N, 
116°  46'  W)  in  northern  Baja  California  (Fig.  I). 
Every  morning  at  about  7:30-9:30  AM  (PDT),  we 
deployed  a  set  of  5-6  drifters,  which  were  circular 
pieces  of  soft  and  flexible,  Scootguard™  non-slip 
surface  (30  cm  in  diameter,  1mm  thickness)  with 
six  small  lead  “split-shot”  lead  weights  (2.5  g  each) 
attached  to  the  edges  (Fig.  2).  When  deployed,  the 
drifters  remained  at  the  air-sea  interface,  with  less 
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Table  1 

Details  on  the  time  (PDT),  number  of  observations,  and  wind  conditions  during  the  drifter  deployments  at  Bahia  Salsipuedes,  Baja  California 
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Fig.  1.  Geographic  location  of  Bahia  Salsipuedes.  The  star 
indicates  the  study  site. 


than  50%  of  their  surface  area  exposed  to  direct 
wind  forcing  (Fig.  2).  Drifters  were  deployed  at  the 
onset  of  the  sea  breeze  events  (i.e.  a  sudden 
increase  in  onshore  winds  that  occurred  at  8:00- 
12:00  each  day),  at  approximately  the  same 
location  every  day  (31°  56.3  F  N;  116°  46.72'  W), 
and  tracked  from  a  small  motor  boat  for  3-7  h, 
depending  on  the  weather  conditions  and  on  the 
boat’s  proximity  to  the  surf  zone.  The  location  of 
each  drifter  was  recorded  about  every  20  min  using 
a  differential  GPS. 

We  assumed  that  the  drifters  were  tracking  the 
motion  of  the  uppermost  layers  of  the  water 
column,  rather  than  their  own  response  to  wind 
forcing.  Our  assumption  was  based  on  the  drifter 
design  (Fig.  2),  and  on  the  statistical  similarity 
between  the  drifter  paths  and  those  followed 
by  pieces  of  kelp  that  were  found  at  the  air-sea 
interface,  collected  in  situ,  and  deployed 
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Fig.  2.  Side  view  picture  (A),  and  top  and  side  view  schematics  (B  and  C)  of  the  drifters  used  in  this  study.  The  arrows  (and  black 
circles  in  B  and  C)  indicate  the  position  of  the  six  small  lead  weights  attached  to  the  drifter’s  edge.  The  edge  areas  (dark  shading  in  B 
and  C)  correspond  to  the  fraction  of  the  drifter’s  surface  that  remained  underwater.  The  central  area  (light  shade  in  B  and  C) 
corresponds  to  the  portion  that  was  exposed  to  the  wind  during  a  deployment. 


concurrently  with  the  drifters  (Solow,  WHOI, 
pers.  comm.). 

Throughout  the  observation  period  we  moni¬ 
tored  the  thermal  structure  of  the  water  column 
using  a  string  of  Onset  StowAway  XTI  loggers 
(response  time  <15  s),  programmed  with  a  sam¬ 
pling  interval  of  40  s  and  located  at  2,  5,  8,  and  1 1 
meters  above  the  bottom  (mab),  plus  one  logger  at 
the  surface.  We  also  monitored  the  vertical 
structure  of  nearshore  currents  with  a  bottom- 
mounted,  upward  looking  1200  kHz  Doppler 
current  meter  (RD  Instruments,  San  Diego,  CA), 
which  measured  at  intervals  of  40  s  (100  pings  per 
ensemble)  and  in  bins  of  0.4  m  within  a  depth 
range  of  0.95-14.55  mab.  Both  the  temperature 
mooring  and  the  current  meter  were  located  at  a 
depth  of  16m  (31°  56.28'  N;  116°  45.68'  W). 
Average  wind  velocity  and  direction  were  recorded 
over  5  min  intervals  using  a  CMB3280  Aanderaa 


Coastal  Monitoring  Buoy,  anchored  at  a  nearby 
location  (31°  55.28'  N;  116°  45.39'  W,  ca.  2  km 
south  of  the  current  meter  mooring).  Prior  to  the 
analysis,  wind  directions  were  rotated  and  aligned 
to  the  average  shoreline  orientation  (355°),  and 
then  decomposed  into  alongshore  and  cross-shore 
components. 

We  studied  the  effect  of  wind  forcing  on  the 
water  column  using  two  different  approaches. 
First,  on  a  time-domain  approach,  we  studied  the 
effect  of  wind  forcing  on  stratification  by  comput¬ 
ing  wind  stress  at  the  surface  as 

M  =  pu*2,  (l) 

where  p  is  air  density  and  u*  is  the  friction 
velocity,  which  was  computed  from  wind  speed 
data  as  in  Large  et  al.  (1995,  Eq.  8).  We  calculated 
stratification  as  the  square  of  the  Brunt-Vaisala 
frequency  (A2)  for  each  one  of  the  four  pairs  of 
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temperature  series  (i.e.  2-5  mab,  5-8  mab,  8- 
1 1  mab,  and  1 1  mab  -surface)  as 


where  g  is  gravity  (9.8  ms-2),  p0  is  the  mean 
density  of  the  water  column,  and  z  is  depth.  We 
then  obtained  depth-averaged  values  for  N2  by 
vertically  averaging  the  four  series  mentioned 
above.  The  temperature  and  wind  time  series  were 
low-pass  filtered  prior  to  the  computations  (cutoff 
frequency  of  6 cycles  h”1),  and  a  constant  salinity 
value  of  34psu  was  assumed.  Second,  on  a 
frequency-domam  approach,  we  performed 
spectral  analyses  on  the  time  series  for  wind 
(~  5,000  observations)  and  temperature 
(~  42,000  observations)  collected  at  the  study  site 
between  7-23  July  and  4—23  July,  respectively. 
Raw  data  were  linearly  detrended  prior  to  the 
computation  of  an  average  Welch  periodogram. 
The  wind  and  temperature  series  were  divided  into 
non-overlapping  segments  of  210  (3.6 d)  and  213 
points  (3.8  d),  respectively. 

To  understand  the  contribution  of  wind  to 
onshore  transport,  we  analyzed  the  effect  of  wind 
stress  on  the  motion  and  trajectory  of  the  drifters. 
We  computed  net  onshore  displacements  for  each 
drifter/day,  and  then  performed  a  nonlinear 
regression  of  net  onshore  displacement  versus  the 
onshore  wind  stress  integrated  over  the  duration  of 
each  deployment.  We  also  analyzed  changes  in  the 
correlation  between  drifter  trajectories  and  wind 
direction  as  a  function  of  wand  speed  for  each 
deployment.  For  n  average  drifter  positions 
recorded  during  a  deployment,  we  computed  n~\ 
velocity  vectors.  Velocity  vectors  were  decom¬ 
posed  into  cross-shore  and  alongshore  compo¬ 
nents  of  the  drifter  trajectories,  which  were  linearly 
regressed  against  the  respective  components  of 
wind  forcing.  We  chose  a  10  min  lag  between  wind 
measurements  and  the  drifter  observations 
to  which  they  were  paired  (i.e.  midpoint  between 
drifter  observations).  We  tested  whether 
drifter  trajectories  can  be  predicted  from  current- 
meter  data  by  comparing  the  observed  trajectories 
with  progressive  vector  diagrams  (PVDs)  from 
our  measurements  of  horizontal  velocities  near  the 
surface.  The  predicted  (*,y)  positions  were 


computed  as 

n  n 

x(tn)  =  A *53  “(O.-K'n)  =  A/  Y)  v0i),  (3) 

/=1  /'= 1 

where  «(/)  and  v(t)  are  the  east  and  north  current 
speeds  measured  by  the  current  meter  at  time  /, 
and  At  is  the  sampling  interval  (Emery  and 
Thomson,  1998,  p.  165).  Due  to  large  errors 
introduced  by  tidal  fluctuations  in  water  depth, 
surface  wave  action,  and  side  lobbing,  data 
collected  from  the  surface  bins  by  the  Doppler 
current  meter  had  to  be  discarded.  Thus,  our 
uppermost  bin  with  reliable  data  was  located  at 
13  mab,  on  average  2  m  below  the  surface. 

Finally,  we  observed  settlement  of  invertebrate 
larvae  on  the  drifters.  Although  these  observations 
were  not  quantitative,  we  use  the  taxonomic 
composition  of  the  settlers  and  the  time  scales 
over  which  settlement  was  observed  to  conjecture 
on  the  potential  effects  of  sea  breeze  events  on  the 
distribution  and  onshore  transport  of  neustonic 
larvae. 


3.  Results 

We  completed  7  days  of  drifter  observations 
during  which  a  total  of  36  drifters  were  deployed 
and  tracked  for  an  average  of  5  h  (Table  1).  Wind 
speeds  measured  during  the  study  were  highly 
variable,  with  most  of  the  variability  explained  by 
changes  in  the  cross-shore  component  (Figs.  3a 
and  4a).  Onshore  winds  dominated  during  the 
second  half  of  the  study,  with  daily  oscillations  of 
up  to  8ms-1  which  extended  for  about  8  h  (Fig. 
3a).  The  dominance  of  a  diurnal  frequency  in  the 
variability  of  cross-shore  winds  was  confirmed  by 
the  computed  power  spectra  (Fig.  4a).  Coincident 
with  these  daily  surges  in  onshore  winds,  we 
observed  short-term  increases  in  the  temperature 
of  surface  waters  (Fig.  3b).  Towards  the  end  of  the 
observation  period,  these  rapid  increases  in  tem¬ 
perature  could  be  detected  throughout  the  water 
column  (Fig.  3b),  which  is  consistent  with  the 
observed  decrease  in  the  depth-averaged  stratifica¬ 
tion  (Fig.  3c).  A  diurnal  frequency  dominated  the 
variability  of  temperature  between  the  surface  and 
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Fig.  3.  Wind  speed  (A),  water  temperature  (B)  and  mean  stratification  of  the  water  column  (C)  during  the  drifter  deployments  at  Bahia 
Salsipuedes  in  July  1999.  Tidal  heights  predicted  for  the  port  of  Ensenada,  Baja  California,  are  also  shown  in  C  (gray  line). 


I  I  mab  (Fig.  4b-c),  and  could  be  detected  at  depth 
to  2  mab  (Fig.  4f).  A  second  peak  in  the  power 
spectra  was  found  at  the  semi-diurnal  frequency 
(i.e.  about  2  cycles  day-1,  Fig.  4c-f)>  which  domi¬ 
nated  the  variability  of  temperature  near  the 
bottom,  but  nonetheless  could  be  detected  up  to 

I I  mab  (Fig.  4c-f).  Large  fluctuations  in  stratifica¬ 
tion  observed  at  the  beginning  of  our  study  (Fig. 
3c)  appeared  to  be  driven  by  the  semi-diurnal 
internal  tides  rather  than  wind  forcing. 

Persistent  wind  forcing  towards  the  end  of  the 
study  caused  a  net  onshore  transport  of  the 
drifters  (Fig.  5a),  which  increased  as  a  function 
of  wind  stress  integrated  over  the  duration  of  the 
deployments  (Fig.  5b).  Drifters  were  transported 
onshore  at  speeds  of  up  to  lOcms-1.  The  intensity 


of  wind  forcing  during  the  experiments  determined 
the  degree  of  correlation  between  drifter  paths  and 
wind  direction  (Fig.  6).  A  better  alignment  of 
drifter  paths  and  wind  on  the  cross-shore  direction 
is  consistent  with  the  higher  wind  speeds  recorded 
along  this  axis  (Fig.  3a).  The  response  of  surface 
waters  (and  therefore  drifter  paths)  to  the  tempor¬ 
al  variability  in  wind  forcing  can  be  clearly  seen  in 
Fig.  7,  in  which  the  wind  vectors  have  been 
superimposed  on  the  drifter  trajectories  for  July  20 
and  July  22. 

The  temporal  variability  in  the  horizontal 
velocities  of  near-surface  currents  (Fig.  8b)  was 
consistent  with  the  timing  and  intensity  of  wind 
forcing  (Fig.  8a).  Near-surface  currents  showed  a 
tendency  to  be  decoupled  from  currents  measured 
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Fig.  4.  Power  spectra  computed  for  the  wind  (A)  and 
temperature  time  series  (B-F)  recorded  at  Bahia  Salsipuedes 
in  July  1999.  Bars  on  the  right-hand  side  of  each  panel 
correspond  to  95%  confidence  intervals.  The  dashed  lines 
indicate  diurnal  variability;  mab:  meters  above  the  bottom. 


deeper  in  the  water  column  (Fig.  8b).  PVDs 
computed  from  horizontal  velocities  measured  at 
13  mab  (Fig.  9)  showed  a  poor  performance  as 
predictors  of  drifter  trajectories.  The  observed 
trajectories  seemed  to  be  the  result  of  an  interac¬ 
tion  between  wind  and  tidal  forcing,  which 
explained  most  of  the  variability  in  the  PVDs 
(Fig.  9).  The  discrepancy  between  predicted  and 
observed  drifter  displacement  varied  as  a  function 
of  the  magnitude  and  variability  of  wind  forcing 
during  the  deployments.  The  largest  discrepancies 
were  observed  on  days  with  strong  and  persistent 
winds,  whereas  smaller  values  tended  to  occur  on 
days  with  more  variable  wind  forcing,  both  in 
terms  of  speed  and  direction  (Table  1).  A  closer 
inspection  of  Fig.  9  and  the  wind  data  series  (Fig. 
3a)  suggests  that  not  only  the  magnitude  and 
variability  of  wind  forcing,  but  also  its  timing 
relative  to  the  tidal  cycle  (Fig.  3c)  ultimately 
determined  drifter  trajectories. 

We  incidentally  observed  settlement  of  neus- 
tonic  zooplankton  on  the  drifters.  The  most 
abundant  taxa  corresponded  to  flatworms  (0.5— 
52.8  indiv  h"1),  brachyuran  and  anomuran  mega- 
lopae  (0. 4-3.0  indiv  h“!),  isopods  (0. 2-2.0  in¬ 
div  h”1),  amphipods  (0.2-1 .7  indiv  h”1),  and 
Iepadid  barnacle  larvae  (0.4-1 .0  indiv  h  *).  The 
highest  settlement  rates  for  these  taxa  were 
observed  during  the  last  four  days  of  the  study 
(19-22  July),  which  exhibited  the  strongest  on¬ 
shore  winds  (Fig.  3a). 

4.  Discussion 

Our  results  indicate  that  neustonic  material  has 
the  potential  to  be  transported  onshore  by  the 
daily  sea  breeze  at  Bahia  Salsipuedes  (Fig.  5).  A 
causal  relationship  is  supported  by  the  observed 
drifter  trajectories  (Figs.  5  and  6)  and  their 
response  to  short-term  fluctuations  in  wind  direc¬ 
tion  (Fig.  7),  as  well  as  by  the  discrepancies 
between  actual  drifter  trajectories  and  those 
predicted  by  progressive  vector  diagrams  (PVDs) 
computed  from  the  horizontal  velocities  of  near- 
surface  currents  (Fig.  9).  These  results  are  relevant 
because  currents  in  the  coastal  ocean  are  generally 
more  energetic  on  the  alongshore  axis  (Winant, 
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Polal  onshore  wind  stress  (N  m*2) 


Fig.  5.  (A)  Average  paths  followed  by  drifters  on  different  days;  the  numbers  are  the  July  date  corresponding  to  each  trajectory, 
whereas  the  arrow  indicates  the  point  of  deployment.  (B)  Net  onshore  displacement  of  drifters  as  a  function  of  onshore  wind  stress 
integrated  over  the  duration  of  each  deployment.  The  large  open  circle  in  (A)  indicates  the  position  of  the  Doppler  current  meter  and 
temperature  moorings. 


Th 
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Fig.  6.  (A)  Changes  in  the  correlation  between  average  drifter  and  wind  direction  as  a  function  of  the  average  wind  speed  during  each 
deployment.  Linear  regression  analyses  were  performed  separately  for  alongshore  (closed  symbols)  and  cross-shore  (open  symbols) 
components  of  wind  and  drifter  direction,  and  for  each  day  with  drifter  deployments  (July  13,  15,  16,  19,  20,  21,  22).  (B  and  C) 
Example  of  how  the  values  indicated  with  squares  (in  A)  were  obtained  for  July  21;  W  is  average  wind  speed  for  the  deployment. 


1980),  with  few  mechanisms  acting  on  the  cross¬ 
shore  dimension  with  sufficient  energy  to  generate 
onshore  transport  at  the  temporal  and  spatial 
scales  observed  here  (e.g.  internal  tidal  bores — see 
Pineda,  1994b). 

The  surface  drifters  used  in  this  study  are  an 
effective  means  to  investigate  the  short-term 


response  of  the  surface  layer  (i.e.  first  few 
centimeters)  to  local  wind  forcing.  Water  move¬ 
ment  at  this  biologically  relevant  layer— where  a 
suite  of  particles  including  neustonic  plankton  can 
be  accumulated  (see  LeFevre,  1986)— is  left 
unresolved  by  most  instruments,  including  acous¬ 
tic  current  meters  and  other  drifter  designs  that 
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Fig.  7.  Example  of  drifter  trajectories  in  relation  to  the 
temporal  variability  of  wind  forcing  on  July  20  and  July  22, 
1999.  Open  circles  correspond  to  drifter  positions,  whereas  the 
thin  lines  are  a  vectorial  representation  of  the  average  wind 
speed  and  direction  measured  during  the  interval  prior  to  each 
drifter  observation.  Closed  circles  indicate  the  point  of  drifter 
deployment. 


can  go  deeper  into  the  water  column  and/or  be 
directly  affected  by  wind  forcing.  One  disadvan¬ 
tage  of  our  approach,  however,  is  the  need  to 
visually  and  continuously  track  the  drifters  from  a 
boat.  Currently  available  technology  should  make 
it  feasible  to  address  such  problem  by  attaching  a 
small  positioning/tracking  device  to  each  surface 
drifter,  so  that  their  position  can  be  recorded 
without  compromising  their  ability  to  track  the 
surface  layer. 

The  discrepancies  we  observed  between  actual 
drifter  paths  and  those  predicted  by  PVDs  should 
be  taken  into  account  before  any  conclusions 
about  larval  transport  are  drawn  from  Eulerian 
current  data.  Our  data  suggest  that  the  observed 
discrepancies  between  the  PVDs  and  drifter 
trajectories  may  be  due  to  an  interaction  between 
the  magnitude  and  variability  of  wind  forcing  and 
tidal  currents.  However,  one  source  of  error  that 
could  explain  part  of  the  discrepancy  is  the  varying 
distance  between  our  top  bin  with  reliable  current 
meter  data  (13mab)  and  the  surface  of  the  water 
column.  Given  that  our  main  objective  was  to 
study  the  effect  of  the  daily  sea  breeze  on  the 
onshore  transport  of  surface  currents,  the  experi¬ 


ments  were  conducted  at  roughly  the  same  time  of 
the  day.  This  means  that  on  some  days  the 
observations  took  place  at  high  tide  (July  13,  15), 
at  low  tide  (July  20,  21),  or  during  the  transition 
from  low  to  high  (July  16,  19).  Unfortunately,  a 
pressure  sensor  was  not  available  and  therefore  we 
lack  accurate  data  on  the  sea  level  at  our  study  site 
during  drifter  deployments.  Tidal  predictions  for 
Ensenada  (ca.  20  km  south  of  Bahia  Salsipuedes), 
indicate  sea-level  fluctuations  of  up  to  2.4  m  during 
the  course  of  the  deployments,  which  would  have 
increased  the  distance  from  the  top  ADCP  bin  to 
the  surface  from  2  to  4.4m.  Such  a  change  could 
certainly  increase  the  discrepancy  between  the 
trajectories  predicted  from  a  PVD  and  those 
actually  followed  by  the  drifters,  especially  on  a 
highly  stratified  water  column  (Figs.  8  and  9). 

Our  observations  were  focused  on  drifter  paths 
and  on  the  physical  variables  associated  with  the 
water  column’s  response  to  wind  forcing,  and  did 
not  include  an  assessment  of  the  spatial  distribu¬ 
tion  and/or  settlement  of  neustonic  larvae  in  the 
area.  Therefore,  a  connection  between  this  trans¬ 
port  mechanism  and  the  recruitment  of  coastal 
species  with  neustonic  larvae  is  yet  to  be  demon¬ 
strated.  If  there  is  a  connection,  however,  and 
larvae  of  coastal  invertebrates  are  in  fact  trans¬ 
ported  onshore  by  the  sea  breeze,  it  would  be 
possible  to  predict  the  spatial-temporal  distribu¬ 
tion  of  larval  settlement  in  some  populations  based 
on  local  wind  patterns  and  shore  orientation.  An 
assessment  of  the  effect  of  sea  breeze  events  on  the 
spatial  distribution  and  timing  of  settlement  is  thus 
warranted. 

A  follow  up  of  our  study  should  include  a 
survey  of  neustonic  meroplankton  in  the  area  with 
concurrent  drifter  deployments  and  settlement 
observations.  Given  the  diurnal  frequency  of  sea 
breeze  events,  any  attempt  to  link  such  events  to 
local  settlement  patterns  will  have  to  consider 
observations  taken  at  shorter  intervals  than  it  is 
customary  (i.e.  at  least  semi-diurnally).  Lack  of 
temporal  resolution  could  result  in  a  misinterpre¬ 
tation  of  the  mechanisms  that  determine  the 
timing  and  intensity  of  settlement  in  populations 
with  neustonic  larvae. 

Wind-driven  onshore  transport  of  neustonic 
invertebrates  has  been  suggested  as  a  cause  of 
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Fig.  8.  Example  of  the  effect  of  wind  forcing  (A)  on  the  circulation  (B)  and  thermal  structure  (C)  of  the  water  column  in  Bahia 
Salsipuedes  on  July  1 9,  1 999.  The  shaded  area  in  A  corresponds  to  the  duration  of  the  drifter  deployment.  The  stickplols  in  B  show  the 
variability  in  direction  and  magnitude  (see  scale  bar  at  lower  left  corner)  of  horizontal  currents  at  the  same  depths  at  which  temperature 
was  recorded.  The  vectors  indicate  the  direction  currents  are  flowing  towards  (up  is  north). 


the  spatial  patterns  observed  in  barnacle  settle¬ 
ment  (Kendall  et  al.,  1982;  Bertness  et  al.,  1996), 
and  in  the  spatial  distribution  of  neustonic  isopods 
(Locke  and  Corey,  1989),  post-larval  crabs 


(Shanks,  1986;  Eggleston  et  al.,  1998)  and  lobster 
megalopae  (Hudon  and  Fradette,  1993).  If  such 
wind  events  are  persistent  and  consistent  enough 
over  demographically  relevant  time  scales,  they 
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Fig.  9,  Comparison  of  observed  versus  predicted  trajectories  of 
drifters  at  Bahia  Salsipuedes,  Baja  California,  in  July  1999 
Open  circles  indicate  the  observed  drifter  trajectories.  Predicted 
trajectories  (closed  circles)  correspond  to  a  PVD,  computed 
from  horizontal  current  velocity  data  collected  at  I3mab.  The 
gray  line  is  a  PVD  computed  using  only  those  frequencies  of  the 
current  data  that  correspond  to  the  barotropic  tide.  The  insets 
show  average  wind  speed  and  direction  (0°  =  N)  over  the 
duration  of  each  deployment.  The  scale  bar  at  the  lower  left 
corner  of  each  panel  corresponds  to  a  distance  of  1  km. 


shoreline  throughout  their  development,  or  to 
repeatedly  reach  it,  hence  increasing  their  prob¬ 
ability  of  being  close  to  adult  habitats  when 
competent  to  settle.  Thus,  the  sea  breeze  could 
constitute  an  effective  transport  and/or  retention 
mechanism  in  areas  with  a  configuration  of 
shoreline  and  wind  forcing  such  as  that  observed 
at  Bahia  Salsipuedes. 

The  effect  of  wind  forcing  on  larval  retention 
and  demographic  isolation  could  be  intensified  if 
larval  behavior  is  invoked,  such  as  the  reverse 
vertical  migration  hypothesized  by  Poulin  et  al. 
(2002)  for  neustonic  larvae  of  the  subtidal  gastro¬ 
pod  Concholepcis  concholepas.  Even  under  condi¬ 
tions  of  strong  upwelling  (i.e.  strong  offshore 
Ekman  transport),  the  larvae  of  C.  concholepas  are 
consistently  found  between  the  shoreline  and  the 
upwelling  front.  Such  observation,  along  with  data 
on  the  temporal  variability  of  larval  distribution, 
suggests  that  a  reverse  vertical  migration  behavior 
might  be  allowing  these  larvae  to  stay  close  to 
subtidal  habitats  (see  Fig.  6  in  Poulin  ct  al.,  2002). 
The  interaction  between  physical  processes  like  the 
sea  breeze  and  behavioral  mechanisms  such  as 
vertical  migrations  could  increase  the  local  reten¬ 
tion  of  neustonic  larvae  and  the  demographic 
isolation  of  adjacent  populations. 

Finally,  these  results  might  also  be  relevant  for 
understanding  the  onshore  transport  of  pollutants, 
often  discharged  a  few  kilometers  from  shore  near 
coastal  cities  (e.g.  Boehm  et  ah,  2002).  Such 
discharges  are  likely  to  contain  buoyant  particles 
than  can  reach  surface  waters.  Our  observations 
suggest  that  such  particles  could  be  rapidly 
transported  onshore  by  the  sea  breeze. 


could  have  an  impact  on  the  connectivity  and/or 
demographic  isolation  among  adjacent  coastal 
populations.  At  sea  level,  the  sea  breeze  is  highly 
asymmetrical,  with  a  marked  predominance  of 
onshore  over  offshore  winds  (Fig.  3a).  Hence,  the 
sea  breeze  can  have  a  cumulative  effect  (i.e.  over 
more  than  one  diurnal  cycle)  on  the  cross-shore 
distribution  of  neustonic  larvae.  A  persistent  sea 
breeze-mediated  onshore  transport  could  allow 
neustonic  larvae  to  either  remain  close  to  the 
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